Evaluation of one-dimensional in situ leaching processes

A methodology to characterizing processes of in situ leaching is developed to investigate various parameters that may affect the recovery of a valuable mineral and assure a successful application of the in situ leaching technology. The governing equations of in situ leaching processes for the one-dimensional case are solved both analytically and numerically for both the consumption rate of a lixiviant and the production rate of a target mineral. A numerical simulator, which is developed to evaluate coupled effects among leaching kinetics, solution flow and transport of the dissolved mineral species in saturated ore deposits, is validated against the steady-state solutions and applied to investigate the transient effects of various parameters on the mineral recovery. Results from the evaluation indicate that there exists an optimal flow velocity range of leach solution for the effective leaching of a particular ore deposit. The determination of this optimal leaching velocity may become a key to the design of a real in situ leach mine. Results of the parameter sensitivity study illustrate the relative importance of other parameters such as porosity, ore grade, and reaction rate. Copyright © 1999 John Wiley & Sons, Ltd.     

Solutions for the deformations and stability of elastoplastic hollow cylinders subjected to boundary pressures

In this paper, a closed-form solution is presented for the stress and displacement distributions throughout a hollow cylinder subjected to uniform pressures acting on its internal and external boundary surfaces under plane strain conditions. The material is assumed to be elastoplastic, obeying a Mohr–Coulomb failure criterion, and exhibiting dilatant plastic deformation according to a non-associated flow rule. The newly developed analytical solution is verified through comparison with the solutions obtained from an infinite boundary problem (for which a closed-form solution exists), and numerical analyses using the program FLAC. The solution is also compared with the results of a borehole collapse test on a thick-walled hollow cylinder of synthetic shale. The analytical solution can be used to calculate the stress and displacement distributions around boreholes and other cylindrical cavities under both infinite and finite boundary conditions under both drained and undrained conditions. Copyright © 1999 John Wiley & Sons, Ltd.     

A semi‐distributed parallel‐type linear reservoir rainfall‐runoff model and its application in Taiwan

The main purpose of this paper is to introduce a semi‐distributed parallel surface rainfall‐runoff conceptual model. In this paper, a general solution of the instantaneous unit hydrograph (IUH) has been derived successfully for N linearly connected reservoirs, each having a different storage constant. The solution is a function of geomorphologic parameters, meteorologic factors and roughness coefficients. The model also takes into account the hydrologic response which is influenced by outflow downstream of a reservoir. For calibration, the shuffled complex evolution (SCE) algorithm is used to search for the global optimal parameters of the model. Because of the parallel structure, the mean roughness parameter of the channel becomes a “conceptual” parameter without a real physical meaning. To evaluate the adaptability of the model adopted, three watersheds around the city of Taipei in Taiwan were chosen to test the effectiveness of the model. The study provides an appropriate rainfall‐runoff model for planning flood mitigation in Taiwan. Copyright © 1999 John Wiley & Sons, Ltd.     

Elastically supported cantilever beam subjected to nonstationary seismic excitation

The dynamic behaviour of bridge piers under seismic load is studied here in the context of random vibration. The earthquake excitation is modelled as white noise filtered by the Clough–Penzien filter in cascade with modulation accounting for intensity non-stationarity. The bridge pier modelled as an elastically supported cantilever beam witha lumped mass at the top. An analytical solution is presented for the response statistics, which may be used to develop probabilistic seismic response spectra for design. It is found that the first two modes of the pier approach to rigid-body motion when the stiffness of the elastic support decreases. Seismic responses increase with the top mass, resultingin significantly high displacement and shear but negligible moment at the top, and higher shear and moment at thebase. Lower stiffness of the elastic support increases the pier top displacement and moment responses, but may increase or reduce shear responses. The probabilistic spectrum of the relative displacement between the bridge superstructure and the pier top may depend on the two systems’ relative modal properties. © 1998 John Wiley & Sons, Ltd.     

渤海湾盆地西部凹陷南段成岩相分析与优质储层预测

储层实测孔隙度、普通薄片、铸体薄片和扫描电镜的资料表明,渤海湾盆地辽河坳陷西部凹陷南段新生界储层在纵向上发育三个次生孔隙带。从而在一个普遍低孔渗背景下,形成了孔隙度相对较高的优质储层。优质储层的形成和分布主要受沉积微相和成岩作用的影响与控制。他们形成于河口坝、辫状分流河道以及心滩沉积微相,目前处于早成岩阶段B期—中成岩阶段A1亚期,发育溶蚀成岩相。该文通过模拟古地温、镜质组反射率、甾烷异构化率和自生石英含量随时间的变化规律,预测了西部凹陷南段成岩阶段和成岩相的横向展布。通过成岩相图和沉积相图的叠合,预测了沙三中亚段优质储层的分布。预测结果表明,优质储层主要发育于斜坡区的中部。     

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。分析模拟结果可知火山岩相、断层和构造起伏三者对火石岭组现今地应力分布影响较大。其中水平主应力的方向主要受断层和近火山口相分布的影响,而水平主应力的大小则是受三者综合作用。在构造低部位,近火山口相组发育处,断层上盘及断层端部皆为主应力的集中区域。依据现今地应力研究成果可为徐家围子断陷下一步开发井网部署、压裂改造方案和水平井的设计以及注水管理提供重要指导。      

盐溶液对膨润土膨胀性的影响

膨润土因其优良的膨胀性和低渗透性常被用作深层地质处置库的缓冲和回填材料,但其膨胀性会受到地下流体溶液浓度的影响。针对这种现象,利用单轴固结仪研究了压实商用膨润土在不同上覆荷载和不同浓度NaCl溶液作用下的膨胀变形特性。结果表明：随着NaCl溶液浓度升高,该膨润土的膨胀性显著降低,且最大膨胀应变与上覆荷载之间在双对数坐标系下存在良好的线性关系。利用有效压力的概念和计算方法,实现了用唯一一条曲线来描述不同盐溶液浓度作用下膨润土的膨胀变形行为,同时证明了有效压力理论的合理性。     

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

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