地震信号广义时频分析及其数值实现

基于分数阶Fourier变换的信号时频分析称为广义时频分析,将其引入地震信号处理中,是对现有基于传统变换时频分析方法的广义化拓展。这里着重数值实现,探索基于分数阶Fourier变换的时频分析,并应用了其在物探技术中新思路和比较的优势。在分数阶Fourier变换基本定义的基础上,数值实现离散分数阶Fourier变换;研究三类广义时频分析方法及其数值实现;针对地震信号,以伪Wigner分布为例,展示广义时频分析效果。经研究表明:实现的数值算法正确有效;地震信号新的广义时频分布能量更加集中,时频分析的聚焦度和信息特征更加清晰。由此可见,广义时频分析为地震信息特征分析、预测和提取提供新的视角和工具,值得进一步研究和推广。     

Interface error analysis for numerical wave propagation

The numerical error associated with finite-difference simulation of wave propagation in discontinuous media consists of two components. The first component is a higher-order error that leads to grid dispersion; it can be controlled by higher-order methods. The second component results from misalignment between numerical grids and material interfaces. We provide an explicit estimate of the interface misalignment error for the second order in time and space staggered finite-difference scheme applied to the acoustic wave equation. Our analysis, confirmed by numerical experiments, demonstrates that the interface error results in a first-order time shift proportional to the distance between the interface and computational grids. A 2D experiment shows that the interface error cannot be suppressed by higher-order methods and indicates that our 1D analysis gives a good prediction about the behavior of the numerical solution in higher dimensions.      

微区-微量样品Rb-Sr同位素分析技术及其应用前景

利用微钻取样技术和微量样品Rb-Sr同位素分析方法,本文对出露在东秦岭造山带的中生代合峪花岗岩的自形钾长石巨晶进行微区-微量样品Rb-Sr同位素组成分析。分析结果表明,钾长石斑晶具有显著的Rb/Sr比值和Sr同位素组成变化,斑晶和基质钾长石均构成年龄为132~133Ma的Rb-Sr等时线,代表岩浆的后期冷却时代。钾长石晶体的初始⁸⁷Sr/⁸⁶Sr比值由边缘相到中心相没有明显的变化,代表花岗质岩浆结晶阶段的Sr同位素组成,暗示合峪花岗岩的钾长石巨晶为原生成因。以高空间分辨率为特征,微区取样技术已经广泛地应用在变质岩和深成岩浆岩的同位素年代学和成因研究。结合微量样品同位素分析技术,微区-微量样品Rb-Sr同位素方法有望在火山岩的成因和年代学方面得到应用。     

In situ production of terrestrial cosmogenic helium and some applications to geochronology

The concentrations of cosmogenic ³He have been measured in a series of basaltic drill core samples from Hawaiian volcanoes Haleakala and Mauna Loa. The ³He concentration in the surface of a radiocarbon dated Mauna Loa flow (20,000 years) gives reasonable agreement with a theoretical production rate of 140 atoms g−1yr−1 and suggests that the uncertainty in this rate is roughly 10%. The results illustrate the feasibility of using ³He to measure exposure ages of young basaltic lava flows and for measuring erosion rates. Erosion rates calculated from the three Haleakala cores range from 7 to 11 meters/million years. The drill core data demonstrate that accurate depth control is crucial to the use and evaluation of cosmogenic helium. Depth profiles from several of the older cores display a non-exponential depth dependence of ³He_c below 170 g cm⁻², which is attributed to the contribution from ⁶Li(n, )T, where the neutrons are from stopped muons. This has important implications for depth dependence of cosmogenic ³He because muons are weakly attenuated compared to the nucleonic component that produces spallation.     

广义混合律及其在地球材料流变学中的应用

