Na+,K+,Rb+,Cs+//SO2-4-H2O体系固液相平衡热力学模拟

锂云母、铯榴石等矿物的硫酸盐焙烧浸出,是从中提取锂铷铯钾等碱金属元素的重要步骤,所得硫酸盐浸出液可近似看作Li⁺,Na⁺,K⁺,Rb⁺,Cs+ // SO₄2--H₂O六元水盐体系;其在浓缩过程呈现复杂的成盐特征,增加了各碱金属元素分离提取的难度,研究其相关水盐体系相图和相平衡性质对理解浸出液中碱金属的分离具有重要的理论指导意义。本文采用Pitzer-Simonson-Clegg超额Gibbs自由能方程和具有内在热力学一致性的CALPHAD方法,构建了Na⁺,K⁺,Rb⁺,Cs⁺//SO₄2--H₂O五元体系的多温固液相平衡热力学模型。借助该模型:1)实现了对Na₂SO₄₊Rb₂SO₄₊H₂O、Na₂SO₄₊Cs₂SO₄₊H₂O和K₂SO₄₊Cs₂SO₄₊H₂O三元体系的多温溶解度等温线的准确描述,并进一步构建了这些体系在240.15~380.15 K温度范围内的多温相图;2)实现了对K₂SO₄+Rb₂SO₄₊H₂O和Rb₂SO₄₊Cs₂SO₄₊H₂O体系中存在的固溶体相—水溶液平衡关系的准确描述,确认了(K,Rb)₂SO₄连续固溶体、以及(Rb,Cs)₂SO₄(α)和(Rb,Cs)₂SO₄(β)间断固溶体存在的合理性;3)预测了上述五元体系中所包含的全部4个四元子体系(Na₂SO₄₊ K₂SO₄₊ Rb₂SO₄+H₂O、Na₂SO₄ + K₂SO₄₊Cs₂SO₄₊ H₂O、Na₂SO₄₊ Rb₂SO₄₊Cs₂SO₄ + H₂O和K₂SO₄₊Rb₂SO₄₊Cs₂SO₄₊ H₂O)的相图,明确了四元相图中固溶体的类型及形成区域。然而,模型预测结果与现有四元体系实验数据间存在普遍的不一致性,尤其是固溶体形成区域,但模型预测结果从非实验角度给出了重要参考,对现有实验提出了挑战,未来对这些体系开展更为精细的实验研究非常必要。     

徐家围子断陷深层火石岭组致密火山岩储层地应力分布规律研究

徐家围子断陷位于松辽盆地北部,其深层中生代火石岭组为致密火山岩气藏,天然裂缝发育,地应力分布非均质性强。结合诱导裂缝法与井径崩落法进行单井现今地应力方向分析,利用声波测井法计算现今地应力大小的纵向分布。依据火石岭组构造顶底面图、火山岩相组及断层分布特征建立非均质三维地质模型;在动、静态岩石物理参数拟合校正的基础上,结合实验测试及已有研究成果,确定不同岩相组和断层岩石物理参数,建立三维力学模型;利用Ansys有限元数值模拟软件建立火石岭组三维数学模型并进行相关运算获得三维现今地应力分布模型。计算结果表明在火石岭组地层中,水平最大主应力方向主要为东西向,应力值范围在86~110 MPa;水平最小主应力方向主要为南北向,应力值范围在67~84 MPa。分析模拟结果可知火山岩相、断层和构造起伏三者对火石岭组现今地应力分布影响较大。其中水平主应力的方向主要受断层和近火山口相分布的影响,而水平主应力的大小则是受三者综合作用。在构造低部位,近火山口相组发育处,断层上盘及断层端部皆为主应力的集中区域。依据现今地应力研究成果可为徐家围子断陷下一步开发井网部署、压裂改造方案和水平井的设计以及注水管理提供重要指导。      

渤中凹陷BZ1探区Qa构造速度异常分析

