Centennial to millennial geomagnetic secular variation

A time-varying spherical harmonic model of the palaeomagnetic field for 0–7 ka is used to investigate large-scale global geomagnetic secular variation on centennial to millennial scales. We study dipole moment evolution over the past 7 kyr, and estimate its rate of change using the Gauss coefficients of degree 1 (dipole coefficients) from the CALS7K.2 field model and by two alternative methods that confirm the robustness of the predicted variations. All methods show substantial dipole moment variation on timescales ranging from centennial to millennial. The dipole moment from CALS7K.2 has the best resolution and is able to resolve the general decrease in dipole moment seen in historical observations since about 1830. The currently observed rate of dipole decay is underestimated by CALS7K.2, but is still not extraordinarily strong in comparison to the rates of change shown by the model over the whole 7 kyr interval. Truly continuous phases of dipole decrease or increase are decadal to centennial in length rather than longer-term features. The general large-scale secular variation shows substantial changes in power in higher spherical harmonic degrees on similar timescales to the dipole. Comparisons are made between statistical variations calculated directly from CALS7K.2 and longer-term palaeosecular variation models: CALS7K.2 has lower overall variance in the dipole and quadrupole terms, but exhibits an imbalance between dispersion in   g ¹₂  and   h ¹₂  , suggestive of long-term non-zonal structure in the secular variations.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.     

大湾铀矿田铀-镭平衡系数及迁移富集特征研究

从分析大湾铀矿田各矿床(点)的铀镭平衡系数变化特征出发,结合放射性同位素比值特征分析铀镭迁移富集行为,指出了各矿床(点)氧化带的深度,以了解各铀矿床(点)潜在的工业价值.     

考虑滑脱效应和温度场的煤层气渗流数值模拟

由于煤层气的解吸热效应,煤层气的运移过程是一个非等温过程。因此,温度场对煤层气渗流有着重要的影响。前人在研究煤层气渗流规律的研究中并没有同时考虑滑脱效应和温度场因素的影响,而在实际的深部开采中综合考虑滑脱效应和温度场的影响对研究深部煤层气运移规律有着重要的意义。因此本文建立了考虑滑脱效应及温度场的煤层气渗流数学模型,利用有限元数值方法研究了考虑滑脱效应和温度场耦合的煤层气运移规律;研究了考虑滑脱效应和温度场对压力场分布的影响;对考虑滑脱效应及温度场因素的煤层气产量进行了预测。得出了随着温度的升高煤储层的压力在降低,温度的升高对煤层气的产量有着负面的影响这一重要结论。     

USB总线与CAN总线协议转换器设计

现场总线(F ield bus)是将自动控制系统中底层的现场控制器和现场智能仪表设备互连的实时控制通讯网络,是5C(Computer,Control,Commun ication,CRT,Change)技术相结合的产物。CAN(Control Area Network)总线是现场总线的一种,广泛应用于汽车行业、机械工业、家用电器、传感器等领域,已经形成国际标准,是被公认的几种最有前途的现场总线之一。提出了一种通过计算机USB接口实现现场总线CAN与计算机(服务器)之间通信的方案。具体介绍了使用USB总线接口芯片CH372和独立的CAN控制器SJA1000实现的硬件电路设计方法、本地端软件的编写方法,给出了主要程序模块的流程。     

隧道锚杆锚固质量无损检测技术

采用声波检测仪,应用声频应力波在不同波阻抗面反射的能量和相位的变化原理,对隧道的锚杆锚固质量进行了无损检测试验研究,并对检测技术进行了探讨.结果表明,作为一种工程质量管理辅助手段,采用应力波对锚杆锚固质量进行无损检测,丰富了隧道围岩锚固质量检测方法,为隧道工程建设提供更好的质量保障.     

吉林台爆破料应力应变关系的三轴试验研究

通过室内大型三轴实验，研究了吉林台水库爆破料在不同级配下的应力应变关系，得到了在一定击实功下爆破料的最大干密度随细料含量变化的规律，分析了爆破料在不同级配和不同围压下应力与应变的变化规律、轴向应变与体积应变的关系及抗剪强度变化特性．从微观的角度说明了变化规律产生的原因，得出爆破料的抗剪强度随级配的变化而变化的规律．     

青藏高原中部的东西向扩张构造运动

系统分析了1933~2003年间青藏高原及其周缘发生的745个中、强地震的震源机制解,研究了高原地壳构造运动及其动力学特征。结果表明,大量正断层型地震集中发生在青藏高原中部海拔4000m以上的地区,其中许多地震是纯正断层型地震。震源机制结果显示,该区正断层型地震的断层走向多为南北方向,断层位错矢量的水平分量均位于近东西方向,这表明青藏高原高海拔地区存在着近东西方向的扩张构造运动。地震震源应力场的研究结果表明,在高原中部高海拔地区,E-W向或WNW-ESE向的水平扩张作用控制着该区的地壳应力场。青藏高原高海拔地区近东西方向的扩张构造运动是该区引张应力场的作用结果,其动力学原因可能与持续隆升的高原自重增大引起的重力崩塌及其周边区域构造应力状况有关。而青藏高原周缘地区,除了东部边缘外,南部的喜马拉雅山前沿以及青藏高原的北部、西部边缘所发生的绝大部分地震都是逆断层型或走滑逆断层型地震。在青藏高原周缘地区,北东或者北北东方向水平挤压的构造应力场为优势应力场。在中国西部的大范围内,主压应力P轴水平分量位于NE-SW方向,形成了一个广域的NE-SW方向的挤压应力场。青藏高原及其周缘应力场特征表明,印度板块的北上运动以及它与欧亚板块之间的碰撞所形成的挤压应力场是高原强烈隆起的直接原因。在青藏高原中南部形成了近东西向引张应力场为主的区域,并以东西向扩张构造运动部分释放其应力积累。研究高原高海拔地区的引张应力场和近东西向扩张构造运动的特征,对于认识青藏高原强烈隆起的地球动力学过程与机制,有着重要的理论意义。     

中国大陆科学钻探主孔声发射法现今地应力状态的确定

对中国大陆科学钻探主孔的岩心进行了声发射测量,确定了301～1531m深度的最大主应力。并与钻孔崩落法(深度1269～1655m)测量结果进行了对比,结果表明,声发射测量所得测值基本上落在钻孔崩落法测值的趋势线上,两种方法所得结果一致,说明测量结果可信。测量结果表明科学钻探主孔地应力大小随深度增加,在浅部301m最大主应力为13.4MPa,在深部1655m为55.2MPa。随深度的增加率为0.0279MPa/m。最大主应力方向为N54°±3.3°E,且方向不随深度变化。