基于三维构造张量的位场边界识别滤波器

地质体构造边界位置的确定是位场数据解释中的一项重要工作,现有很多基于位场梯度张量数据的边界检测滤波器,但存在识别边界位置模糊且无法均衡深浅地质体异常的缺点.本文定义了位场数据的三维构造张量,并提出基于位场构造张量的边界滤波器.为了同时显示不同振幅异常的边界位置,对新定义的滤波器进行归一化处理.在高阶均衡滤波器的计算中需要计算位场的垂向高阶导数,本文引入一种计算的稳定算法,基于拉普拉斯方程利用位场水平导数求解垂向导数,可减小垂向导数计算中产生的误差.将定义的滤波器应用到合成的重磁数据中证明了新方法相比传统的滤波器能更加清晰、准确地圈定边界位置,而且针对同时含有正负异常的地质情况,可避免产生额外的错误边界.最后将新的滤波器应用到实测的重磁数据的解释中,结果显示基于构造张量的滤波器可更准确清晰地划分出断裂的边界位置,发现更多的构造细节.     

The elastic constants of monoclinic single-crystal chrome-diopside to 1,300 K

Values of the complete adiabatic elastic tensor for single-crystal chrome-diopside (a monoclinic pyroxene mineral) are presented from 298 to 1,300 K. The data were obtained using resonant ultrasound spectroscopy (RUS). They are the first published results for the temperature T dependences of the 13 individual elastic constants C _ij of any clinopyroxene mineral. Each C _ij is appropriately described by a linear function in T throughout the range of T. Values for each (∂C _ij /∂T) _P in GPa K⁻¹ are as follows: C ₁₁, −0.0291; C ₂₂, −0.0248; C ₃₃, −0.0179; C ₄₄, −0.0103; C ₅₅, −0.0077; C ₆₆, −0.0152; C ₁₂, −0.0119; C ₁₃, −0.0064; C ₂₃, 0.0000; C ₁₅, 0.0025; C ₂₅, 0.0022; C ₃₅, −0.0046; and C ₄₆, 0.0026. Values of (∂M/∂T) _P in GPa K⁻¹, where M represents an isotropic bulk property calculated from the C _ij data, are as follows: adiabatic bulk modulus K _S , −0.0123; isothermal bulk modulus K _T , −0.0178; and shear modulus G, −0.00998. Some diopside derivatives, notably (∂K _S /∂T) _P , (∂K _T /∂T) _P , and (∂V _P /∂T) _P , where V _P is the compressional wave velocity, have smaller magnitudes than all other minerals of importance in Earth’s mantle, thus, confirming predictions from systematics studies. We find several dimensionless quantities for this monoclinic mineral have normal values compared to other mantle minerals. Further, αK _T (α is the volume coefficient of thermal expansion) for diopside is approximately independent of both T and volume V at elevated temperature, so its equation of state is accurately expressed in simplified form.     

Free surface flow impacting on an elastic structure: Experiment versus numerical simulation

The purpose of this study is to investigate the phenomenon of free surface flow impacting on elastic structures, which is a research topic of great interest in ship and ocean engineering. A series of quasi two-dimensional experiments on dam-break with an elastic plate are conducted. The main features of free surface flow impacting on elastic structures including large impacting force, structural vibration, violent free surface flow, are investigated. The coupled FDM–FEM method developed by the authors is applied for numerical simulation of such dam-break problem. Extensive analysis and discussion based on the comparisons between experimental data and numerical results are made and presented in this paper.     

航空重力梯度数据的模拟研究

针对提高模拟航空重力梯度的精度研究不足的现状,文章运用地球重力场模型构建了计算单点的扰动重力梯度分量的数学模型,推导了Legendre函数及其导数的无奇异性递推计算公式,克服了计算高阶次Legendre函数Pnm(cosθ)的下溢问题,随之对数学模型进行优化;最后利用EGM2008地球位模型模拟生成4km飞行高度处的扰动重力梯度张量作为航空重力梯度观测值。模拟试算表明,该模型能够满足当前航空重力梯度测量的精度要求。     

High pressure elasticity and phase transformation in brucite, Mg(OH)2

The first pressure derivatives of the second-order elastic constants have been calculated for brucite, Mg(OH)₂ from the second- and third-order elastic constants. The deformation theory and finite strain elasticity theory have been used to obtain the second- and third-order elastic constants of Mg(OH)₂ from the strain energy of the lattice. The strain energy ϕ is calculated by taking into account the interactions up to third nearest neighbors in the Mg(OH)₂ lattice. ϕ is then compared with the strain dependent lattice energy from continuum model approximation to obtain the expressions of elastic constants. The complete set of six second-order elastic constants C _IJ of brucite exhibits large anisotropy. Since C ₃₃ (= 21.6 GPa), which corresponds to the strength of the material along the c-axis direction, is less than the longitudinal mode C ₁₁ (= 156.7 GPa), the interlayer binding forces are weaker than the binding forces along the basal plane of Mg(OH)₂. The 14 nonvanishing components of the third-order elastic constants, C _IJK , of brucite have been obtained. All the C _IJK of brucite are negative except the values of C ₁₁₄ (= 230.36 GPa), C ₁₂₄ (= 75.45 GPa) and C ₁₃₄ (= 36.98 GPa). The absolute values of the C _IJK are, in general, one order of magnitude greater than the C _IJ ’s in the Mg(OH)₂ system as usually expected for a crystalline material. To our knowledge, no previous data are available to compare the pressure derivatives of brucite. The pressure derivatives of the two components viz., C ₁₄ and C ₃₃ become negative indicating an elastic instability in brucite while under pressure. This may be related to the phase transition of brucite largely involving rearrangements of H atoms revealed in the Raman spectroscopic, powder neutron diffraction and synchrotron X-ray diffraction studies.     

大地电磁阻抗张量正则分解的地质含义分析

由于常规的张量旋转分析方法的局限性，在处理三维阻抗问题时不能得到令人满意的结果.最新发展的三维阻抗张量的正则分解方法由于不需约束条件，可以自然的从三维阻抗张量中提取完整的参数体系，因而具有很好的应用前景.本文对正则分解的方法从理论上进行了分析和证明，对原有参数计算公式做适当的修改，将正则分解与坐标旋转结合起来，分析二者之间的联系与区别，并探讨参数的地质含义.     

1990年10月天祝—景泰6.2级地震孕震过程剖析及其意义

本文依据天祝—景泰6.2级地震的地震构造、应力场条件及前兆分布,分析研究了该次地震的孕震过程。认为天祝一景泰6.2级地震的力源主要是NE向区域应力场的挤压,垂直差异运动对地震起了促震作用。孕震过程形成的震源应力场在空间上有四象限分布特征,时间上经历了孕育、形成、增强等阶段,相应地震源区则经历了弹性形变、塑性形变、塑性硬化及失稳错动的过程。     

弹性波边界元法正演模拟

弹性波边界元地震模型方法（ＢＥＥＳＭ），实现了二维和三维问题的纵、横波及转换波的同时模拟，并且能模拟任意复杂构造的地震声波正演模型．根据地震模型的特点，本文发展了数值积分计算与矩阵消元同步进行的块状高斯消元法；用解析法处理奇异积分；用无限元法处理边界吸收问题；采用单元长度随介质速度和计算频率变化的变单元算法，及自动剖分单元等技术，提高了计算精度，节省了内存，缩短了计算时间．