控压钻井重力置换与溢流气侵判断准则分析

考虑重力置换的恒定性、气侵的可变性,提出了重力置换与溢流气侵的判断三准则,实现了计算机自动判断气体来源。借助计算机编程对准则求解,并将判断准则应用于某微流量控压钻井的溢流气侵和重力置换判断。结果表明:重力置换发生后,一定量的置换气体侵入环空,气体的侵入量主要受岩层渗透率、钻井液-岩层接触面积的影响,而溢流气侵发生后,气体的侵入量主要受地层欠压差影响;为达到衡压钻井目的,溢流气侵发生后不但要抑制气体滑脱产生的压降,更要平衡地层欠压差,而重力置换只需抑制气体滑脱产生的压降;溢流气侵发生后,环空中流体的流动先由0气流转换为混合流,经回压控制后,再由混合流转换为0气流,重力置换发生后,环空中流体的流动由0气流转换为混合流。     

不同重力条件下管内冷凝现象研究进展

本文对不同重力环境中水平管内冷凝气液两相流动与传热现象的研究现状进行了全面评介, 重点关注于航天应用中的小管径、低流量和以氨为工质的情形,以及部分重力条件下冷凝现象中的重力效应. 冷凝现象中管壁四周液膜的存在,导致水平管内冷凝气液两相流与沸腾或绝热气液两相流在构型上存在着明显的不同, 表面张力的作用增大,向波状和分层流转换所对应的临界Bond数也随之增大, 进而使小管径、低流量水平管内冷凝气液两相流动与传热现象中的重力效应减弱,甚至可以忽略不计. 实验发现该情形中摩擦压降要小于基于常规气液两相流实验数据的经验关联式的预测结果, 而更接近光滑环状流情形; 对冷凝两相传热系数的预测,只有采用基于空隙率与界面剪切率相互耦合的模型, 才能够给出较好的结果.      

柔性网架结构的超弹性力学行为分析

由于地层漏失压力存在差异性,钻完井液的漏失现象是工程中常遇问题。井壁柔性贴补网架是一种解决地层钻井溃漏的方案,其中核心涉及到网架结构超弹性的几何大变形力学行为。本文就柔性网架结构建立了子结构的非线性力学模型,推导了结构位移形变方程,证明了网架结构的几何非线性,给出了钻井堵漏的力学行为算法。这种几何非线性力学模型,通过数值方法得到了验证。该柔性网架的力学分析模型适用于钻井堵漏,也可以为心脑血管SMA支架、气管支架等领域提供参考。     

不同重力下90°弯管内气液两相流流型及流动特性研究

研究了不同重力条件下90°弯管内气液两相流流型分布形态及流动特性. 通过建立90°弯管内气液两相流流动的三维数学物理模型,采用VOF 方法,对10-6g₀, 10-4g₀, 10-2g₀, 1g₀ (g₀= 9.8m/s2) 重力下的90°弯管内气液两相流流型分布特征、截面空隙率、滑速比及气相尾部最大斜向角进行了比较分析. 研究结果表明:所建立的模型能够正确模拟不同重力条件下90°弯管内气液两相流流型和截面空隙率,并得到气液两相弯管二次流与单相二次流的不同特性. 随着重力水平的提高,90°弯管对气相流型的影响作用减弱,气相整体向弯管内侧积聚靠拢,弯管对尾部的斜向作用减弱.      

矩形弹性壳液耦合系统中的重力波分析

根据非线性动力学理论,建立了矩形壳液耦合系统的非线性振动方程组,通过数值求解,发现当激振频率为壳体固有频率与重力波频率之和,且激振力足够大时,会产生大幅低频重力波,通过实验验证,发现了壳液耦合系统中存在的大幅低频重力波现象,实验结果与理论结果基本吻合。     

微重力条件下气/液两相流流型的研究进展

气／液两相流流型是两相流研究领域最基本的课题之一,至今已有数十年的研究历史．但是,由于气／液两相流动现象极为复杂,目前还没有得到一致的结论．近十多年来,利用微重力环境减弱甚至完全消除重力的影响,简化流动中各种不同因素间的相互作用及流型特征,大大促进了对气／液两相流动特征及流型产生与转换机理的研究．同时,微重力条件下的气／液两相流动是空间技术领域必须解决的关键技术问题之一,具有重要的学术意义和重大的应用价值．本文简要总结了微重力条件下气／液两相流流型研究的基本方法以及实验结果和理论进展,指出今后研究中应该注意的一些方向．      

