防砂井挡砂介质堵塞规律实验及堵塞程度定量预测模型

防砂井挡砂介质堵塞造成降产已成为困扰油气井防砂的关键问题。已开展的大量堵塞规律研究主要涉及机理和定性描述,尚未形成可用的堵塞程度定量预测方法。针对这一问题,本文在广泛应用的砾石充填防砂条件下,利用挡砂介质堵塞评价驱替实验装置,进行了砾石层堵塞机理和规律实验,着重考察了砾石层堵塞渗透率比随驱替时间、砾砂中值比、流体流速、流体粘度、粘土含量、细质含量等生产条件的定量变化规律。根据实验数据,首先拟合砾石充填层堵塞渗透率比随时间的定量关系,引入综合堵塞系数和单项堵塞系数表征各生产条件对堵塞规律的影响,然后拟合了堵塞系数与各生产参数的经验关系,整合得到一个根据防砂参数和生产条件预测砾石充填层堵塞程度随时间变化规律的定量模型。结果表明,不同生产条件砾石层渗透率在投产早期下降较快,然后下降速度减缓并趋于基本不再变化,达到堵塞平衡状态。具体变化规律由防砂参数和生产条件控制。使用实际防砂井资料进行了模型应用案例分析,结果和规律与现场实际基本相符。该模型进一步修正后可用于防砂井动态产能预测与评价。     

以新型空壳颗粒复合材料为分配层的成层防护结构抗爆性能试验研究

以成层式人防工程为试验研究背景,分别采用黄沙和新型空壳颗粒复合材料做分配层,开展大比尺集团装药化爆模拟试验,考察其衰减爆炸冲击波的能力。试验结果表明,采用颗粒复合材料构建的分配层,在相同装药量下,遮弹层的变形和破坏明显大于黄沙分配层,应力波上升沿加大,波形脉宽增加,平均峰值应力大概是黄沙分配层的0.65左右,除装药正下方的个别颗粒遭受粉碎性破坏外,绝大部分空壳颗粒依然完好,只是表层泡沫陶瓷略有破损。这说明新型空壳颗粒复合材料不但对爆炸波有非常明显的衰减弥散作用,而且具备承受多次打击的能力。用这种材料做分配层,可大幅提高地下人防工程的抗爆能力。     

颗粒膜的结构及对乳状液稳定性的影响机制

以颗粒膜为基础的乳液因广泛应用于油田、造纸、食品、化妆品、医药等领域而备受关注。本文归纳了在油水界面形成颗粒膜的颗粒密度、粒径及润湿性等特征,并阐述了颗粒在界面吸附扩散的行为。重点从颗粒垂直于界面的分配及方向,颗粒在界面内的排列及进一步形成的空间结构综述颗粒膜结构。总结了颗粒在界面的状态及颗粒膜结构的影响因素,并从能量及力学角度进行分析。颗粒膜对乳状液稳定性影响机制主要体现在颗粒膜的结构及界面黏弹性：颗粒膜的结构阻隔液滴之间的碰撞聚并,这是乳状液稳定的基础;同时,颗粒膜通过改变界面黏弹性使得液滴在运动、碰撞、絮凝时不会轻易崩溃,从而强化了乳状液的稳定性。固体颗粒为乳化剂形成颗粒膜稳定乳状液的机制探究为稳定乳状液制备以及乳状液的破乳提供理论依据,具有现实意义。最后,本文就颗粒膜稳定乳状液的机制研究展望了其未来发展方向及可能的解决途径。     

Tomographic PIV investigation on coherent vortex structures over shark-skin-inspired drag-reducing riblets

Nature has shown us that the microstructure of the skin of fast-swimming sharks in the ocean can reduce the skin friction drag due to the well-known shark-skin effect.In the present study,the effect of shark-skin-inspired riblets on coherent vortex structures in a turbulent boundary layer(TBL) is investigated.This is done by means of tomographic particle image velocimetry(TPIV) measurements in channel fl ws over an acrylic plate of drag-reducing riblets at a friction Reynolds number of 190.The turbulent fl ws over drag-reducing riblets are verifie by a planar time-resolved particle image velocimetry(TRPIV) system initially,and then the TPIV measurements are performed.Two-dimensional(2D) experimental results with a dragreduction rate of around 4.81% are clearly visible over triangle riblets with a peak-to-peak spacing s+of 14,indicating from the drag-reducing performance that the buffer layer within the TBL has thickened;the logarithmic law region has shifted upward and the Reynolds shear stress decreased.A comparison of the spatial topological distributions of the spanwise vorticity of coherent vortex structures extracted at different wall-normal heights through the improved quadrant splitting method shows that riblets weaken the amplitudesof the spanwise vorticity when ejection(Q2) and sweep(Q4) events occur at the near wall,having the greatest effect on Q4 events in particular.The so-called quadrupole statistical model for coherent structures in the whole TBL is verified Meanwhile,their spatial conditional-averaged topological shapes and the spatial scales of quadrupole coherent vortex structures as a whole in the overlying turbulent fl w over riblets are changed,suggesting that the riblets dampen the momentum and energy exchange between the regions of near-wall and outer portion of the TBL by depressing the bursting events(Q2 and Q4),thereby reducing the skin friction drag.     

A symplectic integrator for dynamic coupling between nonlinear vessel motion with variable cross-section and bottom topography and interior shallow-water sloshing

