高速铁路新型路堑基床结构动力响应分析

研制出一种新的防排水结构层用以控制膨胀土路基的含水率变化,从而达到对新建云—桂高速铁路膨胀土地段整治的目的。本文采用FLAC~(3D)有限差分软件研究该新型路堑基床结构(基床中设置新型防排水结构层)在列车荷载作用下动力响应规律,并利用现场试验结果对数值研究结果进行验证。研究表明:新型基床中道砟和新型防排水结构层对动应力衰减贡献较大;新型防排水结构层对竖向动应力、剪应力衰减作用明显,增强了基床结构的动力稳定性;新型基床结构能够更有效地控制基床的动位移;新型基床中振动速度和加速度随深度的增加而减小,道砟对振动速度和加速度的衰减作用明显。新型路堑基床结构振动速度和动应力的现场实测与数值分析规律基本一致,且实测值与计算值大小相近。研究成果可为特殊土地区高速铁路基床的设计、施工及其动力响应研究提供参考。     

Geochemical features of the silt sediment of the Novosibirsk Reservoir

Earlier data on heavy metal concentrations in bottom sediments of the Novosibirsk Reservoir are supplemented by new geochemical characteristics of fine silt sediment in the near-dam part of the reservoir with the list of chemical elements extended to 55. After the beginning of the Ob’s runoff regulation by HPP dam, a new geochemical type of sediment started forming in the new reservoir. This type of sediment features higher concentrations of chemical elements (including heavy metals), reaching the values typical of Baikal silt recorded in the central hollow of the lake.     

中国大陆东南沿海地震带地震迁移规律

用一种新思路、新方法研究了中国大陆东南沿海地震带地震的迁移，结果表明其具有由东向西定向迁移的规律，其各个迁移过程的迁移速度也相对稳定，平均值为４２ｋｍ／ａ，这为该带地震活动的动力源和目前可能的地震危险区提供了一定的依据。     

Influence of the topography of the surface of the Earth on vertical sounding of the temperature profile

Atmospheric temperature and humidity fields as well as information on other meteorological parameters are nowadays retrieved from radiance measurements recorded by operational meteorological satellites. Up to now, the inversion procedures used only take into account crude information on the topography of the Earth’s surface. However, the applied radiative transfer codes have to consider the Earth’s surface as the lower boundary of the atmospheric model and, therefore, need a more precise mean elevation and a classification of the roughness of the Earth’s surface. The influence of the topography of the Earth surface on retrieved temperature profiles is studied by using a physico-statistical inversion method. An objective analysis is made of the more precise mean elevation and derivation of roughness parameters using a new high-resolution digital elevation model (DEM) with a resolution of 500 m×500 m. By means of a geomorphological process and a newly developed topography rejection test, areas with a high surface roughness are localized and singled out. The influence of topography on the retrieved temperature profiles is illustrated by case studies. Changes are found predominantly in areas with a high variation of topography. Using the new high-resolution DEM and the topography rejection test, the geographical position of the calculated temperature profiles tends to be shifted towards areas with a small vertical variation of topography. The mean elevation determined by the new elevation model better characterizes the area observed. Hence, the temperature profiles can be calculated down to lower atmospheric levels. Furthermore, a guess profile better describing the atmospheric situation is selected by the more precise elevation. In addition, the temperature profiles obtained near the coast are improved considerably by the more precise determination of the surface property ‘sea’ and ‘land,’ respectively. Integration of an independent physical information such as topography leads, on average, to a slight improvement of the results of the physico-statistical inversion procedure. In some cases, however, significant improvements have been achieved regarding the desired accuracy of temperature profiles of the order of 1 K. In future, the spatial resolution of new high-resolution sounding instrumentation on the next generation of operational meteorological satellites will be increased. To exploit the resolving power of this new instrumentation, the different variation of the topography of the Earth surface, especially in regions with a high variation of topography, can be taken into account more precisely by using a high-resolution DEM.     

The Rhythm of Earthquake Activity and an Analysis on the Earthquake Tendency of the World, the Chinese Continent and North China

Based on the Morlet complex wavelet transformation, the authors put forward a kind of new method for distinguishing periods of seismic activity and quietude and a new physical thought on the time-dependent wavelet accumulation energy spectrum with periods, the time-frequency distribution of wavelet vibration period spectrum and period-specific wavelet vibration spectrum. By applying the above methods to a time series which is composed of earthquake accumulation energy per year for the world, the Chinese continent and North China,respectively, we obtained some new information about the rhythm of shallow earthquake activity. Considering the historic earthquakes and the rhythm characteristics of current strong earthquake activity, the earthquake tendency in the next years is discussed.     

新国家标准震级标度与传统震级标度对比研究