渤中凹陷BZ1探区Qa构造两口钻井在相距1.3 km的实钻过程中出现了速度异常,三维地震资料时间剖面地层深度与钻井地质分层深度出现了不一致的矛盾现象,即新钻井Qa 3井的钻井预测深度相对于本井地质分层深度更深,速度相对于Qa 2井变小,因而两口井地震速度在横向分布上出现了明显差异。分析表明,速度异常主要由沉积和构造共同作用所引起。利用地震叠加速度谱建模,通过井点速度场、趋势种子点进行校正,获取到更加准确的地下速度,并进行高精度变速成图,能够真实反映地下构造。对比单井速度预测钻井深度与单井速度场校正速度谱建模速度场预测钻井深度,并对比其误差精度,可以看出后者有效地解决了实际钻探过程中因速度异常而导致预测深度误差偏大的问题。因此钻前进行速度场建模与变速构造成图,对今后勘探井或评价井的部署及钻井深度精确预测具有重要的指导作用。     

多点地质统计学建模研究进展

从精细油藏描述中地质建模的意义和现状入手,介绍了多点地质统计学建模研究现状及其与传统地质建模方法的差 异。以辽河盆地西部凹陷某蒸汽驱试验区为例,分析了多点地质统计学建模中建模的基础、训练图像的建立、多点地质统计学 建模与传统地质建模相比所具有的优势等内容。指出多点地质统计学在井间预测方面具有明显优于其他传统建模方法的特 点。在文献调研基础上,结合自身工作实践,探讨了多点地质统计学建模目前存在的问题和未来的发展方向。指出未来多点 地质统计学建模的发展方向主要包括多信息综合地质成因分析基础上的训练图像获取、多点地质统计学算法进行改进和完善 和多点地质统计学建模方法应用领域的扩大等。     

粒度分布的端元建模分析及检验:以“吉兰泰—河套”盆地西部DK-12钻孔晚第四纪沉积物为例*

端元建模分析能够从复杂的多峰分布特征的沉积物中提取出具有不同沉积动力过程的端元,但是,由于沉积物的粒度分布还受到沉积环境等多种因素的影响,该方法的有效性及获得的端元组分的地质意义有待其他环境代用指标的进一步检验。以位于“吉兰泰—河套”盆地西部磴口次级隆起区的DK-12钻孔晚第四纪沉积物为例,采用BEMMA算法对该钻孔沉积物的粒度资料进行了端元建模分析,并以黏土矿物组合和前人的孢粉组合数据作为检验指标,结合该地区的区域地质背景,对获得的4个端元进行了综合检验分析,认为获得的沉积物粒度端元具有明确的地质意义,其中EM 1为远源粉尘、EM 2为近源的风成沙、EM 3和EM 4为河流冲积沙。     

海洋三维VC观测系统优化设计

观测系统设计在地震勘探中起着至关重要的作用,最佳观测系统不仅可以提高资料品质,还能够降低采集成本。通过对海洋垂直缆（vertical cable,VC）进行正演模拟我们发现,随着偏移距的增大,同相轴会出现交叉、合并现象,地层顺序也会发生变化。我们针对单VC和多VC的特点,分别设计了相应的观测系统。然后对所设计的观测系统进行评价与优化,并得出了一些有益的结论：采用增加激发点密度的方法对面元覆盖次数的改善效果要好于采用增加激发面积的方法对面元覆盖次数的改善效果;当目标层存在倾角时,通过在构造走向上增加激发线条数,同时在下倾方向上增加激发线长度可以补偿照明损失;当目标层为背斜或向斜时,通过增大最大环半径来增加背斜和向斜照明范围的方法奏效甚微;当激发面积和接收面积相等时,通过同时增大激发面积和接收面积来提高中心区域面元覆盖次数的方法行不通,而当激发面积大于接收面积时则该问题得以解决。     

Investigation of drag crisis phenomenon using CFD methods

