基于单能窄束伽马射线算法计算高分辨率屏蔽自然伽马探测器的三维测井响应

本文应用单能窄束伽马射线理论建立高分辨率屏蔽自然伽马探测器三维测井响应的数值模拟算法.首先利用单能窄束原理,用稳态扩散方程描述伽马光子密度的空间分布,根据扩散方程基本解、放射源空间分布、探测器位置和屏蔽情况,将屏蔽探测器上的总伽马通量表示成放射性地层中的有效探测区域上的体积分和晶体表面上有效接收面的面积分形式.并根据伽马射线传播路径和探测器屏蔽情况,给出有效探测区域的解析表达式,通过数值积分法计算任意复杂情况下探测器上的自然伽马通量,获得伽马测井响应的3D数值模拟算法.并通过数值模拟结果与模型井数据的对比验证了该算法的有效性.最后通过3D数值模拟算法系统研究考察晶体形状（圆柱形晶体、方形晶体）、晶体长度、仪器在井轴中的位置（居中屏蔽探测器、贴井壁屏蔽探测器）、以及测速等不同情况下自然伽马测井响应,设计出新型高分辨率自然伽马测井仪器.     

辅助微分方程完全匹配层在声波方程数值模拟中的应用

在地震波数值模拟中,为提高算法精度,需要使用高阶时间更新格式,而普通的非分裂完全匹配层（PML）吸收边界局限于低阶时间格式。辅助微分方程完全匹配层（ADE PML）是一种可以适应任意阶时间格式的非分裂完全匹配层技术,且可以直接应用复频移拉伸算子以提高PML在高角度入射时的效果。作者将ADE PML应用于声波方程四阶Runge Kutta时间格式的数值模拟中,对其吸收效能进行了检验。数值模拟表明,复频移ADE PML在高角度入射时表现优于非复频移ADE PML。另外,不同辅助变量更新格式的吸收效果存在微小差异,显格式下计算结果与解析解吻合较好。长时间能量衰减计算表明ADE PML可以稳定至2 × 10⁵时间步。     

南黄海中部海底沉积物声阻抗特性

基于在南黄海中部海域获得的大量海底沉积物样品及其测试数据,分析讨论了该区域海底浅表层沉积物声阻抗特性及其与沉积物物理力学性质参数的相关关系。研究表明:研究区浅表层沉积物声阻抗呈现出分区域分布的特征,可分为西南部声阻抗高值区、东北部声阻抗低值区以及北西-东南向声阻抗梯度带3个区域。通过回归分析,建立了声阻抗与沉积物物理力学参数的经验回归方程。结果表明:声阻抗与密度、含水率、孔隙比、孔隙度、液限、塑限、塑性指数、压缩系数等参数的相关性较好,判定系数R²一般大于0.8,而与砂粒含量、黏粒含量、中值粒径、抗剪强度等参数的相关性相对较差,判定系数R²一般小于0.7。将声阻抗与物理力学性质的相关性和声速与物理力学性质的相关性进行了对比,结果表明,声阻抗与上述参数的判定系数R²均大于声速与上述参数的判定系数,从而表明声阻抗是一个能够更好地反映声学特性与物理力学性质关系的沉积物声学特性参数。     

CO_2 Density-Raman Shift Relation Derived from Synthetic Inclusions in Fused Silica Capillaries and Its Application

Abstract: The densities of CO2 inclusions in minerals are commonly used to determine the crystallizing conditions of the host minerals. However, conventional microthermometry is difficult to apply for inclusions of small size (< 5–10 μm) or low density. Raman analysis is an alternative method for determining CO2 density, provided that the CO2 density–Raman shift relation is known. This study aims to establish this CO2 density–Raman shift relation by using CO2 inclusions synthesized in fused silica capillaries. By using this newly-developed synthetic technique, we formed pure CO2 inclusions, and their densities were determined by microthermometry. The Raman analysis showed that the relation between CO2 density (D in g/cm3) and the separations (Δ in cm?1) between the two main bands (i.e. Fermi diad bands) in CO2 Raman spectra can be represented by a cubic equation: D (g/cm3)=0.74203(?0.019Δ3+5.90332Δ2?610.79472Δ+21050.30165)?3.54278 (r2=0.99920). Our calculated D value for a given Δ is between those obtained from two previously-reported equations, which were derived from different experimental methods. An example was given in this study to demonstrate that the densities of natural CO2 inclusions that could not be derived from microthermometry could be determined by using our method.     

