第 8 章 岩石流体力学.ppt

第8章岩石流体力学,8.1概论8.2岩石流体力学特性8.3岩石流体力学基础8.4岩体--液体耦合力学模型及解法8.5岩体--气体耦合力学模型及解法,8.1概述,岩石流体力学固流耦合）RockFluidmechanics[LouisC.1974提出岩石水力学,RockHydraulics主要研究固流耦合作用下岩体变形、破坏与流体渗流的规律及其工程效应。,第8章岩石流体力学,8.1概述,研究领域油气开采、煤层气开采、水文地质、地下水动力学、环境科学、水利水电工程、水库诱发地震、海水入侵、地下水引起的地面变形、矿井突水与瓦斯突出、水上水下采煤、常规地热开发、HDR开发、水压致裂等极为广泛的工程领域。特点通过考虑岩体与流体的相互作用（含物理与化学作用）给出工程实际性态的描述。而单纯的渗流力学或固体力学则是通过估计另一方的影响而给出描述，显然后两者存在较大误差，甚至会引起错误。,第8章岩石流体力学,8.1概述,社会实践推动了相关科学研究的发展1）石油、天然气、煤层气开采2）环境工程核废料、高污染化学废料3）地震科学4）水利、水电工程5）煤矿开采瓦斯问题，水突出问题6）地热开采常规地热、高温岩体地热（HDR）,第8章岩石流体力学,8.1概述,研究方法角度划分1）研究描述各类问题的耦合数学模型及其解法岩体与水、岩体与粘性流体、岩体与气体、固气液、固流热的耦合数学模型。2）研制新型实验设备，以揭示流体与岩体相互作用的内在规律有效应力、化学流体作用下岩体特性、岩体渗流特性等。3）研究岩石流体力学的新工艺或指导工程设计。,第8章岩石流体力学,8.2岩石流体力学特性8.2.1流体与骨架的性质,1流体的基本性质A密度定义为流体单位体积的质量.随压力p热力学温度T变化.状态方程p,T或f,p,T0单位物理单位制ML-3工程单位制FT2L-4真实气体p/Zp,TRT理想气体p/RT,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质重率定义为流体单位体积的重量。g单位物理单位制M,L,T,的量纲ML-2T-2工程单位制F,L,T的量纲FL-3比重或相对密度液体比重为4C时液体密度与纯水密度之比.气体规定温度下气体密度与氢气或空气密度之比.单位cm.g.s制中,密度单位g/cm3重率dyn/cm3,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质B黏度牛顿流体,第8章岩石流体力学,8.2.1流体与骨架的性质,非牛顿流体不服从方程1的流体都是非牛顿流体.所有气体及其简单液体都是非牛顿流体.ABingham塑性流体B假塑性流体C胀流性流体D触变性流体E流变性流体,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质黏度的单位动力黏度量纲ML-1T-1,或FL-2T厘米克秒制中1泊P1达因秒/厘米2克厘米/秒2*秒/厘米克厘米-1秒-11厘泊cP0.01P,20C水的动力黏度为1cP.国际标准单位帕斯卡.秒Pa.s,换算关系1P0.1Pa.s,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质黏度的单位运动黏度/量纲L2T-1,或FL-2T厘米克秒制中单位厘米2/秒cm2/s1cm2/s称为1斯托克斯,常用单位为厘斯托克斯.,第8章岩石流体力学,8.2.1流体与骨架的性质,1流体的基本性质流体的压缩系数压缩系数就是当物质承受法向压力或法向张力变化时,其体积和密度变化的度量.压缩系数的倒数称为弹性模量.,第8章岩石流体力学,8.2.1流体与骨架的性质,2孔隙岩体的基本性质孔隙率与有效孔隙率有效孔隙率定义为介质中相互连通的孔隙即有效孔隙的体积与介质总体积比值.neUve/Ub,UveUvneUv,第8章岩石流体力学,8.2.1流体与骨架的性质,2孔隙岩体的基本性质比面多孔介质的比面M定义为单位体积多孔介质Ub内所有颗粒的总表面积AsMAs/Ub,[M]L-1基本特征细粒的表面比粗粒的大的多.砂岩的比面为1500cm2/cm3,无烟煤最大为287m2/g.比面的测定方法统计法;吸附法;湿热法比面特征的模拟门杰海绵,席尔宾斯基地毯,第8章岩石流体力学,8.2.1流体与骨架的性质,第8章岩石流体力学,8.2.2岩体渗流物性方程,1线性渗流的物性方程Darcy定律QKAh1-h2/LqKJJh1-h2/LDarcy定律的推广,第8章岩石流体力学,8.2.2岩体渗流物性方程,1线性渗流的物性方程Darcy定律的适用范围雷诺数表示Reqd/层流区1--10低流速,粘滞力过渡区100惯性力起主要作用紊流区,第8章岩石流体力学,8.2.2岩体渗流物性方程,1线性渗流的物性方程水力传导系数与渗透率Darcy定律中的比例系数K称为水力传导系数或渗透系数。在各向同性介质中，水力传导系数定义为单位水力梯度的比流量，它是一个表示多孔介质输运流体能力的标量量纲L/T。特征与流体及骨架的性质有关。流体性质为密度，动力黏度或运动粘度；骨架性质为粒径分布，颗粒、形状、比表面、弯曲率、孔隙率等。水力传导系数Kk/kg/k（量纲L2）叫多孔骨架的渗透率或内在渗透率，仅与骨架性质有关,第8章岩石流体力学,8.2.2岩体渗流物性方程,1线性渗流的物性方程水力传导系数与渗透率单位在米制单位中渗透率单位cm2或m2，对于摄氏20度的水，K1cm/s相当于k1.02*10-5cm2工程中渗透率单位达西（d）公式k（Q/A）/△p/△x对于摄氏20度的水,可做如下换算1d9.869710-9CM29.61310-4CM/S许多情况下用md,1md10-3d,第8章岩石流体力学,8.2.2岩体渗流物性方程,2非线性渗流规律水Jones1975提出碳酸钙岩石裂隙渗透系数公式LouisPeuga1976提出应力对渗透系数的影响赵阳升等式中a,b,c为常数.p为平均孔隙压,MPa,K的单位取cm/s,或对于水可换算为md.,第8章岩石流体力学,8.2.2岩体渗流物性方程,2非线性渗流规律,第8章岩石流体力学,,2非线性渗流规律3D应力下,孔隙裂纹气体渗流规律,阳泉3煤层煤样瓦斯渗透系数与体积应力、孔隙压的关系图,第8章岩石流体力学,8.2.2岩体渗流物性方程,8.2.2岩体渗流物性方程,3裂隙介质中的渗流物性方程发展50年代苏联学者研究Pomm1966Snow1986Louis1974Witherspoon[1982]基本搞清了裂隙岩体的渗流规律,第8章岩石流体力学,8.2.2岩体渗流物性方程,3裂隙介质中的渗流物性方程1单一裂缝水力学模型等宽度裂隙模型对于壁面光滑的细小裂隙,运动方程为Louis1974qKfJf层流渗透系数为Kfgd2/12对于裂隙面粗糙不平Kfgd2/12C为连通系数,连通面积与总面积的比;C为裂隙面相对粗糙度系数C18.8△/2d1.5沟槽流模型Tsang1987Kfgd3/12,第8章岩石流体力学,8.2.2岩体渗流物性方程,3裂隙介质中的渗流物性方程裂隙岩体介质的渗透规律Romm1966假定裂隙岩体介质为拟连续介质体,水在整体岩体中的流动用Darcy定律的形式qkfd/LJ裂隙组的等效渗透系数Kegd3/12LLouis认为裂隙组的渗透系数为Kd/LKfKm层流紊流均适用,第8章岩石流体力学,8.2.2岩体渗流物性方程,3裂隙介质中的渗流物性方程三维应力下,裂缝的渗透规律,第8章岩石流体力学,3D应力下,裂缝水渗流规律,第8章岩石流体力学,3D应力下,裂缝气体渗流规律,第8章岩石流体力学,3D应力下,裂缝气液二相流体渗流规律,Sgc40.725,,第8章岩石流体力学,8.2.3有效应力规律,1太沙基原理1923年,太沙基在研究土的固结时提出Biot1940’s发展为三维有效应力规律价值1）饱和土中两个受力体系的相互关系；2）土体变形与两种应力的关系；3）土体强度与两种应力的关系。1950‘s年代,土压力三轴仪研制成功,Skempton实验发现孔隙压力不仅受法向应力作用,同时受剪应力影响,第8章岩石流体力学,8.2.3有效应力规律,2修正的太沙基原理1960年,Skempton证明太沙基原理对工程应用足够精确.1960年,Skempton发现太沙基有效应力原理对有些土工工程不适用,提出修正形,Robinson1959、Handin1963研究才形成如下形式01称为有效应力系数.岩石的有效应力是岩体不均匀分布的孔隙裂隙在孔隙压的作用下的宏观响应的平均表示.近年来,赵阳升证明有效应力系数是孔隙压和体积应力的双线性函数.,第8章岩石流体力学,有效应力系数规律赵阳升,1995年a1a2a3pa4p,第8章岩石流体力学,8.2.3有效应力规律,2修正的太沙基原理有效应力系数的分区特征,第8章岩石流体力学,8.2.3有效应力规律,2修正的太沙基原理对于裂缝而言法向应力为,第8章岩石流体力学,8.2.4流体作用下的岩体特性,1问题孔隙裂隙流体的存在对天然岩体有很强的破坏作用,其影响程度取决于岩体矿物组成与流体组成的化学特性及其组合作用.2水的影响Rca-bWc赵阳升,靳钟铭Ea/Wc-b李成江Ea-bp胡耀青,赵阳升3气体的影响日本氏平增之1986姚宇平1988,赵阳升贺军1991,梁冰19944热流体与化学流体的影响,第8章岩石流体力学,8.2.4流体作用下的岩体特性,第8章岩石流体力学,8.2.4流体作用下的岩体特性,第8章岩石流体力学,8.3岩石流体力学基础8.3.1不变形多孔介质中的质量守恒,1基本连续性方程近似地在P点按泰勒展开,并忽略高级项得控制体表面流入超出的数量质量守恒或,第8章岩石流体力学,8.3.1不变形多孔介质中的质量守恒,第8章岩石流体力学,8.3.1不变形多孔介质中的质量守恒,质量通量质量守恒方程对于不变形介质n常数将Darcy定律代入方程对于不可压缩流体Laplace方程,第8章岩石流体力学,8.3.1不变形多孔介质中的质量守恒,可压缩流体渗流控制方程理想气体ZP1介质均质各向同性,第8章岩石流体力学,8.3.2可压密介质中的质量守恒,第8章岩石流体力学,质量守恒方程固体骨架的可压缩性岩体骨架的变形由固体颗粒体积Us与孔隙体积Up的变形组成相应的压缩系数为UbUsUpUs1-nUbUpnUb一般故,8.3.2可压密介质中的质量守恒,第8章岩石流体力学,,只有垂向压密的问题岩体的压缩系数为，有效应力骨架的变形为流体的压缩系数为当const定义Sop为比质量储存系数对于均质各向同性体有,8.3.3承压含水层和越流含水层的流动,第8章岩石流体力学,,承压含水层忽略厚度的影响导水系数TKb量纲L2/T,m2/h含水层储水系数为含水层垂线上的平均测压水头下降或上升一个单位时,从单位水平面积含水层中释放或储存的水的体积比储水系数SSSb承压含水层渗流控制方程,8.3.3承压含水层和越流含水层的流动,第8章岩石流体力学,,8.3.3承压含水层和越流含水层的流动,第8章岩石流体力学,,越流含水层,8.3.3承压含水层和越流含水层的流动,第8章岩石流体力学,,越流含水层主要考虑沿弱透水层的垂向运动.设垂向比流量为q’v,则有对于均质各向同性的主含水层有对于弱透水层有当弱透水层上面存在一保持固定势的含水层式中,8.3.3承压含水层和越流含水层的流动,第8章岩石流体力学,,越流含水层引入称为越流因素的另一个越流参数,8.3.4流函数与势函数,第8章岩石流体力学,,流线如果瞬时曲线上每一点的切线方向也是流体在该点的速度方向,则此曲线为流线,qdr0另一形式q用势函数代替在二维情况下qydx-qxdy0解为x,yconst称为流函数.对于均质各向同性岩体q-grad,k,Cauchy--riemann条件,8.3.4流函数与势函数,第8章岩石流体力学,,对于均质各向同性体,等势线处处垂直于流线非均质各向同性体,8.3.5初边值条件,第8章岩石流体力学,,解的适定性存在;唯一;稳定边界条件一类边界给定势边界Dirichlet条件二类边界给定通量qnq*lnqnx,y,zNewmann条件三类边界混合边界Cauchy条件,8.3.5初边值条件,第8章岩石流体力学,,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,研究历史1960‘s开始,理论相对成熟.裂隙岩体的渗流理论分为两大类1裂隙网络渗流模型,核心思想裂隙和孔隙看成单独的两套系统,先单独研究,再考虑二者的水力交换.Barenblantt;Warren;Root,Wittke1969,19702拟连续介质模型,核心思想认为孔隙渗流很小,可忽略.考虑裂隙分布特征,确定渗透系数,建立理论模型.代表Romm1966,Snow1969,1972,Louis,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,1裂隙网络渗流模型Barebblalltt假设a重叠的孔隙裂隙连续体;随机分布的裂隙系统,存在质量交换;b根据量纲分析导出的质量交换量QpfCKfP1-P2/c假设裂隙岩体对于裂隙和孔隙二类介质仍然是均质各向同性的.控制方程,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,1裂隙网络渗流模型,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,1裂隙网络渗流模型Wittke模型假设裂隙岩体的渗流空间由裂隙个体组成的裂隙网络所组成.根据质量守恒,可以给出三组方程a在裂隙交点处,流入与流出节点的流量代数和为0.B分割各基质岩块的裂隙看成一多边形的闭合回路,而各各回路的压差的代数和为0.C单一的裂隙形成的弯曲通道,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,8.3.6裂隙岩体的渗流模型,第8章岩石流体力学,,2拟连续介质模型Romm模型假设裂隙均布于岩体中,孔隙很不发育.按等效的连续介质处理.关键是寻求等效渗透系数.对于一组裂隙,8.4岩体--液体耦合力学模型及解法,第8章岩石流体力学,,8.4.1连续介质模型基本假设1固体骨架为弹性均质各向同性介质.2水渗流在微段压力梯度上服从达西定律,整体上非线性分布;3介质为单相水所饱和;4服从修正的太沙基有效应力规律;5多孔介质的体积变形等于孔隙变形.控制方程,8.4.2块裂介质模型块裂介质固体变形与液体渗流耦合数学模型及解法基质岩块液体渗流方程裂缝液体渗流方程基质岩块变形方程裂缝变形方程,裂缝有效应力规律基质岩块有效应力规律,8.4岩体--液体耦合力学模型及解法,第8章岩石流体力学,8.5岩体--气体耦合力学模型及解法,第8章岩石流体力学,,8.5.1连续介质模型基本假设1固体骨架为弹性均质各向同性介质.2水渗流在微段压力梯度上服从达西定律,整体上非线性分布;3介质为单相气体所饱和;4服从修正的太沙基有效应力规律;5多孔介质的体积变形等于孔隙变形.6气体为理想气体;７气体含量服从拉格缪尔公式.控制方程,8.5岩体--气体耦合力学模型及解法,第8章岩石流体力学,,块裂介质固体变形与气体渗流耦合数学模型基质岩块气体渗流方程裂缝气体渗流方程岩块变形方程裂缝变形方程,岩块有效应力规律裂缝有效应力规律,8.5岩体--气体耦合力学模型及解法,第8章岩石流体力学,,8.5.3求解策略按时间序列线性近似的非线性方程迭代求解方法,,,块裂介质固气耦合的煤层气开采理论模型,第8章岩石流体力学,瓦斯抽放仿真图（3）抽放钻孔位于裂缝交会点,第8章岩石流体力学,