I-II复合型裂缝断裂角的分层剪滞模型与分析

采用修正的剪滞理论建立了岩石、混凝土等准脆性材料的I-II复合型裂缝在单向拉伸荷载作用下的计算模型,得到了与实验相吻合且优于传统S判据的断裂角。通过对远场应力、斜裂缝区应力以及子层位移的合理简化,得到了求解剪滞分析模型的边界条件,进而得到了含斜裂缝的各子层位移分布函数。引入最大应力集中因子,对I-II复合型裂缝前缘应力场进行简化;基于斜裂缝沿最大应力集中因子方向扩展,得到裂缝的断裂角。根据斜裂缝的应力分布,设置不同的子层分区,得到了更为细化的位移分布模式。通过对计算数据的分析,针对单向拉伸荷载作用下的I-II复合型裂缝,建立了按应力场分区设置子层的分层剪滞模型,得到更为精确的斜裂缝断裂角。     

Ⅰ-Ⅱ复合型裂缝断裂角的分层剪滞模型与分析

采用修正的剪滞理论建立了岩石、混凝土等准脆性材料的Ⅰ-Ⅱ复合型裂缝在单向拉伸荷载作用下的计算模型.得到了与实验相吻合且优于传统S判据的断裂角.通过对远场应力、斜裂缝区应力以及子层位秽的合理简化,得到了求解剪滞分析模型的边界条件,进而得到了含斜裂缝的各子层位移分布函数.引入最大应力集中因子,对Ⅰ-Ⅱ复合型裂缝前缘应力场进行简化;基于斜裂缝沿最大应力集中因子方向扩展,得到裂缝的断裂角.根据斜裂缝的应力分布,设置不同的子层分区,得到了更为细化的位移分布模式.通过对计算数据的分析,针对单向拉伸荷载作用下的Ⅰ-Ⅱ复合型裂缝,建立了按应力场分区设置子层的分层剪滞模型,得到更为精确的斜裂缝断裂角.     

自适应Kalman和零相移滤波算法在重力信号处理中的对比

为了有效消除海洋重力仪测量信号的噪声,提高重力数据的获取精度,根据随机过程理论,借助基于二阶高斯一马尔可夫异常位模型的重力异常协方差函数,得到海洋重力测量中重力异常信号的状态方程.对sage-husa滤波算法和零相移算法进行了理论分析,为了抑制零相移滤波器的首尾数据畸变,求解了滤波器的初始状态.根据实测重力数据,进行了去噪仿真试验.理论分析和仿真结果表明,sage-husa与零相移滤波算法均能较好地抑制采样重力数据中的噪声干扰,但零相移滤波算法的性能优于sage-husa滤波器.     

基于SPSS的重力匹配区域选择算法

为了提高重力辅助导航系统的定位精度,对重力匹配区域的选择进行了研究,提出一种基于SPSS回归分析的重力匹配区域选择方法。通过对重力数据进行回归分析,得到匹配区域的回归函数关系式,以此作为重力匹配判断准则对重力匹配区域进行选择。通过惯导/计程仪/重力组合导航系统的仿真试验,证明基于SPSS的重力匹配区域选择方法能找出重力匹配较好的区域,能提高重力辅助导航系统的定位精度。     

Modeling of Solute Transport in a 3D Rough-Walled Fracture–Matrix System

Fluid flow and solute transport in a 3D rough-walled fracture–matrix system were simulated by directly solving the Navier–Stokes equations for fracture flow and solving the transport equation for the whole domain of fracture and matrix with considering matrix diffusion. The rough-walled fracture–matrix model was built from laser-scanned surface tomography of a real rock sample, by considering realistic features of surfaces roughness and asperity contacts. The numerical modeling results were compared with both analytical solutions based on simplified fracture surface geometry and numerical results by particle tracking based on the Reynolds equation. The aim is to investigate impacts of surface roughness on solute transport in natural fracture–matrix systems and to quantify the uncertainties in application of simplified models. The results show that fracture surface roughness significantly increases heterogeneity of velocity field in the rough-walled fractures, which consequently cause complex transport behavior, especially the dispersive distributions of solute concentration in the fracture and complex concentration profiles in the matrix. Such complex transport behaviors caused by surface roughness are important sources of uncertainty that needs to be considered for modeling of solute transport processes in fractured rocks. The presented direct numerical simulations of fluid flow and solute transport serve as efficient numerical experiments that provide reliable results for the analysis of effective transmissivity as well as effective dispersion coefficient in rough-walled fracture–matrix systems. Such analysis is helpful in model verifications, uncertainty quantifications and design of laboratorial experiments.     

