渗流力学历年习题集锦

渗流力学历年习题集锦渗流力学作为研究流体在多孔介质中运动规律的科学，是许多工程学科的理论基础，其应用广泛涉及石油与天然气开发、地下水利用、水利工程、环境工程等多个领域。掌握渗流力学，不仅需要深刻理解基本概念与原理，更离不开大量的习题练习。通过习题，可以将抽象的理论具体化，将复杂的公式应用于实际问题，从而真正内化知识，提升分析与解决问题的能力。本文旨在整理一份渗流力学历年习题集锦，这些习题均选自不同时期、不同层次的教学与考核，具有一定的代表性和典型性。我们将按照渗流力学的知识模块进行分类，并辅以简要的解析思路，希望能为各位同学和研究者提供一个系统的练习与回顾资料。请注意，习题的价值不仅在于答案本身，更在于解题过程中的思考与方法的锤炼。一、渗流基本概念与渗流基本定律说明：本部分主要考察对渗流力学核心概念的理解，以及达西定律的掌握程度。习题1：概念辨析与基本计算题目：简述渗流速度与真实渗流速度的区别与联系。若某多孔介质的孔隙度为n，渗流速度为v，试写出真实渗流速度u的表达式。在实验室中，通过一根长度为L、横截面积为A的岩心进行渗流实验，测得入口压力为p₁，出口压力为p₂（p₁>p₂），实验流体粘度为μ，单位时间内流出岩心的流体体积为Q。试根据达西定律计算该岩心的渗透率k，并说明达西定律的适用条件。解析思路：*首先明确渗流速度（又称表观速度、假想速度）是流体通过多孔介质整个横截面积（包括孔隙和骨架）的宏观平均速度；而真实渗流速度是流体仅在孔隙通道内流动的平均速度。*二者联系：u=v/n。此式表明，真实速度大于渗流速度，孔隙度越小，差距越大。*渗透率计算：根据达西定律Q=kAΔp/(μL)，其中Δp=p₁-p₂。因此，k=QμL/(AΔp)。注意各物理量的单位需统一（如采用SI单位制或油气行业常用单位制）。*达西定律适用条件：主要适用于层流渗流；流体为不可压缩牛顿流体；多孔介质为均质各向同性；渗流速度较小，惯性力可以忽略。超出此范围，可能进入非线性渗流区域。习题2：达西定律的修正与适用界限题目：当渗流速度较大时，达西定律不再适用，此时渗流速度v与压力梯度Δp/L的关系呈现非线性。试写出一种常见的描述这种非线性渗流规律的经验公式，并指出该公式中各参数的含义。如何通过实验方法判断渗流是否遵循达西定律？解析思路：*描述高速非线性渗流的经验公式，最著名的是福希海默（Forchheimer）方程：Δp/L=aμv+bρv²。其中，a和b为实验常数，与多孔介质和流体性质有关；μ为流体粘度；ρ为流体密度。等式右边第一项为粘性阻力项，第二项为惯性阻力项。当v较小时，第二项可忽略，方程退化为达西定律（a=1/k）。*判断达西定律适用性的实验方法：通过改变流速Q（进而改变v），测量对应的Δp。然后绘制v-Δp/L曲线。若为过原点的直线，则符合达西定律；若为曲线或不过原点的直线，则偏离达西定律。或者绘制Q（或v）与Δp的关系曲线，线性关系为达西流特征。二、单相不可压缩液体稳定渗流说明：本部分重点考察单相不可压缩流体在均质、各向同性或各向异性介质中，以稳定流方式运动时的压力分布、流量计算等问题，主要涉及拉普拉斯方程的建立与求解。习题3：平面单向稳定渗流题目：某均质各向同性无限大含水层，其中有一完整直线供水边界，边界上压力为pₑ。在距边界距离为r₀处有一口生产井，以恒定产量Q开采。假设流体为不可压缩液体，渗流为平面单向稳定渗流，地层厚度为h，渗透率为k，流体粘度为μ。试推导井壁处（r=r₀）的压力p_w的表达式。若已知pₑ、Q、k、μ、h、r₀，如何计算从边界到井壁的压力降落值？解析思路：*对于平面单向稳定渗流，基本微分方程（在一维直角坐标系下，x方向）为d²p/dx²=0。其通解为p=C₁x+C₂。*但本题描述的是“距边界距离为r₀处有一口生产井”，此处“r₀”可能为“x₀”更为贴切，因为单向流通常对应线性坐标。假设边界位于x=0，压力为pₑ；井位于x=x₀处，井壁压力为p_w。*渗流速度v=Q/(A)=Q/(h*1)（取单位宽度）。根据达西定律v=-k/μdp/dx。*分离变量积分：∫(p_wtopₑ)dp=-(μv/k)∫(x₀to0)dx。*解得：pₑ-p_w=(μv/k)x₀=(μQx₀)/(khx₀)？不，此处需仔细。v=Q/(h*B)，B为宽度，在单位宽度下B=1，故v=Q/h。则pₑ-p_w=(μv/k)*x₀=(μQx₀)/(kh)。因此，井壁压力p_w=pₑ-(μQx₀)/(kh)。*压力降落值Δp=pₑ-p_w=(μQx₀)/(kh)。*（注：若题目本意是平面径向流，但描述为“单向”，则可能存在表述混淆。需仔细审题。若为径向流，则公式为Δp=(μQ)/(2πkh)ln(rₑ/r_w)，其中rₑ为供给半径，r_w为井半径。此处按题目明确的“平面单向稳定渗流”作答。）习题4：平面径向稳定渗流与影响半径题目：一口生产井位于均质各向同性圆形油藏中心，油藏供给边界半径为rₑ，边界压力为pₑ。井半径为r_w，井底流压为p_wf。油藏厚度为h，原油粘度为μ，地层渗透率为k。原油在地层中作平面径向稳定渗流。(1)试推导该井产量Q的表达式。(2)若已知Q、k、h、μ、pₑ、p_wf，如何定义和计算该油藏的“影响半径”？影响半径的物理意义是什么？解析思路：*(1)平面径向稳定渗流的基本微分方程为(1/r)d/dr(rdp/dr)=0。其通解为p=C₁lnr+C₂。*应用边界条件：r=rₑ时，p=pₑ；r=r_w时，p=p_wf。*代入通解得：C₁=(p_wf-pₑ)/ln(r_w/rₑ)=(pₑ-p_wf)/ln(rₑ/r_w)。*根据达西定律，渗流速度v_r=-k/μdp/dr。产量Q=2πrhv_r=-2πrh(k/μ)dp/dr。*将dp/dr=C₁/r代入，得Q=2πkh(pₑ-p_wf)/(μln(rₑ/r_w))。此即平面径向稳定渗流的产量公式，是油藏工程中的核心公式之一。*(2)影响半径（rₑ）通常指在定产量生产时，从井中心到压力基本未受扰动（仍为原始地层压力pₑ）的距离。在稳定渗流中，它是一个确定的边界。*若已知Q等参数，可通过产量公式反求影响半径：ln(rₑ/r_w)=2πkh(pₑ-p_wf)/(μQ)，进而rₑ=r_wexp[2πkh(pₑ-p_wf)/(μQ)]。*物理意义：影响半径表示了渗流场的作用范围。在影响半径之外，地层压力不受井的生产活动影响。它综合反映了地层渗透性、流体粘度、产量大小以及生产压差等因素的影响。三、单相微可压缩液体不稳定渗流说明：本部分主要考察对不稳定渗流基本概念、压力传导方程、以及典型解（如瞬时点源解、Theis公式、无限大地层定产条件下的压力降落解）的理解与应用。习题5：不稳定渗流的压力传导方程题目：试推导单相微可压缩液体在均质各向同性多孔介质中流动时的压力传导方程（连续性方程与运动方程、状态方程相结合）。指出方程中导压系数（压力传导系数）α的表达式及其物理意义。解析思路：*推导过程是渗流力学的经典内容，需综合运用：1.连续性方程（质量守恒）：∂(ρn)/∂t+∇·(ρv)=0。2.运动方程（达西定律）：v=-(k/μ)∇p。3.状态方程（液体微可压缩）：ρ=ρ₀[1+C_L(p-p₀)]，其中C_L为液体压缩系数；多孔介质压缩系数C_f：n=n₀[1+C_f(p-p₀)]。通常忽略孔隙度随压力的变化，或合并为综合压缩系数C_t=C_L+C_f。*对上述方程进行线性化和化简（利用小扰动假设，忽略高阶小量），最终可得到以压力p为变量的扩散型方程，即压力传导方程：∇²p=(1/α)∂p/∂t其中，拉普拉斯算子∇²在不同坐标系下有不同表达形式（直角、柱坐标、球坐标）。*导压系数α=k/(μC_t)。其物理意义是：表征压力扰动在多孔介质中传播速度的快慢。α值越大，压力波传播得越快，地层压力达到稳定的时间越短。它综合反映了岩石渗透性、流体粘度和系统压缩性的影响。习题6：无限大地层定产条件下的不稳定渗流题目：均质各向同性无限大地层中存在一口生产井，假设地层原始压力为p_i，井以恒定产量Q投产。试写出描述井底压力p_wf随生产时间t变化的表达式（即不稳定渗流的井底压力降落公式），并说明公式中各参数的含义及所采用的假设条件。若已知某时刻t的井底压力p_wf，如何利用该公式估算地层渗透率k？解析思路：*无限大地层定产条件下的不稳定渗流，其典型解是基于瞬时点源解通过积分得到的。对于均质各向同性地层，考虑井径无限小的线源解，井底压力降落公式（以无因次形式或有量纲形式表示均可，此处以常用的有量纲形式为例，采用SI单位制时需注意常数项）为：p_i-p_wf(t)=(μQ)/(4πkh)*Ei(-r_w²/(4αt))其中，Ei(-x)为指数积分函数，x=r_w²/(4αt)=μC_tr_w²/(4kt)。当时间t足够大，x较小时，Ei(-x)≈-γ-lnx（γ为欧拉常数，约为0.5772），此时公式可简化为：p_i-p_wf(t)=(μQ)/(4πkh)[ln(4αt/r_w²)-γ]或p_wf(t)=p_i-(μQ)/(4πkh)ln(t/t₀)（其中t₀为包含常数项和已知参数的组合）。这就是著名的霍纳（Horner）近似或常规压降公式的简化形式。*参数含义：p_i为原始地层压力，p_wf(t)为t时刻井底流压，μ为流体粘度，Q为产量，k为地层渗透率，h为地层厚度，r_w为井半径，α为导压系数，t为生产时间。*假设条件：地层均质、各向同性、无限大；原始压力均匀；流体为微可压缩牛顿流体；达西定律适用；井为完善井（裸眼完井，钻开全部地层厚度），产量恒定；忽略井储效应和表皮效应（或认为表皮系数S=0）；井以定产量投产，投产前地层处于平衡状态。*估算地层渗透率k：若已知某时刻t的井底压力p_wf，可利用简化后的对数形式公式。通过测量不同时刻的井底压力数据，绘制p_wf与lnt的关系曲线（MDH曲线），其斜率m=(μQ)/(4πkh)。因此，k=(μQ)/(4πhm)。实际应用中，常通过直线段斜率计算。四、多相渗流基础说明：多相渗流比单相渗流更为复杂，本部分主要考察相对渗透率的概念、影响因素以及两相渗流的基本规律。习题7：相对渗透率曲线及其应用题目：什么是相对渗透率？简述油、水两相相对渗透率曲线的一般特征。若某油藏中同时存在油和水两相渗流，已知油相粘度为μ_o，水相粘度为μ_w，油相相对渗透率为k_ro，水相相对渗透率为k_rw，总压力梯度为Δp/L。试写出油相渗流速度v_o和水相渗流速度v_w的表达式。解析思路：*相对渗透率：多相流体共存并渗流时，某一相流体的有效渗透率与该介质绝对渗透率的比值，称为该相流体的相对渗透率。它无量纲，取值范围在0到1之间。*油水两相相对渗透率曲线特征：*共存水饱和度S_wc（束缚水饱和度）：当含水饱和度S_w=S_wc时，k_rw=0，k_ro=k_ro_max（最大油相相对渗透率）。*残余油饱和度S_or：当含水饱和度S_w=1-S_or时，k_ro=0，k_rw=k_rw_max（最大水相相对渗透率）。*中间饱和度范围：随着含水饱和度S_w的增加，k_ro逐渐减小，k_rw逐渐增加。两条曲线通常会有一个交点，交点处k_ro=k_rw。*相对渗透率是饱和度的函数，且与润湿性、孔隙结构、流体分布等因素有关。*渗流速度表达式：油相：v_o=-(kk_ro/μ_o)(Δp/L)水相：v_w=-(kk_rw/μ_w)(Δp/L)总渗流速度v=v_o+v_w。注意此处压力梯度Δp/L对两相而言是相同的，即两相在同一压力梯度下流动。五、井流问题与试井分析初步说明：试井分析是渗流力学理论在现场应用的重要手段，本部分考察对试井基本概念和方法的理解。习题8：稳定试井与产能方程题目：简述稳定试井的基本原理和主要目的。某油井进行稳定试井，测得不同工作制度下的产量Q和对应的井底流压p_wf数据如下表（假设地层压力p_R为常数）：产量Q(m³/d)井底流压p_wf(MPa)----------------------------