矿体三维数据模型及品位插值方法研究

传统的储量计算方法采用几何平均法计算矿体的体积,用部分化验数据的平均品位代替矿块的整体品位,其计算精度很难满足需要.以三维矿体模型为基础计算矿体体积能够提高体积计算精度,以空间插值技术对矿体任意空间位置的品位进行估值,有助于提高品位估值的精度.提出了一种适于矿体储量计算的混合数据模型（线框块段模型）,详细论述了根据剖面数据生成矿体线框模型进而在线框模型基础上生成块段模型的方法.介绍了距离反比法和克立格法的原理,利用距离反比法结合某金矿体数据进行了储量估算,结果表明以线框块段模型为基础,采用空间插值技术计算矿体储量具有一定的实用价值.     

软土蠕变固结特性研究

通过大量对比试验系统研究了广东典型软土的蠕变固结特性和计算模型。试验结果表明：软土的固结和次固结在变形中所占比例和发生受多种因素影响,包括初始荷载、加荷比、土体固结度、含水率、孔隙比;次固结系数与土体的应力历史、当前固结状态及加荷比有关,在预压荷载下会明显减小;次固结系数与压缩指数之比一般为常数,其值在0.025?0.10之间;剪切蠕变特性与土体所承受的固结压力有关,土样蠕变规律呈应变硬化特征;三轴蠕变特性与排水条件关系密切,固结作用可以削弱蠕变效应;同等加荷条件下,偏应力-轴向应变关系在排水剪切中呈线性关系。而在不排水剪切中呈非线性,并有明显的屈服特征;不排水剪切条件下持续加荷的加荷方式可以适当提高不排水强度。根据上述试验成果分别从考虑非线性本构关系和考虑排水条件分析了蠕变固结模型的建立路径。     

液化土层上建筑物倾斜计算要点初步分析

地震中饱和砂土地基液化会引起结构物倾斜并导致其功能丧失,但目前缺乏相应的数值模拟方法。通过振动台简单模型试验,寻找输入波-基底竖向动应力-基底孔压-结构震陷之间的关系,提出发展液化土层上建筑物倾斜数值模拟方法所需要考虑的要点和应满足的条件。结果表明：（1）分析方法中孔压增长模型应适于不规则波计算,能准确地计算出孔压增长过程,准确地计算出峰值一样但不同波形下孔压增长的差别;（2）孔压增长模型应能反映土的各向异性特性对孔压增长的影响,能合理地计算出拉、压不同应力作用下孔压增长的差异;（3）孔压增长模型应能反映非均等固结条件对孔压增长过程的影响,能合理地计算出结构底部不同固结比土体中孔压增长过程;（4）分析方法中土体变形的计算应能跟踪液化过程中的变形发展,且具备大变形计算能力。     

一种真空预压区地下水位计算方法

根据抽真空时加固区内真空渗流场的形成机制和变化特点,把真空渗流场的形成过程分为非稳定渗流和稳定渗流两个阶段。利用渗流量计算公式,计算加固区内渗流量随时间变化过程曲线。结合两个阶段加固区地下渗流场的变化特性,分别利用Theis和Dupuit计算模型计算不同时刻的地下水位,并用实测进行验证。     

曲线拟合在核探测器信号幅度提取中的应用

研究数字化核探测器输出脉冲信号的幅度提取问题。根据闪烁探测器的电路特性，推导了其脉冲信号理想曲线方程。综合利用三和值法、最小二乘法、试探法等方法，设计了该理想曲线拟合算法。由脉冲信号采样点数据，通过曲线拟合，计算出曲线的极大值，该值即为闪烁探测器输出脉冲信号的幅度。正反演数据结果对比表明，采用曲线拟合算法提取闪烁探测器输出脉冲信号幅度，其方案可行，计算结果更准确。相较于数字滤波和数字脉冲成形方法提取核探测器输出脉冲幅度，该算法更简单，对仪器的硬件处理速度要求更低。该方法开辟了核探测器脉冲幅度提取的新途径。     

热力学计算模拟对初始月幔结构的约束SCIEI北大核心CSCD

月球岩浆洋结晶形成的初始月球内部结构是其后续演化过程的开端,其结晶过程受月球岩浆洋的初始深度和物质组成这两个参数的制约。由于缺少直接来自月球深部的岩石样品,目前关于月球岩浆洋演化过程的探讨主要依赖实验和计算模拟手段。岩浆洋模型中形成的月壳厚度是否与探测结果一致是月球岩浆洋演化模型合理性的重要约束。最新的GRAIL(Gravity Recovery and Interior Laboratory)探测数据推算月壳厚度为34~43km,低于阿波罗时期认为的约70km,这对已有的月球岩浆洋演化模型提出了挑战。本文采用并修正FXMOTR程序包,针对月球岩浆洋在不同的初始深度和物质组成情况下的结晶过程,进行了一系列热力学计算模拟。通过量化月球岩浆洋的初始深度和物质组成对月壳厚度的影响,结合关于月球内部微量元素分配的研究结果,对比了月球岩浆洋结晶后期的残余熔体与原始克里普组分(urKREEP)的成分。本文的模拟结果显示,一个全月幔熔融且初始成分为月球初始上月幔组成(LPUM)的岩浆洋将在其深部结晶2.5%石榴子石,形成的月壳厚度符合GRAIL的约束,并且结晶出了合适的urKREEP成分。在此模型的基础上获取了月球初始的内部成分和密度结构,并对后期月幔翻转(Overturn)的程度进行了探讨。     