本文总结了广义混合律的基本特征，并以具弱相支撑结构的固液相悬浮液（岩浆）和具强相支撑结构的多孔材料（玄武岩、砂岩或部分熔融岩石）为例，系统地介绍了广义混合律在预测多矿物岩石或复合材料弹性力学或流变学性质随组分性质、体积分数和显微构造而变化的应用。广义混合律与现存理论的紧密吻合以及对多矿物岩石和复合材料力学性质描述的准确性绝非偶然，隐藏其后必有尚未人知的物理基础。严格的理论证明极可能是数学力学上富于挑战性的课题。     

河津凝灰岩广义应力松弛特性及数值模拟研究

在很多地下工程实践中,岩体既不是纯蠕变,也不是纯应力松弛,而是随着时间的增加,应力和应变同时发生变化,表现出时间依存性,导致岩体最终破坏,这种现象用一般的蠕变和应力松弛很难解释清楚。全面系统地阐述了广义应力松弛特性及其工程物理意义,借用非线性Maxwell模型,求解了广义应力松弛条件下可变模量本构方程解析解,特别提出了通过交替荷载速率试验和加、卸载试验,精确求解峰前区域荷载速率依从性系数n_1的新方法。进行了70%和90%应力水平下广义应力松弛试验,并讨论和分析了不同应力水平、不同方向系数下广义应力松弛试验的时间依存性,最后用基于非线性Maxwell模型的可变模量本构方程对70%和90%应力水平下的广义应力松弛试验进行了数值模拟,结果表明,计算和试验结果一致性较好,很好地解释了河津凝灰岩广义应力松弛特性。     

油气成藏过程中的同位素测年方法评述

笔者对国内外同位素测定年在确定油气成藏方面的应用现状进行评述的基础上,重点分析了K-Ar法,Ar-Ar法,U-Pb法,Rb-Sr法等方法所存在的问题,认为自生伊利石⁴⁰Ar-³⁹Ar法在不断克服技术难题之后,仍是油气成藏年代学不可缺少的重要手段;激光微区探针⁴⁰Ar-³⁹Ar法与流体包裹体技术相结合,是油气多期成藏⁴⁰Ar-³⁹Ar年代学研究取得成功的突破点;利用沥青和干酪根中微量金属U-Pb,Pb-Pb,Rb-Sr,Sm-Nd和Re-Os体系的同位素分析方法获得油气生成、运移的年龄,是成熟的放射性同位素方法在石油地质学中的应用;晶洞充填物U-Pb测年,可能是碳酸盐岩储层油气成藏年代学研究又一新的突破点。     

碳酸盐矿物激光原位U-Pb定年基本原理、分析方法与地学应用

碳酸盐矿物作为地壳中沉积作用与流体活动的直接产物,记录了沉积、热事件的全过程,是进行同位素年代学研究的理想矿物。碳酸盐矿物传统的U系列定年法、同位素稀释法定年体系成功率低,耗时长,导致定年难度大,限制了碳酸盐矿物地质年代学的发展。近年来,碳酸盐矿物激光原位U-Pb定年技术取得重要进展,使低U碳酸盐矿物地质年代的精确测定成为了可能。碳酸盐矿物LA-ICP-MS U-Pb年代学具有空间分辨率高、耗时短等显著优势,被广泛应用于地学研究中。文章归纳总结了碳酸盐矿物同位素定年体系的基本原理及分析方法,回顾了近5年碳酸盐矿物LA-ICP-MS U-Pb年代学应用的研究进展,重点论述该方法在确定脆性变形时代、岩石破裂以及盆地流体流动时限、碳酸盐地层的成岩时代和成矿热液活动时限4个方面的应用,并梳理了该方法尝试解决的科学问题和取得的新认识。碳酸盐矿物定年方法的发展和应用,增强了人们对地壳变形和演化的理解,解决了部分迄今为止难以确定的地质体历史演化问题,为未来地质年代学的研究提供了新思路,在地球科学研究中具有重要潜力。     

土力学的现状及其数值分析方法中某些问题的讨论

