细胞自动机地震波模拟的并行化算法

采用区域分解的办法对格子气自动机地震波传播模拟的并行化算法进行了研究,提出了两种区域分解算法,并在SGIPowerChallenge并行计算机上进行了实际地质模型中平面波反射地震响应的计算,获得了较高的并行计算效率．     

External atmospheric corrections in geodetic very-long-baseline interferometry

Three methods to correct for the atmospheric propagation delay in very-long-baseline interferometry (VLBI) measurements were investigated. In the analysis, the NASA R&D experiments from January 1993 to June 1995 were used. The methods were compared in correcting for the excess propagation delay due to water vapour, the “wet” delay, at one of the sites, the Onsala Space Observatory on the west coast of Sweden. The three methods were: (1) estimating the wet delay using the VLBI data themselves; (2) inferring the wet delay from water vapour radiometer (WVR) data, and (3) using independent estimates based on data from the global positioning system (GPS). Optimum elevation cutoff angles were 22^∘ and 26^∘ when using WVR and GPS data, respectively. The results were found to be similar in terms of reproducibility of the estimated baseline lengths. The shortest baselines tend to benefit from external measurements, whereas the lack of improvement in the longer baselines may be partly due to the large amount of data thrown away when removing observations at low elevation angles. Over a 2 week period of intensive measurements, the two methods using external data showed an overall improvement, for all baseline lengths, compared to the first method. This indicates that there are long-term systematic errors in the wet delay data estimated using WVR and GPS data. Received: 27 October 1998 / Accepted: 20 May 1999     

The application of a shift theorem analysis technique to multipoint measurements

A Fourier domain technique has been proposed previously which, in principle, quantifies the extent to which multipoint in-situ measurements can identify whether or not an observed structure is time stationary in its rest frame. Once a structure, sampled for example by four spacecraft, is shown to be quasi-stationary in its rest frame, the structure’s velocity vector can be determined with respect to the sampling spacecraft. We investigate the properties of this technique, which we will refer to as a stationarity test, by applying it to two point measurements of a simulated boundary layer. The boundary layer was evolved using a PIC (particle in cell) electromagnetic code. Initial and boundary conditions were chosen such, that two cases could be considered, i.e. a spacecraft pair moving through (1) a time stationary boundary structure and (2) a boundary structure which is evolving (expanding) in time. The code also introduces noise in the simulated data time series which is uncorrelated between the two spacecraft. We demonstrate that, provided that the time series is Hanning windowed, the test is effective in determining the relative velocity between the boundary layer and spacecraft and in determining the range of frequencies over which the data can be treated as time stationary or time evolving. This work presents a first step towards understanding the effectiveness of this technique, as required in order for it to be applied to multispacecraft data.     

Method discussion for quick response grey prediction of stronger aftershocks of an earthquake sequence

In this paper, we take occurrence process of early strong aftershocks of a main-after shock type′s earthquake sequence as a complex grey system, and introduce predicting method for its stronger aftershocks by grey predicting theory. Through inspection prediction for 1998 Zhangbei MS=6.2 earthquake sequence, it shows that the grey predicting method maybe has active significance for the investigation of quick response prediction problems of stronger aftershocks of an earthquake sequence.     

A Probabilistic Model for Fatigue Crack Propagation Analysis

A simple probabilistic model for predicting crack growth behavior under random loading is pre-sented.In the model,the parameters c and m in the Paris-Erdogan Equation are taken as random variables,and their stochastic characteristic values are obtained through fatigue crack propagation tests on an offshorestructural steel under constant amplitude loading.Furthermore,by using the Monte Carlo simulation tech-nique,the fatigue crack propagation life to reach a given crack length is predicted.The tests are conducted toverify the applicability of the theoretical prediction of the fatigue crack propagation.     

南海西北部水域斑节对虾资源的调查研究

1989－1991年对南海西北部水域进行斑节对虾Penaeusmonodon资源专项调查。结果表明,斑节对虾的密集分布区在海南岛东部和南部海区,其现存资源量为82．6万尾,周年有体长30mm以下的幼虾出现。主要产卵期在9－12月。产卵场在海南岛东部和南部水深30-60m的海区。     

完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

柴达木盆地东部地震地面运动放大效应

柴达木盆地是青藏高原东北部大型断陷山间盆地,该地区的流动观测记录了2008年11月10日发生于大柴旦附近的M_W6.3地震。和附近的基岩上的记录相比,盆地内部的记录显示出非常显著的地面运动放大效应,表现为峰值速度的增大、持续时间的延长,其呈现出长持续时间的后续震相。傅里叶频谱分析表明盆地内部显著的后续震相的频率和直达波相比较低,地面质点运动轨迹图显示后续震相为面波运动特征。为了解释地面运动的差异,构建二维模型,通过交错网格高阶有限差分方法计算了地震波在盆地内部的传播过程,结果显示盆地内部低速层的存在造成直达波的放大以及多次反射与转换,盆地边缘结构造成的波的相干叠加产生了强烈的次生面波,其低频、大振幅、长持续时间的特征是盆地内部地面运动放大的主要原因。     

云南禄劝地震震源处的应力场求解

采用前人反演得到的云南禄劝地震的地震矩张量和主震及部分余震的震源机制解,以该主震震中为中心在全球CMT目录中查询到的部分地震的震源机制解,先将地震矩张量转化为震源机制解,运用精细网格搜索反演方法将震源机制解反演得到禄劝地震和其周边地区的应力场。对反演得出的两个应力场进行差异性对比研究,结果表明禄劝地震震源处主压应力场为NNW—SSE向,主压应力与主张应力相当,周边地区的主压应力方向为NW—SE,断层破裂面倾角大,以走滑正断层为主,主压应力占优势。但是由于云南地区主要由NNW—SSE和NW—SE方向的主压应力控制,并受本文所选的周边地区的经纬度及所处地区的控制,所以该区域在总体上受NW—SE方向的压应力控制,局部地区受NNW—SSE方向的压应力控制。该结果可以用来分析该地区的地震地质背景和断层形成条件,对地球动力环境的研究有一定的意义。     

Experimental study of aeolian sand ripples in a wind tunnel

The topographic parameters and propagation velocity of aeolian sand ripples reflect complex erosion, transport, and deposition processes of sand on the land surface. In this study, three Nikon cameras located in the windward (0–1 m), middle (4.5–5.5 m), and downwind (9–10 m) zones of a 10 m long sand bed are used to continuously record changes in sand ripples. Based on the data extracted from these images, this study reaches the following conclusions. (1) The initial formation and full development times of sand ripples over a flatbed decrease with wind velocity. (2) The wavelengths of full development sand ripples are approximately twice the wavelengths of initially formed sand ripples. Both wavelengths increase linearly with friction velocity. During the developing stage of sand ripples, the wavelength increases linearly with time. (3) The propagation velocity of full development sand ripples is approximately 0.6 times that of the initially formed sand ripples. The propagation velocity of both initial and full development of sand ripples increase as power functions with respect to friction velocity. During the developing stage of sand ripples, the propagation velocity decreases with time following a power law. These results provide new information for understanding the formation and evolution of aeolian sand ripples and help improve numerical simulations. Copyright © 2017 John Wiley & Sons, Ltd.