Numerical study on the shear-flow behavior and transport process in rough rock fractures

This paper presents a systematic research for understanding mechanical shearing effects on the fluid flow and the solute transport behavior of rough fractures through a numerical simulation approach. The aperture fields were modeled based on a real rock fracture geometry and the normal displacement obtained from the shear-flow test. The fluid flow through the rough fracture under shear was simulated using a finite element code that solves the Reynolds equation, and the transport behavior through the rough fracture under shear was simulated calculating the advection–dispersion equation. The results show that the fracture apertures increase as the shear displacement increases, with a few major flow channels detected through the fracture. The shear-induced flow channels increase both flow connectivity and transport connectivity, which accelerate the movement of solutes in a particular direction and lead to early breakthrough of the contaminants. Adsorption, acting as a retardation term, has a decisive influence on the transport process. These results can give a basic knowledge of the hydromechanical and solute transport progress through fracture, and will be helpful to safety assessment for high-level radioactive waste disposal facilities.     

微重力条件下部分充液贮箱气液界面波动特性的数值模拟

微重力环境中部分充液贮箱内气液界面和气、液两相介质在残余重力或加速度干扰下的运动特征是先进空间流体管理技术的基础. 本文针对空间贮箱常用构型和实际尺寸, 基于邦德数相似准则设计了3个缩比模型, 数值模拟了原型贮箱和缩比模型中加速度变化引起的贮箱内气液两相流动及气液界面上界面波的传播. 数值模拟结果验证了原型和模型间的运动相似性, 发现在满足邦德数相似准则的前提下, 系统还近似满足韦伯数相似准则, 或等价地, 近似满足弗劳德数相似准则. 此外, 数值模拟结果也表明原型和模型间的运动存在细微偏差, 这主要源于黏性耗散作用的差异. 由韦伯数相似准则可知, 缩比加大, 贮箱尺寸减小, 重力突变后由表面张力释放出来的驱动力增强, 相同韦伯数下流动速度增大, 黏性耗散作用随之增强, 本文的数值模拟结果证实了该结论. 相关结果可以用于指导空间贮箱流体管理技术的地面模拟试验的方案设计等.      

近场水下爆炸瞬态强非线性流固耦合无网格数值模拟研究

近场水下爆炸涉及多相流体的掺杂耦合以及结构的大变形、损伤和断裂等瞬态强非线性现象, 传统的网格算法在模拟近场水下爆炸时面临结构网格畸变、多相界面捕捉精度不足等难题, 鉴于此, 本文建立了完全无网格的近场水下爆炸冲击波和气泡全物理过程瞬态强非线性流固耦合动力学模型. 流体采用基于黎曼求解器的光滑粒子流体动力学(SPH)方法求解, 结构采用重构核粒子法(RKPM)求解, 并基于法向通量边界条件实现流固耦合. 为提高SPH对流场间断的求解精度, 引入黎曼问题思想并结合MUSCL重构算法, 为解决流场粒子体积变化剧烈导致的精度下降问题, 应用了自适应粒子分割与合并方法. 为模拟水下爆炸对结构造成的损伤断裂, 基于退化实体几何表述, 采用Lemaitre损伤算法, 建立了RKPM壳结构断裂损伤模型. 依据所建立的SPH-RKPM流固耦合模型, 对近场水下爆炸冲击波传播、气泡脉动与射流以及结构毁伤进行了模拟, 将得到的冲击波载荷、气泡演化以及结构响应与实验值和其他数值解对比, 验证了当前建立的SPH-RKPM流固耦合模型的有效性和精度, 并给出了水下爆炸载荷特性及其对结构的流固耦合毁伤机制与规律, 旨在为近场水下爆炸载荷预报提供理论和基础性技术支撑, 为毁伤威力评估和舰船防护结构设计提供参考.      

A relationship between the fractal dimension and scaling groups of unstable miscible displacements

A detailed two-dimensional flow visualization study was performed to examine the dynamics of viscous fingering in miscible displacements. Detailed quantitative miscible displacement experiments using a microcomputer-based imaging workstation on a variety of oil recovery fluid systems were performed. The effect of two dimensionless scaling groups, namely gravity number and viscosity ratio, on the displacement behavior was investigated. Based on image analysis, the irregular fingering patterns of the flow visualization experiments were analyzed for fractal characteristics. Results indicate that the areal sweep efficiency of unstable miscible displacement follows a fractal scaling law with a fractal dimension and proportionality constant related to the gravity number and the viscosity ratio. The study shows that the fractal dimension decreases with decreasing gravity number and increasing viscosity ratio. This relationship was mapped by an artificial neural network model, which can be used to estimate the fractal dimension and the proportionality constant of miscible displacements as functions of the two scaling groups. These results have potential application in the mathematical modeling of unstable EOR displacements and in the scaling of laboratory displacements to field conditions. Received: 19 December 1999/Accepted: 26 January 2001     