概括性地论述了土力学目前的状况,指出土力学仍然处于半理论、半经验的发展水平。论述了土力学数值分析方法在工程应用中存在的问题,并对如何正确地使用土力学的数值分析方法和改变土力学数值分析在工程应用中的不利状况进行了讨论。另外,就如何改善这种状况提出了一些措施。从这些讨论可以得到以下结论,土力学还处于发展的初期阶段,随着土力学的理论及其数值分析方法的不断发展,以及其使用经验的不断积累,计算土力学在工程实践中将会得到越来越多的应用,它必将成为岩土工程分析中的有力工具。     

Ferric iron partitioning between plagioclase and silicate liquid: thermodynamics and petrological applications

A series of Fe and Mg partition experiments between plagioclase and silicate liquid were performed in the system SiO₂-Al₂O₃-Fe₂O₃-FeO-MgO-CaO-Na₂O under oxygen fugacities from below the IW buffer up to that of air. A thermodynamic model of plagioclase solid solution for the (CaAl,NaSi,KSi)(Fe³⁺,Al³⁺)Si₂O₈-Ca(Fe²⁺,Mg)Si₃O₈ system is proposed and is calibrated by regression analysis based on new and previously reported experimental data of Fe and Mg partitioning between plagioclase and silicate liquid, and reported thermodynamic properties of end members, ternary feldspar and silicate liquid. Using the derived thermodynamic model, FeO^t, MgO content and Mg/(Fe^t+Mg) in plagioclase can be predicted from liquid composition with standard deviations of ǂ.34 wt% (relative error =9%) and ǂ.08 wt% (14%) and ǂ.7 (8%) respectively. Calculated Fe³⁺-Al exchange chemical potentials of plagioclase, m_{Fe^{3 +} ( Al )- 1} ^Pl{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Pl}} agree with those calculated using reported thermodynamic models for multicomponent spinel, m_{Fe^{3 +} ( Al )- 1} ^Sp{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Sp}} and clinopyroxene, m_{Fe^{3 +} ( Al )- 1} ^Cpx{\rm \mu }_{{\rm Fe}^{{\rm 3 + }} \left( {{\rm Al}} \right)_{{\rm - 1}} }^{{\rm Cpx}} . The FeO^t content of plagioclase coexisting with spinel or clinopyroxene is affected by Fe³⁺/(Fe³⁺+Al) and Mg/(Fe+Mg) of spinel or clinopyroxene and temperature, while it is independent of the anorthite content of plagioclase. Three oxygen barometers based on the proposed model are investigated. Although the oxygen fugacities predicted by the plagioclase-liquid oxygen barometer are scattered, this study found that plagioclase-spinel-clinopyroxene-oxygen and plagioclase-olivine-oxygen equilibria can be used as practical oxygen barometers. As a petrological application, prediction of plagioclase composition and f_O2 are carried out for the Upper Zone of the Skaergaard intrusion. The estimated oxygen fugacities are well below QFM buffer and consistent with the estimation of oxidization states in previous studies.     

基于实时更新数值模型的高陡边坡动态仿真分析方法及应用

在水电工程中,河谷两遍的工程边坡具有尺度大、高而陡的特点。同时这些边坡服务年限长,稳定性要求高。对于耦合考虑施工进度和工程安全的动态分析还较少涉及。由于地质参数、边坡结构和开挖支护的实时变化,对边坡工程安全的实时评价已成为制约水利水电工程建设安全的关键技术问题。考虑边坡施工过程中开挖进度的调整、新地质信息的揭露、支护方案的修改及计算参数的修正等信息实时变化的客观实际,重点针对目前较少涉及的施工期高边坡工程安全的实时分析和控制问题,提出了基于数值模拟技术的边坡实时安全仿真方法;基于数据库技术和ABAQUS二次开发技术实现施工进度信息、地质信息及支护信息到数值计算模型的动态映射,开发了边坡安全实时分析系统,最终实现了与施工进程相适应的边坡安全状态实时仿真分析和预测。     