胶东金矿成矿深度的构造校正测算及成矿预测

成矿深度测算对于矿床学理论研究和深部找矿都有重要意义。经典的成矿深度“压力/比重”计算方法,缺乏考虑构造应力在成矿过程中的影响。前人按照“压力/比重”的计算方法,提出胶东蚀变岩金矿是6 000~8 000 m深的元古宙成矿,石英脉型金矿是深度在3 000 m左右的中生代成矿,并据此建立了金矿垂直五层楼的分带模式。依据该模式指导的深部预测勘查效果不好。“成矿深度构造校正测算”是近几十年逐渐成长起来的一个新方法,即先减去“构造附加压力”后再进行成矿深度测算。本文介绍“成矿深度构造校正测算”的计算方法,指出其应用条件和预测意义。开展成矿深度构造校正测算需要以下条件:(1)确定成矿模式;(2)开展野外构造变形岩相测量;(3)测量岩石矿物应变,恢复成矿构造应力场,计算构造附加静水压力;(4)测算成矿深度。根据“成矿深度构造校正测算”方法已经获得胶东多个典型金矿成矿深度的测算结果:(1)夏甸金矿成矿深度为-1 979.51~-3 014.72 m;(2)焦家金矿成矿深度为-1 632.4~-2 331.6 m;(3)大尹格庄金矿成矿深度为-2 775.4~-4 164.5 m;(4)新城金矿成矿深度为-1 781.29~-2 750.0 m;(5)玲珑金矿成矿深度为-720.55~-3 454.97 m。根据以上典型金矿成矿深度的测算结果,本文认为胶东金矿属于深-1 000~-4 500 m的浅成热液剪切带型矿床,由此推断胶东典型金矿矿体主要部分仍然赋存在深部。按照“构造校正测算”方法得到的成矿深度,结合地质、物探和化探信息,预测金矿发育“深部第二富集带”,已经得到胶东金矿勘查工作的证实。     

以高斯理论为基础计算弯曲钻孔节点空间坐标

为了拟合出更加接近钻孔实际的曲线,结合钻孔测斜实际,以严密的高斯理论为基础,给出了以钻孔天顶角和倾角为参数的弯曲钻孔的参数方程。并在此基础上,推出了钻孔曲线各点的计算公式,更准确地描述了钻孔的实际情况,为后续工作剖面图的绘制和储量计算提供了基础。这里使用高斯和折线法作为一个实例,对比计算了赤土店铅锌矿区ZK13002钻孔。计算结果表明:在钻孔曲线形状、累计弧长、插值点特征等多个方面的高斯法都优于折线法,而因算法的改变引起的工作量又是很小的。     

Hydrothermal alteration processes in the giant Dahutang tungsten deposit,South China: Implications from litho-geochemistry and mass balance calculation

The giant Dahutang tungsten (W) deposit has a total reserve of more than 1.31 Mt WO₃. Veinlet-disseminated scheelite and vein type wolframite mineralization are developed in this deposit, which are related to Late Mesozoic biotite granite. Four major types of alterations, which include albitization, potassic-alteration, and greisenization, and overprinted silicification developed in contact zone. The mass balance calculate of the four alteration types were used to further understanding of the mineralization process. The fresh porphyritic biotite granite has high Nb, Ta, and W, but low Ca and Sr while the Jiuling granodiorite has high Ca and Sr, but low Nb, Ta, and W concentrations. The altered porphyritic biotite granite indicated that the Nb, Ta, and W were leached out from the fresh porphyritic biotite granite, especially by sodic alteration. The low Ca and Sr contents of the altered Neoproterozoic Jiuling granodiorite indicate that Ca and Sr had been leached out from the fresh granodiorite by the fluid from Mesozoic porphyritic biotite granites. The metal W of the Dahutang deposit was mainly derived from the fluid exsolution from the melt and alteration of W-bearing granites. This study of alteration presents a new hydrothermal circulation model to understand tungsten mineralization in the Dahutang deposit.     

Heat capacity and phase equilibria of wadeite-type K2Si4O9

The low-temperature heat capacity (C _p) of Si-wadeite (K₂Si₄O₉) synthesized with a piston cylinder device was measured over the range of 5–303 K using the heat capacity option of a physical properties measurement system. The entropy of Si-wadeite at standard temperature and pressure calculated from the measured heat capacity data is 253.8 ± 0.6 J mol⁻¹ K⁻¹, which is considerably larger than some of the previous estimated values. The calculated phase transition boundaries in the system K₂O–Al₂O₃–SiO₂ are generally consistent with previous experimental results. Together with our calculated phase boundaries, seven multi-anvil experiments at 1,400 K and 6.0–7.7 GPa suggest that no equilibrium stability field of kalsilite + coesite intervenes between the stability field of sanidine and that of coesite + kyanite + Si-wadeite, in contrast to previous predictions. First-order approximations were undertaken to calculate the phase diagram in the system K₂Si₄O₉ at lower pressure and temperature. Large discrepancies were shown between the calculated diagram compared with previously published versions, suggesting that further experimental or/and calorimetric work is needed to better constrain the low-pressure phase relations of the K₂Si₄O₉ polymorphs. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.