水平井的水力特征及其解析解的适用条件

利用自制的水平井砂槽模型, 进行了一系列不同流量条件下的水平井抽水试验, 结果表明: (1)在不同出流条件下, 水平井井管中可以同时出现层流-粗糙紊流多种不同流态; (2) 水平井出流条件下, 井管中的水头损失既不能忽略, 也不服从线性变化规律.它与井管中的水流流态有关.用“等水头井壁”或“等强度线汇”来刻画水平井井壁边界条件是不全面的.根据“等强度线汇”理论得到的解析解与试验结果对比发现, 本试验条件下解析解的近似适用条件是: 水平井管中的水流全部为层流(Re < 2 320)或者层流和层流-光滑紊流过渡区(Re < 4 000)同时并存的情况.当水平井管中出现光滑紊流区(Re> 4 000), 即同时有层流、层流-光滑紊流过渡态和光滑紊流或更多种流态时, 解析解已不再适用, 此时必须用新的层流-管流耦合模型来求解.      

裂缝性低渗透油藏流-固耦合数值模拟

根据低渗裂缝性砂岩油藏的储层特征,建立适合裂缝性砂岩油藏渗流的等效连续介质模型;通过实验得出有效压力改变对储层孔隙度、渗透率影响的规律;考虑低渗透油藏渗流时启动压力梯度和低渗储层的流-固耦合特性,将渗流力学与弹塑性力学相结合,建立起适合低渗透油藏的流-固耦合渗流数学模型,并给出其数值解.在黑油模型和弹塑性有限元程序的基础上,编制了计算低渗透油藏流-固耦合渗流的计算软件.通过数值模拟和不考虑流-固耦合时的计算结果相比,可以看出低渗油藏中流-固耦合效应十分明显.     

单向流边界层泥沙起动规律

Shields曲线常用于表示泥沙起动的临界条件,基于边界层理论,对Shields曲线各个流区的线型进行了推证;考虑粘结力的作用,对Shields参数及Shields曲线进行了修正,并给出修正Shields曲线表达式;在此基础之上,从边界层角度重新阐述了Shields曲线。结果表明:Shields曲线在光滑紊流及层流区呈直线分布,在过渡区与阻力系数线型保持一致,在粗糙紊流区呈水平直线分布;修正后Shields曲线与原始Shields曲线在形式上保持一致,修正Shields曲线表达式与实测数据吻合较好,适用于粗、细颗粒泥沙起动条件的计算;Shields曲线事实上代表了Shields参数与沙粒周围绕流流态的关系,同一颗粒处于不同流区起动时,其起动切应力不同。     

土质边坡的饱和−非饱和渗流分析及特殊应力修正

基于连续多孔介质原理与混合体原理引入流?固耦合模型,采用FISH语言编制程序,建立了初始渗流场设置函数、饱和?非饱和渗流分析模块、特殊应力修正模块以及非饱和单元抗剪强度修正模块,实现了利用FLAC3D进行饱和?非饱和渗流分析,同时基于饱和?非饱和渗流场修正非饱和区土体单元的有效应力以及抗剪强度,完成了将FLAC3D中饱和土流?固耦合计算原理扩展到非饱和土中。在二次开发的分析模块基础上对Liakopoulos砂柱排水试验进行数值模拟,验证了计算模型和计算程序的正确性。同时也研究了降雨入渗对土质边坡的渗流场、位移场以及稳定性演变过程的影响,揭示了降雨诱发边坡失稳的机制。     

地下水封油库围岩地下水渗流量计算

为准确地计算特大型地下水封储油岩洞库地下水的渗流量,根据现场调查、测试与监测以及室内外物理力学试验,采用了三维多孔连续介质流-固耦合有限差分数值模拟计算分析软件模拟地下水渗流场。提出了一种适用于特大型地下水封储油岩洞库地下水渗流量计算的合理方法。经与法国专家经验法、大岛洋志经验式、《铁路工程水文地质勘测规范》（TB10049－96）等方法的计算结果进行比较,结果比较接近。该计算方法适用于试验数据和资料匮乏的可行性研究阶段的特大型地下水封储油岩洞库围岩地下水渗流量或涌水量的计算;考虑地下水流-固耦合分析时,单洞室的涌水量比只考虑地下水流动作分析时稍小;考虑流-固耦合作用比不考虑流-固耦合作用时围岩地下水的最大渗流速度稍小。     

