Singularity-free expression of magnetic field of cuboid under undulating terrain

Most of the current computing methods used to determine the magnetic field of a uniformly magnetized cuboid assume that the observation point is located in the upper half space without a source. However, such methods may generate analytical singularities for conditions of undulating terrain. Based on basic geomagnetic field theories, in this study an improved magnetic field expression is derived using an integration method of variable substitution, and all singularity problems for the entire space without a source are discussed and solved. This integration process is simpler than that of previous methods, and final integral results with a more uniform form. ΔT at all points in the source-free space can be calculated without requiring coordinate transformation; thus forward modeling is also simplified. Corresponding model tests indicate that the new magnetic field expression is more correct because there is no analytical singularity and can be used with undulating terrain.     

Newmark design spectra considering earthquake magnitudes and site categories

Newmark design spectra have been implemented in many building codes, especially in building codes for critical structures. Previous studies show that Newmark design spectra exhibit lower amplitudes at high frequencies and larger amplitudes at low frequencies in comparison with spectra developed by statistical methods. To resolve this problem, this study considers three suites of ground motions recorded at three types of sites. Using these ground motions, influences of the shear-wave velocity, earthquake magnitudes, source-to-site distances on the ratios of ground motion parameters are studied, and spectrum amplification factors are statistically calculated. Spectral bounds for combinations of three site categories and two cases of earthquake magnitudes are estimated. Site design spectrum coefficients for the three site categories considering earthquake magnitudes are established. The problems of Newmark design spectra could be resolved by using the site design spectrum coefficients to modify the spectral values of Newmark design spectra in the acceleration sensitive, velocity sensitive, and displacement sensitive regions.     

Identification of modal parameters of a 600-m-high skyscraper from field vibration tests

Accurate prediction of the dynamic responses of a high-rise building subjected to dynamic loads such as earthquake and wind excitations requires the information of its structural dynamic properties such as modal parameters including natural frequencies and damping ratios. This paper presents the identification results of the modal parameters based on field vibration tests on a 600-m high skyscraper. A set of tests, including ambient vibration test (AVT) and free vibration test (FVT), were conducted on the skyscraper to identify its modal parameters. Firstly, this paper presents and discusses the modal parameters of the skyscraper assessed by several identification methods applied to the AVT measurements. These methods include the wavelet transform (WT) method, the stochastic subspace identification (SSI) method, and the random decrement technique (RDT). Secondly, an active mass damper (AMD) system with total mass 1000 tons equipped into the skyscraper was used to excite the building for estimation of the modal parameters by FVT. Thirdly, this paper presents observations on the structural dynamic behavior of the skyscraper with the operation of the AMD system during a typhoon event. The field measurement results show that the AMD system functioned efficiently for suppression of the wind-induced vibrations of the skyscraper during the typhoon. This paper aims to further understand the structural dynamic properties of super-tall buildings and provide useful information for structural design and vibration control of future skyscrapers.     

Seismic amplitude inversion for the transversely isotropic media with vertical axis of symmetry

Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data.     

小台站的大数据和云计算设想

地球科学是一门观测科学,相关科学研究及公益服务,如气象、地震等,都依赖于野外台站的观测。目前,由于多种因素的作用,导致台站观测和后续分析应用严重脱节。本文基于大数据和云计算的思想,将台站工作和目前后续研究中心的工作进行有机结合,提出了一套新的工作思路,即将数据计算分析和区域监测预报的工作由台站来完成,然后由一个综合的数据中心负责存储和管理台站提供的"大数据"信息,并根据台站需要和反馈情况研发相关软件,为台站的分析计算提供云计算的平台。     

陇东黄土高原下垫面不均匀性指标的建立及其对大孔径闪烁仪(LAS)观测感热通量的影响

目前定量研究下垫面不均匀性对大孔径闪烁仪(LAS)、涡动相关仪(ECS)观测感热通量差异的影响还比较少.本文利用黄土高原庆阳观测站2012年6、7月典型晴天两主风向范围E-SE和SW-W的陆面过程数据,在建立了下垫面不均匀性定量化指标的基础上,分析了下垫面不均匀性与地表温度变率的关系及其对LAS、ECS观测感热通量差异的影响.结果表明:庆阳站下垫面不均匀性大小η和地表温度变率г具有很好的相关性,相关系数达到0.566以上,证明了新建立指标的合理性.下垫面不均匀性大小η和LAS、ECS观测的感热通量差异HLAS-HECS具有很好的一致性,相关系数达到0.634.下垫面不均匀性对LAS和ECS的感热通量差异影响显著,下垫面越不均匀LAS和ECS测量的感热通量差异越大.对两主风向分别进行分析,在E-SE风向范围η和HLAS-HECS的相关系数为0.430,HLAS和HEC拟合的线性趋势系数为1.279,在SW-W风向范围η和HLAS-HECS的相关系数为0.680,HLAS和HEC拟合的线性趋势系数为1.297.下垫面不均匀性的影响程度越大,LAS和ECS观测的感热通量差异越大.     

南海海底地震仪异常数据的分析和处理

海底地震仪(Ocean Bottom Seismometer, OBS)数据处理至关重要,是获取深部地壳结构的基础与前提.2006年实施OBS2006-2测线时,有2台OBS(OBS03,OBS06)数据出现异常,无法使用.由于海上航次花费巨大,采集到的数据弥足珍贵.本文采用数据格式检查、邻近台站对比分析、重采样等方法,成功地对这2台OBS数据进行了解编处理,得到了这两个台站的综合地震记录剖面;利用上述方法对2011年实施的OBS973-3测线中的异常台站OBS03进行了分析处理,同样得到了OBS03台站的综合地震剖面;通过查看两次海上实验班报发现,OBS2006-2测线之OBS06与OBS973-3测线之OBS03内部Sedis编号相同,为同一台记录仪器,再一次验证上述处理方法正确可行;然后对OBS2006-2测线2个台站进行震相识别与走时拾取后,利用前人纵波速度模型开展了射线追踪与走时模拟.此次对异常OBS数据的重新处理工作,不仅为OBS探测提供了宝贵的数据处理经验,而且将提高OBS2006-2测线地壳结构的可靠性和约束性,具有重要的研究意义.