2003年北黄海5.1级地震的地震活动特征

对2003年北黄海5．1级地震的地震活动背景、地震序列、震中周围中小地震活动图像及地震学参数等进行了详细分析。认为该次地震为前-主-余型，经历了一个活跃-平静-发震的过程；地震前震中附近出现了孕震空区和一些地震活动性参数异常。     

Simple formulas for the dynamic stiffness of pile groups

Simple formulas are derived for the dynamic stiffness of pile group foundations subjected to horizontal and rocking dynamic loads. The formulations are based on the construction of a general model of impedance matrices as the condensation of matrices of mass, damping, and stiffness, and on the identification of the values of these matrices on an extensive database of numerical experiments computed using coupled finite element–boundary element models. The formulations obtained can be readily used for the design of both floating piles on homogeneous half‐space and end‐bearing piles and are applicable for a wide range of mechanical and geometrical parameters of the soil and piles, in particular for large pile groups. For the seismic design of a building, the use of the simple formulas rather than a full computational model is shown to induce little error on the evaluation of the response spectra and time histories. Copyright © 2009 John Wiley & Sons, Ltd.     

Electrical aerosol spectrometer of Tartu University

The electrical aerosol spectrometer (EAS) of the parallel measuring principle at Tartu University is an efficient instrument for rapid measurement of the unstable size spectrum of aerosol particles. The measuring range from 10 nm to 10 μm is achieved by simultaneously using a pair of differential mobility analyzers with two different particle chargers. The particle spectrum is calculated and measurement errors are estimated in real time by using a least-squares method. Experimental calibration ensures reliability of measurement. The instrument is well suited for continuous monitoring of atmospheric aerosol.     

井中激电地-井方式井旁球体正反演

以激发极化法中的地-井方式中的体极化球形体作为研究对象,采用解析法对地下半空间的电位场进行求解,利用提取出的井中二次场异常电位差进行反演,得出极化球体的球心离井的距离和球心的埋深。通过计算机进行正反演模拟,验证了进行正反演计算的解析式的正确性,并从物理意义上对正反演过程中出现的现象进行分析。最后给出该方法的适用条件以及注意事项。研究对激发极化法地下勘探目标的定量分析具有明确意义。     

2011年江苏南部两场大暴雨对比分析

利用高分辨率的加密气象自动站资料、FY2D卫星资料、多普勒雷达资料、常规观测资料以及6 h 1次的NCEP再分析资料等,对2011年6月18日和2011年7月18日江苏地区分别发生在梅雨期开始阶段和结束阶段的两场暴雨进行中尺度天气系统演变和雷达回波参数等特征的对比分析。结果表明:(1)6月18日的天气形势是典型的梅雨期降水形势,在梅雨锋附近产生了区域性暴雨。水汽输送主要是对流层中低层的西南暖湿气流。7月18日的局地暴雨则是出现在低压倒槽顶端右侧的偏东气流中。(2)两次暴雨过程强降水发生前都存在对流层低层辐合快速增强的过程。7月18日暴雨强降水发生前散度值下降则更为迅速。(3)两次暴雨过程中强降水区都出现在地面辐合系统附近的东北气流中,且随着地面辐合系统移动。(4)两次暴雨过程都出现了T_BB低于-62℃的强对流云团。(5)6月18日,与多个线性排列的"逆风区"对应的强回波中心形成了"列车效应";7月18日,对流回波带上单体不断流入,在低空急流左前端合并成团状强对流区,分别是形成两次暴雨的重要原因。     

Simultaneous identification of parameter,initial condition,and boundary condition in groundwater modelling

In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

论多种数据联合反演的模式及算法

针对日益丰富的数据资料，综述了大地测量反演的发展进程，论述了多种数据联合反演的意义和必要性，给出了兼顾权比确定的大地测量、地震和地质三类数据的联合反演模式及采用水准测量、GPS和重力三种数据联合反演断层参数的具体公式，特别介绍了求解优化问题所有全局最优解的区间算法，并给出了基本算例，结果体现了区间算法无可比拟的优越性及其在大地测量反演领域中的应用前景。     

Intel MKL在Bernese GNSS数据处理中的应用

为提高Bernese GNSS software数据处理效率，将英特尔数学核心函数库(math kernel library，MKL)应用于Bernese 精密定轨数据处理，对比分析多个MKL矩阵求逆函数与Bernese 原有程序的计算效率。使用2019-03全球200个测站北斗/GNSS数据进行实验分析，结果表明，采用参数预消除策略时，参数预消除步骤消耗时间明显大于矩阵求逆，使用MKL处理数据效率提升不明显；而未采用参数预消除策略时，使用MKL矩阵求逆函数可显著提高矩阵求逆效率，其中dpotri函数矩阵求逆计算效率最高，消耗时间平均值为133 s，相比Bernese 原有程序计算速度可提高13倍。     

多面函数参数自适应选取方法在GPS高程拟合中的应用

提出一种综合的多面函数参数自适应选取方法。该方法利用核函数阶数和平滑因子存在最优值以及核函数结点分布应符合高程异常变化趋势的特点,综合利用正交化算法选取核函数结点以及二维粒子群算法确定最佳阶数和平滑因子,实现多面函数参数的完全自适应选取。将该方法应用于地形起伏差异不同的2个测区,结果表明,相比传统经验性方法和部分参数自适应方法,该方法用于GPS高程拟合的精度和可靠性更高。     

青藏高原大尺度地表流体的重力效应特征分析

基于2002~2017年三维大气再分析模型、水文模式及GRACE卫星重力数据,采用负荷格林函数和卫星重力反演方法分析青藏高原大尺度地表流体的重力效应特征。结果表明,青藏高原地区大气扣除季节项前后的重力效应峰对峰变化介于1.9~10.7 μGal和1.1~4.6 μGal,其负荷效应呈现出显著的季节性特征;土壤水和积雪扣除季节项前后的重力效应峰对峰变化介于0.9~9.9 μGal和0.7~8.3 μGal,其负荷效应呈现出显著的季节性和年际变化特征;GRACE时变重力场扣除季节项前后的重力变化峰对峰变化介于4.4~24.1 μGal和3.4~18.3 μGal,考虑观测误差、构造运动和信号泄漏误差的影响,其显著的季节性和年际变化反映了多种地表流体的综合影响。