The coupled motion between shallow-water sloshing in a moving vessel with variable cross-section and bottom topography, and the vessel dynamics is considered, with the vessel dynamics restricted to horizontal motion governed by a nonlinear spring. The coupled fluid and vessel equations in Eulerian coordinates are transformed to the Lagrangian particle path setting which leads to a formulation with nice properties for numerical simulation. In the Lagrangian representation, a simple and fast numerical algorithm based on the Störmer–Verlet method, is implemented. The numerical scheme conserves the total energy in the system, as well as giving the partition of energy between the fluid and vessel. Numerical simulations of the coupled nonlinear dynamics are presented.     

基于DIC的受载砂土颗粒体系的细观参数及运动分析

本文结合DIC方法对砂土体系进行可视化的直剪试验。通过对所采集的数字图片进行预处理,得到利于计算砂土体系细观参数的优质图片,经过图像分割和区域标记,最后利用MATLAB计算得到体系孔隙率、定向率和各向异性率等细观参数。对图像进行DIC计算分析,得到砂土颗粒体系的位移等值线云图,分析了体系内的砂土颗粒在直剪试验中的运动情况,并定性描述了位移突变对体系部分力链传递过程的影响。     

基于无迹粒子滤波的飞行器再入段传递对准方法

针对再入式飞行器再入返回前需要精确和快速的传递对准需求以及空间环境特殊性,提出了一种适用于再入式飞行器空间平台对准方法,建立J2000坐标系下光纤陀螺惯组的系统状态方程,利用星敏感器输出的姿态和光纤陀螺惯组输出的姿态差值作为观测量,实现飞行器的"姿态"匹配传递对准算法。针对系统状态方程和量测方程的非线性特性以及单独使用无迹卡尔曼滤波(UKF)方法存在滤波不稳定性和单独使用粒子滤波(PF)方法存在粒子退化导致估计误差变大的问题,提出了无迹粒子滤波(UPF)算法,采用卫星工具包(STK)软件模拟飞行器在轨运行进行仿真验证。仿真结果表明:该方法可以在3s内使得失准角的估计误差小于0.6',估计精度提高了30%,收敛时间缩短了70%,相比传统方法提高了对准的精确性和快速性。     

三维溃坝流的光滑粒子动力学方法模拟

采用光滑粒子动力学SPH(Smoothed Particle Hydrodynamics)方法对三维溃坝流问题进行了数值模拟。为了逼真地模拟出坝内水体与壁面间相互作用而产生的水花飞溅、融合以及近壁面流动等现象,加入了混合长度形式的湍流模型。为了有效地防止粒子穿透固壁,提出了一种新型的适合三维数值模拟的固壁边界处理方法。应用SPH方法对三维溃坝流进行了数值模拟,并分别考虑了未添加障碍物和添加圆柱障碍物两种情形。计算结果表明,改进SPH方法能够精细地捕捉溃坝流在不同时刻的自由液面,并获得稳定而精确的数值结果。     

块石形状对土石混合体力学行为影响的颗粒流模拟

考虑块石形状为球体、正方体和长方体三种情况,通过正方体与球体相比较来探究块石棱角度不同对土石混合体力学特性的影响,通过长方体与正方体相比较来探究块石球度不同对土石混合体力学特性的影响。首先,提出特定形状块石三维离散元精细建模的方法;接着建立含石量为30%和80%的块石形状分别为球体、正方体和长方体的土石混合体三维颗粒流数值模型;然后,对土石混合体大三轴试验进行颗粒流模拟,获得了不同含石量及不同块石形状的土石混合体试样的宏观力学特征;最后,详细分析了块石形状对土石混合体力学行为影响的细观机理。结果表明,块石形状会影响土石混合体的力学行为,其影响的程度与含石量密切相关;配位数、块石颗粒平均旋转量和摩擦功的演化规律能够很好地从细观水平上反映块石形状的作用机理。     

Improvement of boundary conditions for non-planar boundaries represented by polygons with an initial particle arrangement technique

The boundary conditions represented by polygons in moving particle semi-implicit (MPS) method (Koshizuka and Oka, Nuclear Science and Engineering, 1996 Koshizuka, S., and Y. Oka. 1996. “Moving-Particle Semi-Implicit Method for Fragmentation of Incompressible Fluid.” Nuclear Science and Engineering 123: 421–434.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) have been widely used in the industry simulations since it can simply simulate complex geometry with high efficiency. However, the inaccurate particle number density near non-planar wall boundaries dramatically affects the accuracy of simulations. In this paper, we propose an initial boundary particle arrangement technique coupled with the wall weight function method (Zhang et al. Transaction of JSCES, 2015 Zhang, T.G., S. Koshizuka, K. Shibata, K. Murotani, and E. Ishii. 2015. “Improved Wall Weight Function With Polygon Boundary in Moving Particle Semi-Implicit Method.” Transaction of JSCES. No. 20150012. [Google Scholar]) to improve the particle number density near slopes and curved surfaces with boundary conditions represented by polygons in three dimensions. Two uniform grids are utilized in the proposed technique. The grid points in the first uniform grid are used to construct boundary particles, and the second uniform grid stores the same information as in the work by Zhang et al. The wall weight functions of the grid points in the second uniform grid are calculated by newly constructed boundary particles. The wall weight functions of the fluid particles are interpolated from the values stored on the grid points in the second uniform grid. Because boundary particles are located on the polygons, complex geometries can be accurately represented. The proposed method can dramatically improve the particle number density and maintain the high efficiency. The performance of the previously proposed wall weight function (Zhang et al.) with the boundary particle arrangement technique is verified in comparison with the wall weight function without boundary particle arrangement by investigating two example geometries. The simulations of a water tank with a wedge and a complex geometry show the general applicability of the boundary particle arrangement technique to complex geometries and demonstrate its improvement of the wall weight function near the slopes and curved surfaces.