底水油藏水平井分段完井产能预测

在分析分段完井渗流机理的基础上,根据质量守恒原理及动量定理推导出了水平井分段完井井筒压降计算模型,并引入了一个修正的壁面摩擦系数,可以准确计算分段后各微元段的压降.针对底水油藏类型,根据势的叠加原理导出了单相原油三维稳态渗流的压力分布公式,再结合分段压降模型,建立了油藏渗流与井筒流动耦合产能预测模型.实例计算表明:用此...     

层流—紊流共存流场中岩溶裂隙网络演化过程的数值模拟方法

研究岩溶水系统的演化过程对许多资源与环境问题有着重要意义。岩溶地区水资源、油气资源的预测与开发、水土流失成因及防治等课题都与此密切相关。岩溶管道是地下水对裂隙的逐渐溶蚀扩展形成的,岩溶演化初期,所有裂隙宽度不大,流场整体成层流状态。随着溶蚀的进行,一部分裂隙优先扩大使其中的水流进入紊流状态,而另一部分裂隙中水流仍呈层流状态。文章提出一种数值方法,能够模拟层流-紊流共存流场的岩溶演化过程。利用蒙特卡洛(Monte-Carlo)方法模拟初始裂隙网络,通过非连续介质方法模拟裂隙网络中的渗流。裂隙扩展的速度通过岩石表面溶蚀速度经验公式计算,使用迭代方法求解层流-紊流共存条件下的水头非线性方程。构建了能够利用解析法求解的模型,把数值解和解析解的结果进行对比,验证了本研究的数值法及软件的可行性。      

裂隙岩体流-热耦合传热的三维数值模拟分析

通过对潘西煤矿水文地质条件的分析,基于裂隙岩体的流-热耦合数学模型,描述了裂隙岩体渗流场分布和水流及岩体的温度场分布,并结合边界条件及计算参数对裂隙岩体的流-热耦合传热进行了数值模拟和分析。数值模拟结果表明,岩体内裂隙水流所引发的热量迁移,对裂隙岩体的温度场分布有重要影响。断裂带及地下水流的存在改变了岩体的原有温度场分布。在渗流初期,温度梯度矢量沿渗流方向向两侧岩体方向流动,由于两侧岩体的渗透性系数低于断裂带处的渗透性系数,右侧等温线及温度梯度矢量方向逐渐向渗流方向移动,改变了两侧岩体的温度场分布。通过对断裂带内裂隙水流渗透性系数的折减,分析渗透性系数发生变化时对岩体温度场分布的影响,渗透性系数越大,伴随的热量迁移增大,对岩体的温度场分布的影响也越大。     

钻井液流态对孔壁稳定性的影响

沉积地层中孔壁的稳定性与钻井液流态有着密切关系。结合河南叶舞凹陷盆地2 386 m盐矿钻探工程,采用不同的泵量和上返流速在不同孔段进行了钻探试验。试验和理论计算表明:在水敏性地层（黏土、泥岩、砂质泥岩等）浅层岩石孔隙率和塑性指数较高时,在钻井液层流流态情况下,主要以缩径为主,最大缩径率达26.6%;中部和深部同等地层,在紊流流态情况下,缩径和超径同时存在。其中,随着孔深的增加,缩径率呈减小趋势,超径率则有增大趋势,泥岩地层最大超径率达157.4%,盐岩地层超径率高达225.7%。      

急倾斜上保护层开采瓦斯越流固-气耦合模型及保护范围