三棱柱形网格下自由表面流的三维数值模拟

为了模拟自由表面流动问题,建立了基于三棱柱形网格的三维Navier-Stokes方程的半隐式动压离散模型。出流边界采用在动量方程右端附加线性衰减项的海绵层处理方式来消减波能。引入k-ε双方程紊流模型求解涡黏性系数,使方程达到封闭。用规则周期波通过淹没障碍物的传播变形和复合明渠流动算例,验证了当遇到短波高频问题时,所述模型模拟结果与测量值吻合良好,体现了模型的精确性和有效性。     

The relationship between nonconservative schemes and initial values of nonlinear evolution equations

For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensionalshallow water wave equation as an example, the necessary conditions of computational stability are given.Based on numerical tests, the relationship between the nonlinear computational stability and the construc-tion of difference schemes, as well as the form of initial values, is further discussed. It is proved throughboth theoretical analysis and numerical tests that if the construction of difference schemes is definite, thecomputational stability of nonconservative schemes is decided by the form of initial values.     

利用GPS监测电离层不均匀结构探讨

利用上海地区GPS综合应用网提供的高时空分辨率的双频GPS观测资料,研究了该区域内一电离层不均匀体的产生、消亡过程.首先,采用Kalman滤波的方法改善双频伪距之差的观测精度,并利用参数估计的方法计算该时段内相应的硬件延迟.再根据电离层单层模型,利用GPS双频观测量、测站位置和GPS精密星历,求出GPS信号穿刺点的坐标和垂直方向电离层的电子含量;然后内插并获取其等值线图.等值线图随时间的变化表明,受等离子体湍流的影响,2003年9月8日北京时间9时40分左右在38°N、118°E左右产生了一电离层不均匀体,其尺度大约在50km左右,生存时间大约为5min.受地球重力场和高空风场的影响,该不均匀体向东北方向扩散.然后,利用大气扩散模型,按扩散方程计算分析了该不均匀体可能发生的电离层层区.理论计算表明,该不均匀体发生在电离层扩展F区,高度在350km左右.     

非结构化网格下大范围波生流模拟和应用

波浪破碎引起的沿岸流是近岸海域的关键水动力因素。利用基于缓坡方程得到的光程函数方程和波作用守恒方程建立了考虑绕射效应的大范围波浪传播模型,模型可以考虑流场的影响;将波浪模型计算得到的辐射应力、波浪紊动系数等参数添加到三维水动力模型中,得到大范围近岸波生流的计算模型。模型中流场和波浪可以共用计算网格,且可同步耦合计算;模型基于非结构化网格,可以拟合复杂岸线的变化。模型对波生沿岸流、环流和逆流进行了验证,同时对实际海域的波生流进行了计算,结果表明:该模型对近岸波浪破碎引起的波生流具有很好的精度和适用性,可用于实际工程的计算。     

相混合计算法确定蒙脱石含量的对比研究

基于X射线粉晶衍射分析、化学成分分析数据,依据物质平衡原理,应用相混合计算法对陕西洋县、河南确山、内蒙古赤峰三地的膨润土进行了定量物相分析。根据电子探针微区分析所得蒙脱石晶粒化学成分数据,采用以阴离子为基准的氢当量法,计算了膨润土中蒙脱石的晶体化学式。应用CELLSR软件,计算得到膨润土中各主要矿物的晶胞参数。通过与传统的确定膨润土中蒙脱石含量的吸蓝量法和X射线衍射多相Rietveld分析法进行定量对比,结果表明,相混合计算法是确定多种矿物共生的膨润土中蒙脱石含量最有效的物相分析方法。尽管受不同分析方法的误差积累影响,但是由该法确定蒙脱石的含量仍能保证±20％的计算精度。     

Calculation of seismic transmission speed of HTI medium

During the seismic wave propagation process, as for the anisotropic lower medium, the speed is a function of the propagating direction. This article focuses on solving the problem how to get the transmittance angle and speed, knowing the upper seismic wave propagation velocity and the angle of incidence conditions. The main theories used Snell law, Christoffel equation and speed characteristics. Taking the HTI medium as an example, the authors give the detailed solving process and draw the space velocity characteristic curve.