质谱计逐级蒸发-沉积测定单颗粒锆石年龄原理及可靠性的证明——兼论锆石同位素地质年代学的适用性

讨论了锆石由于铅丢失造成自晶体表层到晶体内部放射成因207Pb/206Pb的变化趋势，热离子质谱计直进行单颗锆石石蒸发-沉积测年时，锆石在灯丝上加热蒸发和在冷灯丝上沉积时物质的分布状况，锆石中的U-Pb体系是否为封闭体系的判别标志，并对锆石年代学的适用性进行了讨论。     

天目湖热力学状况的监测与分析

根据2001年6月～2002年5月在天目湖进行的每月一次的全湖水温观测以及2002年5月25～26日在1号点的定点24h连续水温观测资料,分析研究了天目湖水温时空分布、垂直分布、温跃层的形成与变化以及湖泊热力学对溶解氧、营养盐的影响,揭示了各自的基本变化规律。结果表明:天目湖水温的日变化、年变化及垂直分布明显,而其水平差异不大;温跃层随着季节的变化而呈现增强-稳定-减弱-消失的周期变化;水温的变化以及温跃层的存在引起了水中溶解氧以及营养盐的变化。     

多沙河流的非恒定一维水沙数学模型及其应用

介绍多沙河流上非恒定一维水沙数学模型采用网格上流量与水位交错布置的方案,应用Total Variation Diminishing(TVD)格式求解基本方程,模型具有守恒性好、无虚假振荡、算法简单的优点,克服了断面急剧变化时计算容易失稳的困难。同时,由于采用了适合于多沙河流的泥沙运动基本方程及相关参数,因此,该一维模型能够适应多沙河流高浓度输沙及大幅度冲淤变化的特征,成功模拟了黄河干、支流上的洪水传播。     

同位素双稀释剂技术的数值模拟与应用

近年来,非传统稳定同位素地球化学已成为新兴发展的研究领域,同位素双稀释剂技术是该领域的主流分析方法之一,但始终存在的困扰问题是如何进行同位素双稀释剂间以及稀释剂与样品间的最优混合。这两种混合比严重制约着同位素比值测试的准确度、精度和重现性。本文从同位素双稀释剂技术的数学原理出发,以Cr、Mo同位素为例,结合MonteCarlo误差模型,深入讨论和分析了两单稀释剂间比值(λ)及双稀释剂与自然样品间比值(p)的最优化数值模拟,得到了以双稀释剂~50(Cr)-~54(Cr)、~97(Mo)-~(100)Mo为例的λ、p值的最优取值范围。该数值模拟可适用于所有能够利用同位素双稀释剂技术的元素。在此基础上,就数值模拟结果如何应用于同位素双稀释剂的选择和组成标定也给予了详细说明,这为我国开展和使用同位素双稀释剂技术提供了基本思路和需要注意的关键问题,有助于推动同位素双稀释剂技术在我国地学界的广泛应用和发展。     

Fuzzy probabilistic seismic hazard analysis with applications to Kunming city,China

China is prone to highly frequent earthquakes due to specific geographical location, which could cause significant losses to society and economy. The task of seismic hazard analysis is to estimate the potential level of ground motion parameters that would be produced by future earthquakes. In this paper, a novel method based on fuzzy logic techniques and probabilistic approach is proposed for seismic hazard analysis (FPSHA). In FPSHA, we employ fuzzy sets for quantification of earthquake magnitude and source-to-site distance, and fuzzy inference rules for ground motion attenuation relationships. The membership functions for earthquake magnitude and source-to-site distance are provided based on expert judgments, and the construction of fuzzy rules for peak ground acceleration relationships is also based on expert judgment. This methodology enables to include aleatory and epistemic uncertainty in the process of seismic hazard analysis. The advantage of the proposed method is in its efficiency, reliability, practicability, and precision. A case study is investigated for seismic hazard analysis of Kunming city in Yunnan Province, People’s Republic of China. The results of the proposed fuzzy logic-based model are compared to other models, which confirms the accuracy in predicting the probability of exceeding a certain level of the peak ground acceleration. Further, the results can provide a sound basis for decision making of disaster reduction and prevention in Yunnan province.     