开采保护层是防治煤与瓦斯突出最有效的措施之一,其关键是保护范围的合理确定。针对急倾斜多煤层上保护层开采有效保护范围的划定问题,基于煤层瓦斯越流理论,根据煤岩层变形与瓦斯渗流的固-气耦合作用,建立了瓦斯渗流场方程和煤岩体变形场方程,得到了急倾斜上保护层开采瓦斯越流固-气耦合数学模型。以南桐矿区某矿上保护层开采为实例,通过多物理场耦合系统,建立了该矿上保护层开采瓦斯越流几何模型并进行数值计算,获得了上保护层工作面开采后被保护层瓦斯压力的分布规律,确定了上保护层开采的卸压保护范围。数值计算与现场考察试验结果具有一致性,由此验证了数值计算的合理性。研究结果可以对现场保护范围的划定及卸压瓦斯抽放等提供理论指导,具有实际工程意义。     

地下水水平井流的模型及数值模拟方法—考虑井管内不同流态

评述用线汇刻画水平井管可能存在的问题, 并用理想模型证明.提出考虑井管不同流态(线性流与非线性流) 水流阻力的水平井流的模型, 并将内边界取在水平井的出水口处, 避免了采用线汇刻画水平井管时假定其流量分布和水头分布的困难.运用作者提出的等效渗透系数的概念与确定方法, 将水平井-含水层系统视为含有圆柱形透镜体的非均质含水层, 给出新的水平井流的数学模型.求解了一个理想模型, 得出河下水平井开采地下水的动态.      

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Fluvial and submarine morphodynamics of laminar and near‐laminar flows: a synthesis

The interaction of flow with an erodible bed in alluvial rivers and deep‐sea channels gives rise to a wide range of self‐formed morphologies, including channels, ripples, dunes, antidunes, alternate bars, multiple‐row bars, meandering and braiding. As the flow is invariably turbulent in field manifestations of these morphologies, there has been a tendency to assume that turbulence is necessary for them to form. While turbulence undoubtedly has an important influence when it is present, it is not necessary for any of these features. Indeed, all of these features can be formed by the morphodynamic interaction of purely laminar or nearly laminar flow with an erodible bed. This paper provides a survey and synthesis of a wide range of laminar or near‐laminar flow analogues of morphologies observed in the field. Laminar‐flow analogues of turbulent‐flow morphologies cannot and should not be expected to satisfy dynamic similarity in terms of all relevant dimensionless parameters. What is of more significance is the convergence of the underlying physics. It is illustrated in this paper that many existing theoretical frameworks for the explanation of turbulent‐flow morphodynamics require only relatively minor modification in order to adapt them to laminar flows.     

高水力梯度条件下颗粒堆积型多孔介质渗流规律试验研究

煤矿开采面临的水文地质条件越来越复杂,尤其是遭遇承压含水层的水压力越来越大,突水灾害发生时必然会带来高水力梯度引起的破碎岩体突水通道内高速非线性渗流问题。据此,研制高水力梯度（最大600）条件下堆积型多孔介质中高速非线性渗流试验装置,采用堆积型钢球模拟破碎岩体,对粒径为1、2、3、4、5、6 mm共6种光滑钢球分别开展了一维均质圆柱渗流试验。试验结果表明：对于由1～6 mm钢球堆积而成的孔隙率为0.44～0.45的多孔介质,当水力梯度大于145时,通过分析水力梯度-平均流速（J-v）曲线和水力梯度-雷诺数（J-Re）关系曲线,将流动状态划分为3个模式：线性层流、非线性层流、紊流,并获得了从线性层流过渡到非线性层流的临界流速为0.23～0.78 cm/s、临界水力梯度为3～8;从层流到紊流转捩的临界流速为1.6～4.8 cm/s、临界水力梯度为90～145。从小粒径多孔介质到大粒径多孔介质的渗流过程中,临界流速越来越大,而临界水力梯度逐渐减小。 渗透率与粒径的平方、非达西流影响系数与粒径的倒数均呈线性正相关,非达西流影响系数随着渗透率的增加呈指数减小。该研究对多孔介质非线性渗流的理论研究以及实际工程中高承压含水层突涌水问题有重要借鉴意义。     

磨料两相圆形射流紊流方程的建立及其求解