使用四川省地震台网2013年4月以来的宽频带数字地震资料,对新国家标准震级标度与传统震级标度中的4类震级进行了对比分析。结果表明:新M_L与传统M_L、M_S（BB）与M_S7之间的差别最小,相关性最高;所有事件的新m_b均大于或等于传统m_b;大部分事件的新M_L大于传统M_L;m_B（BB）大于m_B的事件多于m_B（BB）小于m_B的事件;M_S（BB）小于M_S7的事件数量与M_S（BB）大于M_S7的事件接近;新M_L较传统M_L更符合区域地质实际状况,新m_b较传统m_b更便于测量,M_S（BB）较M_S7、m_B（BB）较m_B更适应当前的宽频带数字地震资料,其测定结果更稳定。新国家标准震级标度充分考虑了我国震级测定的历史连续性,很好地继承和衔接了传统震级标度,在可操作性和资料适用性方面有极大提高。      

设计速度传感反馈地震计的新方案新思路

在深入分析研究了速度传感反馈地震计传统设计方案的基础上，提出了技术性能更优的传递函数两步实现方案。两步方案为反馈环路设计带来的巨大灵活性，又促使我们跳出传统的“延长周期”模式的思路束缚，提出了单纯阻尼反馈环的设计新思路。理论计算表明，两步方案与阻尼环结合的设计可以使速度传感反馈地震计的技术性能达到更高境界。反馈深度大、环路稳定性好，动态范围大和极点稳性高使它明显优于传统设计方案。     

基于二次场的可控源电磁法三维有限元-无限元数值模拟

本文提出了一种新的混合有限元-无限元三维可控源电磁法(CSEM)问题快速高精度正演模拟算法.首先从电场双旋度方程出发,推导了水平电偶极子源的二次场边值问题,采用无限元代替截断边界条件和有限元离散内部计算区域的新策略,达到减小计算区域的目的,基于并行直接求解技术,实现多源CSEM问题的快速精确求解.其次,通过层状解析模型测试,一方面验证了新算法的正确性,另一方面通过与其他三种已知CSEM问题求解策略进行对比,表明了本文提出的基于二次场有限元-无限元算法具有离散区域小、求解速度快和计算精度高等优点.最后,通过3D模型计算,清晰直观地模拟了场源阴影效应,为野外数据的处理与解释提供指导.     

Analytical solutions for the coefficient of variation of the volume-averaged solute concentration in heterogeneous aquifers

Under the assumption that local solute dispersion is negligible, a new general formula (in the form of a convolution integral) is found for the arbitrary k-point ensemble moment of the local concentration of a solute convected in arbitrary m spatial dimensions with general sure initial conditions. From this general formula new closed-form solutions in m=2 spatial dimensions are derived for 2-point ensemble moments of the local solute concentration for the impulse (Dirac delta) and Gaussian initial conditions. When integrated over an averaging window, these solutions lead to new closed-form expressions for the first two ensemble moments of thevolume-averaged solute concentration and to the corresponding concentration coefficients of variation (CV). Also, for the impulse (Dirac delta) solute concentration initial condition, the second ensemble moment of thesolute point concentration in two spatial dimensions and the corresponding CV are demonstrated to be unbound. For impulse initial conditions the CVs for volume-averaged concentrations axe compared with each other for a tracer from the Borden aquifer experiment. The point-concentration CV is unacceptably large in the whole domain, implying that the ensemble mean concentration is inappropriate for predicting the actual concentration values. The volume-averaged concentration CV decreases significantly with an increasing averaging volume. Since local dispersion is neglected, the new solutions should be interpreted as upper limits for the yet to be derived solutions that account for local dispersion; and so should the presented CVs for Borden tracers. The new analytical solutions may be used to test the accuracy of Monte Carlo simulations or other numerical algorithms that deal with the stochastic solute transport. They may also be used to determine the size of the averaging volume needed to make a quasi-sure statement about the solute mass contained in it.     

Dispersion-relationship-preserving seismic modelling using the cross-rhombus stencil with the finite-difference coefficients solved by an over-determined linear system

Finite-difference modeling with a cross-rhombus stencil with high-order accuracy in both spatial and temporal derivatives is a potential method for efficient seismic simulation. The finite-difference coefficients determined by Taylor-series expansion usually preserve the dispersion property in a limited wavenumber range and fixed angles of propagation. To construct the dispersion-relationship-preserving scheme for satisfying high-wavenumber components and multiple angles, we expand the dispersion relation of the cross-rhombus stencil to an over-determined system and apply a regularization method to obtain the stable least-squares solution of the finite-difference coefficients. The new dispersion-relationship-preserving based scheme not only satisfies several designated wavenumbers but also has high-order accuracy in temporal discretization. The numerical analysis demonstrates that the new scheme possesses a better dispersion characteristic and more relaxed stability conditions compared with the Taylor-series expansion based methods. Seismic wave simulations for the homogeneous model and the Sigsbee model demonstrate that the new scheme yields small dispersion error and improves the accuracy of the forward modelling.     

Volcanic rocks and the depth of the Benioff zone in Sud Lipez,Southwest Bolivia

The depth of the Benioff zone has been determined from the new analysis of volcanic rocks and from new seismic data from S.W. Bolivia. Both data give consistently 260 km as a depth of the Benioff zone and for the origin of magma of calc-alkali magma series in SW-Bolivia. These inferences are consistent with the plate tectonics model.     