几种梯度场计算的数值算法及其地学应用

梯度作为标量场在空间上最重要的微观变化特征,能较好地反映其在空间上的变化趋势。本文介绍了中心差分算子、相邻梯度因子和Sobel梯度算子三种梯度数值算法。以点电荷在三维空间中产生的电势场构造了一标量场,通过理论公式求解了其梯度值,并将理论梯度值与三种数值算法的结果进行了对比与统计分析,验证了三种梯度数值算法的有效性和适用性。结果表明,排除场源的影响后,Sobel梯度算子相较其它两种数值算法更逼近于理论梯度值。最后,利用Sobel梯度算子求取了某铜矿区可控源音频大地电磁法(CSAMT)三维视电阻率场的梯度场,用梯度场刻画了视电阻率在空间上的变化趋势,提高了对三维地下空间的电性分布规律和地质情况的认识。     

Large deformation dynamic analysis of saturated porous media with applications to penetration problems

This paper outlines the development as well as implementation of a numerical procedure for coupled finite element analysis of dynamic problems in geomechanics, particularly those involving large deformations and soil-structure interaction. The procedure is based on Biot’s theory for the dynamic behaviour of saturated porous media. The nonlinear behaviour of the solid phase of the soil is represented by either the Mohr Coulomb or Modified Cam Clay material model. The interface between soil and structure is modelled by the so-called node-to-segment contact method. The contact algorithm uses a penalty approach to enforce constraints and to prevent rigid body interpenetration. Moreover, the contact algorithm utilises a smooth discretisation of the contact surfaces to decrease numerical oscillations. An Arbitrary Lagrangian–Eulerian (ALE) scheme preserves the quality and topology of the finite element mesh throughout the numerical simulation. The generalised-α method is used to integrate the governing equations of motion in the time domain. Some aspects of the numerical procedure are validated by solving two benchmark problems. Subsequently, dynamic soil behaviour including the development of excess pore-water pressure due to the fast installation of a single pile and the penetration of a free falling torpedo anchor are studied. The numerical results indicate the robustness and applicability of the proposed method. Typical distributions of the predicted excess pore-water pressures generated due to the dynamic penetration of an object into a saturated soil are presented, revealing higher magnitudes of pore pressure at the face of the penetrometer and lower values along the shaft. A smooth discretisation of the contact interface between soil and structure is found to be a crucial factor to avoid severe oscillations in the predicted dynamic response of the soil.     

Three-dimensional numerical analysis of deep excavations with cross walls

Previous plane strain analysis of a case history has shown that cross walls in an excavation can effectively reduce movements induced by deep excavation. This study performed three-dimensional numerical analyses for 4 deep excavation cases with different installations of cross walls, including different excavation depths, cross wall intervals and cross wall depths. Both the observed and computed wall deflections for the 4 cases were compared with those of the same excavations that were assumed with no cross walls installed to demonstrate the effectiveness of cross walls in reducing lateral wall deflections. The results show that the cross wall also had a corner effect similar to that of the diaphragm wall. The deflection of the diaphragm wall was smallest at the location of the cross wall installed and then increased with the increasing distance from the cross wall, up to the midpoint between two cross walls. Many factors such as in situ soil properties, diaphragm wall properties, construction procedure, cross wall depth and so on may affect the amount of reduction in lateral wall deflections due to the installation of cross walls. Under the same condition, the amount of reduction was highly dependent on the depth of cross walls, distance to the cross walls and the cross wall interval.