根据一定合理假设从单相圆形层流边界方程的建立和求解入手，分别建立了两相圆形层流射流边界层方程，单相圆形紊流边界层方程和两相圆形紊流边界层方程并给出了各种情况下的流速分布等参数的解答，为研究磨料两相射流对材料的打击和切割作用奠定了理论基础。     

An experimental study of single and two-phase fluid flow through fractured granite specimens

Triaxial tests on the two-phase flow of air and water through fractured granite specimens were performed to discover whether the two-phase fluid flow within rock fractures was laminar or turbulent. The two-phase flow characterization was carried out based on the macroscopic two-phase steady state flow model and the homogeneous steady state flow model. Rock specimens with a single natural fracture (joint roughness coefficient, JRC < 10) were tested using two-phase, high pressure triaxial rig. Experimental results show that the estimated Reynolds numbers for various inlet fluid pressures are well below 1000. The findings of this study reveal that both single and two-phase flow through rock fractures (JRC < 10) can be characterized as laminar flows at moderate inlet fluid pressures. However, for single-phase air flow, an increase in inlet air pressures may result in the formation of turbulent flow.     

Depositional processes,bedform development and hybrid bed formation in rapidly decelerated cohesive (mud–sand) sediment flows

Flows with high suspended sediment concentrations are common in many sedimentary environments, and their flow properties may show a transitional behaviour between fully turbulent and quasi‐laminar plug flows. The characteristics of these transitional flows are known to be a function of both clay concentration and type, as well as the applied fluid stress, but so far the interaction of these transitional flows with a loose sediment bed has received little attention. Information on this type of interaction is essential for the recognition and prediction of sedimentary structures formed by cohesive transitional flows in, for example, fluvial, estuarine and deep‐marine deposits. This paper investigates the behaviour of rapidly decelerated to steady flows that contain a mixture of sand, silt and clay, and explores the effect of different clay (kaolin) concentrations on the dynamics of flow over a mobile bed, and the bedforms and stratification produced. Experiments were conducted in a recirculating slurry flume capable of transporting high clay concentrations. Ultrasonic Doppler velocity profiling was used to measure the flow velocity within these concentrated suspension flows. The development of current ripples under decelerated flows of differing kaolin concentration was documented and evolution of their height, wavelength and migration rate quantified. This work confirms past work over smooth, fixed beds which showed that, as clay concentration rises, a distinct sequence of flow types is generated: turbulent flow, turbulence‐enhanced transitional flow, lower transitional plug flow, upper transitional plug flow and a quasi‐laminar plug flow. Each of these flow types produces an initial flat bed upon rapid flow deceleration, followed by reworking of these deposits through the development of current ripples during the subsequent steady flow in turbulent flow, turbulence‐enhanced transitional flow and lower transitional plug flow. The initial flat beds are structureless, but have diagnostic textural properties, caused by differential settling of sand, silt and cohesive mud, which forms characteristic bipartite beds that initially consist of sand overlain by silt or clay. As clay concentration in the formative flow increases, ripples first increase in mean height and wavelength under turbulence‐enhanced transitional flow and lower transitional plug‐flow regimes, which is attributed to the additional turbulence generated under these flows that subsequently causes greater lee side erosion. As clay concentration increases further from a lower transitional plug flow, ripples cease to exist under the upper transitional plug flow and quasi‐laminar plug flow conditions investigated herein. This disappearance of ripples appears due to both turbulence suppression at higher clay concentrations, as well as the increasing shear strength of the bed sediment that becomes more difficult to erode as clay concentration increases. The stratification within the ripples formed after rapid deceleration of the transitional flows reflects the availability of sediment from the bipartite bed. The exact nature of the ripple cross‐stratification in these flows is a direct function of the duration of the formative flow and the texture of the initial flat bed, and ripples do not form in cohesive flows with a Reynolds number smaller than ca 12 000. Examples are given of how the unique properties of the current ripples and plane beds, developing below decelerated transitional flows, could aid in the interpretation of depositional processes in modern and ancient sediments. This interpretation includes a new model for hybrid beds that explains their formation in terms of a combination of vertical grain‐size segregation and longitudinal flow transformation.