超高密度电阻率法在土洞、煤窑采空区和岩溶勘探中应用研究

超高密度电阻率法是Australia ZZ Resistivity Imaging研发中心积近十年的研究成果提出的全新的电阻率法勘探理念,是电法勘探从手工时代向计算机时代的一次跨越式发展,采用一次布极,采集所有可能组合电极间的电位信息,信息量较常规电法在信息量上翻了40倍以上,并且避开视电阻率概念,利用处理软件直接反演获取断面真电阻率值,针对目前还没有行之有效的地面物探方法可以勘探清楚的土洞、煤窑采空区和岩溶等复杂地质问题,我们采用该方法分别对机场岩溶和土洞、煤矿运煤巷道、隐伏小煤矿采空区和复杂地质病害隧道进行了试验勘探,经钻探验证,准确率在60%以上,达到国内领先水平,取得了较好效果,对复杂地质的工程勘探具有一定参考价值.     

复杂山地随机介质GMM-ULTI法射线追踪

对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律.     

顾及电离层变化的层析反演新算法

区别于以往GPS电离层层析研究主要关注迭代模型的思路,本文从两方面入手提高GPS电离层层析迭代算法的反演精度:一方面,顾及传统电离层层析迭代模型仅与对电子密度误差起放大作用的GPS射线截距权重相关的不足,提出考虑层析像素格网中的电子密度对GPS TEC的贡献建立新的迭代模型,在不同电子密度像素格网内重新分配GPS TEC实测值与其反演值之间的差距;另一方面,顾及电离层层析迭代算法中松弛因子对反演结果的影响,提出考虑电子密度变化构造新的松弛因子,抑制传播噪声对电子密度反演精度的影响.实验结果显示,相对于传统代数重构算法(ART),新方法反演的电离层电子密度剖面更接近于电离层测高仪观测的电子密度剖面,提高了电子密度反演精度.     

A length-biased version of the Birnbaum-Saunders distribution with application in water quality

In this article, we develop a new model based on the Birnbaum-Saunders distribution that results to be both useful and practical for environmental sciences. The density, distribution and hazard functions, moments and properties of this new model are presented. A graphical analysis of the density is also provided. Furthermore, we estimate parameters, propose asymptotic inference and discuss influence diagnostics by using likelihood methods for the new distribution. An illustrative example with real data related to water quality indicates the adequacy on the new distribution.     

新一代系统CDSN应用地震数据的格式及其使用方法

ＣＤＳＮ正处在更新换代的重要时期，为了使广大地震研究人员更好地应用ＣＤＳＮ资料，本文详细叙述了更新的ＣＤＳＮ地震事件数据格式，提出了新数据与已发行数据格式相兼容的办法，同时说明了有关格式转换程序的使用方法。     

一种计算固体潮观测数据改正值的新方法

从标定前一天和标定当天的地形变前兆数据变化趋势入手，在分析当天数据结构的基础上采用最小二乘法原理，提出了一种计算数采改正值的新方法。经典型实例计算表明，该方法精度较高，效果明显。并对误差的来源作了进一步的分析。     

为新区划图编制所建立的地震动衰减关系

介绍了建立新一代地震区划图所采用的地震动参数衰减关系的总体思路,并从资料、衰减关系分区、衰减关系模型、回归方法、转换等方面说明了地震动参数衰减关系的建立过程,给出了我国分区地震烈度和地震动参数衰减关系结果。新的地震动衰减关系的建立,具有如下特点:一是基于更加丰富可靠的强震记录和烈度资料;二是采用了具有大震近场饱和特征的地震动衰减模型;三是采用了使结果更加稳定的分步回归方法;四是在地震动衰减关系分区时考虑了地震活动性特征。与第四代地震区划图衰减关系相比,由于地震动衰减模型的变化和高震级强震记录的增加,高震级下的峰值加速度有所降低,而中强地震区的峰值加速度则在低震级时有所提高。     

Features of the solar-type activity of stars with various depths of the convective zone

The levels of chromospheric and coronal activity of the Sun are compared with new vast observations of late-type stars. The solar chromosphere turned to be more powerful than in the main body of stars and the corona is considerably weakened. A wavelet analysis of activity indices and measurement results for the magnetic field of the Sun as a star was performed for several solar cycles. It was obtained that solar activity in the differential rotation differs from phenomena on less massive K stars with cycles and, in contrast to them, the large-scale magnetic field of the Sun is a regulating factor for active processes. These results can be naturally explained with the assumption that the activity of a star with a given mass depends on the depth of the lower base of the convective zone. This seems to require the development of knowledge about the two-level dynamo and a new approach to studying solar-type activity.     

结构地震响应分析新方法

本文通过对传统动力反应分析方法的总结,阐明了建立隐式和显式方法的一般思路及数学本质,提出了使用系统位移反应向量三阶导数的实用显式高阶积分方法,分析了该显式方法的精度和稳定性,并对建立更高阶隐式和显式方法以及方法的精度和稳定性作了初步讨论。结果表明,本文方法具有明显的优点。