Development of hot dry rocks by hydraulic stimulation: Natural fracture network simulation

This paper presents an advanced computational modelling of natural fracture networks in HDR (hot dry rock) reservoirs. The model stochastically simulates discrete properties of natural fractures, utilizing multi-set orientation and fractal mathematics. The simulated fracture networks are essential for further stimulation and fluid flow studies. The model has been verified using the data of actual fracture stimulation programs conducted by Gas Research Institute and Department of Energy at the multi-site. It is validated that the simulated fracture distribution is sufficiently similar to that observed in the reservoir. This paper also examines the detrimental effects of the simulated natural fracture network on the stimulated fluid flow capacity. The effective permeability enhancement (due to hydraulic stimulation) is found almost proportional to the density of the reservoir natural fractures.     

井筒气液两相流对致密气压裂水平井试井的影响

多段压裂水平井技术是目前开采致密气最常用的方法之一,在致密气压裂水平井试井测试中常常伴随着一定的产水量,井筒气液两相流会增加井筒流体的流动阻力,加大井筒流体流动对试井解释的影响.为了明确井筒气液两相流对致密气藏压裂水平井试井的影响,提高产水致密气压裂水平井的试井解释精度,建立了一种井筒气液两相流与地层渗流耦合的试井模型,采用数值方法对模型进行求解,获得了考虑井筒气液两相流的压裂水平井试井理论曲线、压力场分布及裂缝产量分布.研究结果表明:井筒气液两相流会增加试井理论曲线中压力和压力导数值,造成靠近入窗点的压力扩散要快于远离入窗点的压力扩散,引起靠近入窗点的裂缝产量要高于远离入窗点的裂缝产量.现场实例分析进一步说明,不考虑井筒两相流可能会对产水压裂水平井的试井解释结果产生很大误差,主要表现为水平井筒假设为无限大导流能力会使得拟合得到的表皮系数偏大,将测试点视为入窗点会使得拟合得到的原始地层压力偏小.所建立的考虑井筒两相流的压裂水平井试井模型为产水致密气井试井资料的正确解释提供了重要技术保障.      

单组分多相系统驱替过程的格子Boltzmann模拟

结合格子Boltzmann方法中的Shan-Chen单组分多相模型,引入流体相间的内聚力和流体与 固体壁面间的黏附力,对二维孔隙网格中非浸润气相驱替完全浸润液相的过程进行模拟,流 体相间的交界面自然形成,整个驱替过程属于毛细指进. 随着毛细数的增加,黏性力的主导 作用增强,使得气相入侵的孔隙尺度减小,因此驱替形态随毛细数的不同有很大差别. 在微 重力的作用下,整个驱替过程受毛细力、重力和黏性力的共同作用,重力起到了稳定交界面 避免窜流的作用.     

Radial Flow of Non-Newtonian Power-Law Fluid in a Rough-Walled Fracture: Effect of Fluid Rheology

Fluid flow in a single rough-walled rock fracture has been extensively studied over the last three decades. All but few of these studies, however, have been done with Newtonian fluids and unidirectional flow in rectangular fractures. Notwithstanding the importance of such setups for theoretical understanding of fundamental issues in fracture flow, practical applications in drilling and petroleum engineering often involve radial flow of a non-Newtonian fluid. An example is a borehole intersecting a natural fracture during drilling in a fractured rock. In this study, steady-state incompressible radial flow from a circular well into a self-affine rough-walled fracture was simulated numerically using the lubrication theory approximation. The fluid rheology was power law. The flow behavior index was equal to 0.6, 0.8, 1.0 (Newtonian), 1.2, or 1.4. Asperities diverted the flow from an axisymmetric radial pattern that would be observed in a smooth-walled fracture. The extent of the deviation from radial flow was found to increase as the fluid became more shear-thickening. To reveal finer details of the flow, a tracer was introduced at the borehole wall and was transported by the flow. The front of the tracer propagating into the fracture was found to become slightly smoother with a more shear-thickening fluid. In the vicinity of contacts between fracture faces a more shear-thickening fluid could deliver the tracer closer to the contact spots.