When fluid flow passes a cylinder, the drag crisis phenomenon occurs between the sub-critical and the super-critical Reynolds numbers. The focus of the present studies was on the numerical prediction of the drag crisis based on CFD methods. In this work, block structured meshes with refined grids near the cylinder surface and in the downstream were employed. Both 2D and 3D simulations were performed using various turbulence models, including the SST k − ω model, the k − ϵ model, the SST with LCTM, the DES model, and the LES model. In the convergence studies, the effects of the grid size, the time step, the first grid size and the aspect ratio (for 3D simulations) on the solutions were examined. The errors due to spatial and time discretizations were quantified according to a V&V procedure. Validation studies were carried out for various Reynolds numbers between Re = 6.31 × 10⁴ and 7.57 × 10⁵. The averaged drag force, the RMS of lift force and the Strouhal number were compared with experimental data. The studies indicated that standard 2D and 3D RANS methods were inadequate to capture the drag crisis phenomenon. The LES method however has the potential to address the problem.     

A modified method to estimate chain inverse catenary profile in clay based on chain equation and chain yield envelope

Offshore floating facilities are fixed by anchoring systems embedded in seabed soils through chains or ropes. The chain inverse catenary profile embedded in soils influences both the anchor failure mechanism and the anchor holding capacity. The chain mobilizes varying soil normal and tangential resistances during motion, hence it is with difficulty to depict the chain profile. The present work proposed a modified method to estimate the chain inverse catenary profile with high accuracy based on the chain equations and the chain yield envelope. A testing arrangement with three load cells and two MEMS (Micro-electromechanical systems) accelerometers included was designed in model tests. By model tests, the loading combinations of the soil tangential and normal resistances on the chain were obtained and the yield envelopes for both chain and rope were determined. In addition, supplemental model tests were performed to validate the modified method proposed in this study, and the testing results indicated that the estimated chain inverse catenary profile was in good agreement with the actual one. Moreover, the testing arrangement is beneficial in investigating the chain-soil-anchor interaction.     

江苏近海风暴潮增水数值模拟研究

为研究江苏近海海域风暴潮的特性以及为该海域风暴潮增水变化机理及后报做铺垫,本文基于FVCOM(Finite Volume Coast and Ocean Model)海洋模式和Jelesnianski圆形台风风场模型,建立了江苏近海风暴潮数值模型,并对江苏近海的天文潮以及1109号台风和1210号台风引起的风暴潮进行模拟。结合验潮站水位观测,研究了连云港站和吕泗站的天文潮和风暴潮增水过程。我们将风暴潮与天文潮非线性作用下的风暴潮增水和纯风暴潮增水过程进行对比,讨论了天文潮与1109号和1210号台风风暴潮之间的非线性作用引起的增水特征。结果均表明,在天文潮高潮时,天文潮和风暴潮之间的非线性作用可以抑制增水,在天文潮低潮时,天文潮和风暴潮之间的非线性作用有利于增水。除了气象因子以及天文潮和风暴潮之间的非线性作用外,该海区的地理环境也对台风风暴潮增水产生影响。因此对江苏近海的海岸线变化和浅滩地形变化进行敏感性试验,结果表明,本文所设计的海岸线变化对该海域的风暴潮增水影响较小,江苏沿海岸线的向外推移使得江苏海域风暴潮的增水略微上涨,而本文所设计的地形的变化对风暴潮增水影响较大。     

Application of smoothed particle hydrodynamics for modeling the wave-moored floating breakwater interaction

The application of a Smoothed Particle Hydrodynamics (SPH) model to simulate the nonlinear interaction between waves and a moored floating breakwater is presented. The main aim is to predict and validate the response of the moored floating structure under the action of periodic waves. The Euler equations together with an artificial viscosity are used as the governing equations to describe the flow field. The motion of the moored floating body is described using the Newton’s second law of motion. The interactions between the waves and structures are modeled by setting a series of SPH particles on the boundary of the structure. The hydrodynamic forces acting on the floating body are evaluated by summing up the interacting forces on the boundary particles from the neighboring fluid particles. The water surface elevations, the movements of the floating body and the moored forces are all calculated and compared with the available experimental data. Good agreements are obtained for the dynamic response and hydrodynamic performance of the floating body. The numerical results of different immersion depths of the floating body are compared with that of the corresponding fixed body. The effects of the relative length and the density of the structure on the performance of the floating body are analyzed.