Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

An efficient solution of real-time data processing for multi-GNSS network

Global navigation satellite systems (GNSS) are acting as an indispensable tool for geodetic research and global monitoring of the Earth, and they have been rapidly developed over the past few years with abundant GNSS networks, modern constellations, and significant improvement in mathematic models of data processing. However, due to the increasing number of satellites and stations, the computational efficiency becomes a key issue and it could hamper the further development of GNSS applications. In this contribution, this problem is overcome from the aspects of both dense linear algebra algorithms and GNSS processing strategy. First, in order to fully explore the power of modern microprocessors, the square root information filter solution based on the blocked QR factorization employing as many matrix–matrix operations as possible is introduced. In addition, the algorithm complexity of GNSS data processing is further decreased by centralizing the carrier-phase observations and ambiguity parameters, as well as performing the real-time ambiguity resolution and elimination. Based on the QR factorization of the simulated matrix, we can conclude that compared to unblocked QR factorization, the blocked QR factorization can greatly improve processing efficiency with a magnitude of nearly two orders on a personal computer with four 3.30 GHz cores. Then, with 82 globally distributed stations, the processing efficiency is further validated in multi-GNSS (GPS/BDS/Galileo) satellite clock estimation. The results suggest that it will take about 31.38 s per epoch for the unblocked method. While, without any loss of accuracy, it only takes 0.50 and 0.31 s for our new algorithm per epoch for float and fixed clock solutions, respectively.     

Three dimensional large scale earthquake simulator with six degrees of freedom

Conclusions When the aseismic test of large scale structure is performed, it is a major performance that the actual response waveform is identical to the input signal waveform under the condition of analog control, especially when the damping of the structure tested is small. The rotational motions of the simulator are induced due to the structure inertia, so that it is a key technique for the simulator design that realizes independent control of the simulator with six degrees of freedom and suppresses the rotational motions. In the design, the modern control techniques are applied and a high performance of the three dimensional earthquake simulator is developed successfully. The simulator can simulate actual earthquake load and it is an important test device for aseismic research. This paper is a main project of the Eighth Five Plan, State Seismological Bureau.     

Characteristics of particle size for creeping and saltating sand grains in aeolian transport

Creep and saltation are the primary modes of surface transport involved in the fluid‐like movement of aeolian sands. Although numerous studies have focused on saltation, few studies have focused on creep, primarily because of the experimental difficulty and the limited amount of theoretical information available on this process. Grain size and its distribution characteristics are key controls on the modes of sand movement and their transport masses. Based on a series of wind tunnel experiments, this paper presents new data regarding the saltation flux, obtained using a flat sampler, and on the creeping mass, obtained using a specifically designed bed trap, associated with four friction velocities (0·41, 0·47, 0·55 and 0·61 m sec^?1). These data yielded information regarding creeping and saltating sand grains and their particle size characteristics at various heights, which led to the following conclusions: (i) the creeping masses increased as a power function (q = ?1·02 + 14·19u_*³) of friction wind velocities, with a correlation (R²) of 0·95; (ii) the flux of aeolian sand flow decreases exponentially with increasing height (q = a exp(–z/b)) and increases as a power function (q = ?26·30 + 428·40 u_*³) of the friction wind velocity; (iii) the particle size of creeping sand grains is ca 1·15 times of the mean diameter of salting sand grains at a height of 0 to 2 cm, which is 1·14 times of the mean diameter of sand grains in a bed; and (iv) the mean diameter of saltating sand grains decreases rapidly with increasing height whereas, while at a given height, the mean diameter of saltating sand grains is positively correlated with the friction wind velocity. Although these results require additional experimental validation, they provide new information for modelling of aeolian sand transport processes.     

聚合物驱组合段塞分注实验研究

萨尔图油田南二区东部140号断层以东地区葡11—4油层,以多段多韵律和正韵律沉积为主,平均孔隙度30．0％,平均空气渗透率1980×10^-3μm2。水驱结束时该油层平均含水92％,采出程度39．5％。计划对该区块实施聚合物驱开采。以往经验认为,组合段塞调整注聚技术能够有效改善流度控制作用。然而对组合注入情况下如何进一步提高驱油效果、注入井分注时机、层段注入强度对开发效果有何影响却很少探讨。通过室内驱油实验,在聚驱的组合段塞注入不同时机（水驱空白阶段、含水下降阶段、含水稳定阶段、含水回升阶段）下调整注入强度实施分注以及在含水回升分注基础上进一步调整注入强度,同时对比不分注实验,评价各方案驱油效果。实验表明：不实施分注时聚驱采收率在水驱基础上提高了19．07个百分点,而组合段塞注入时实施分注均能够在不分注基础上进一步提高驱油效果;空白水驱时实施分注采收率提高幅度最大,为26．55个百分点;随着分注时机的滞后采收率提高幅度增加值依次递减,分注时机越早越好。     

覆盖中国沿海地区的精细化台风风暴潮模型的研究及适用

精细化风暴潮预报是目前风暴潮预报重点发展方向之一,本文首次建立起了一个覆盖整个中国沿海地区的精细化台风风暴潮数值模型,克服了以往分区域数值模型的不足,该模型在中国沿海地区的分辨率达到300m左右。模型采用了并行计算,并对2012年和2013年灾害性台风风暴潮过程进行了数值检验,计算精度和计算所用时间都能够满足业务化运行的要求。本文同时还根据中国气象局、美国国家气象局等5家主要台风预报机构给出的24h台风预报,对2013年度灾害性台风风暴潮过程进行了24h数值预报检验,检验结果表明:根据中国气象局台风登陆前24h预报可以得到更准确的风暴潮预报结果,其预报结果优于其他各家预报结果。该结论可以为今后的台风风暴潮预报中台风路径的选取提供重要的参考。     

Utilization of water resources,ecological balance and land desertification in the Tarim Basin,Xinjiang

Because of the human exploitation and utilization of water resources in the Tarim Basin, the water resources consumption has changed from mainly natural ecosystem to artificial oasis ecosystem, and the environment has changed correspondingly. The basic changes are: desertification and oasis development coexist, both “the human being advance and the desert retreat” and “the desert advance and the human being retreat” coexist, but the latter is dominant. In the upper reaches, water volume drawing to irrigated agricultural areas has increased, artificial oases have been enlarging and moving from the deltas in the lower reaches of many rivers to the piedmont plains. In the middle and lower reaches of the Tarim River, the stream flow has decreased, old oases have declined, natural vegetations have been degenerating, desertification has been enlarging, and the environment has deteriorated. The transition regions, which consist of forestlands, grasslands and waters between the desert and the oases, have been decreasing continuously, their shelter function to the oases has been weakened, and the desert is threatening the oases seriously.     

大回线源内三维地形对电磁响应影响研究

我国地质地形条件复杂多样,地形效应研究与影响校正是电磁勘探数据处理与解释中的重要问题,忽略地形影响将造成数据解释的较大误差.为了避免复杂的地形效应,实际勘探中尽量将发射源与接收点布置在地形平缓处,但源内仍然可能存在水坑、空洞等小型地形对电磁响应造成一定影响.本文基于有限元法实现了频域、时间域电磁法三维数值模拟,构建了大回线源内含凹陷地形的多种勘探场景,分析了地形尺寸、填充介质变化对电磁响应的影响规律.文中通过与均匀半空间解析解、前人的块状高导体模型的数值结果对比,验证了基于有限元法的频域、时域电磁法三维正演的正确性和精度.在矩形回线源内构造了不同尺寸的凹陷区域,将凹陷区域填充空气、雨水介质模拟空洞和水坑,以均匀半空间响应为参考,定量考察了空洞、水坑对电磁响应的影响程度.研究表明,回线源内地形对电磁响应的影响主要反映在较高频段和时间域早期,水坑对电磁响应的影响显著高于空洞,在实际勘探中,应尽量避开水坑,而存在较小空洞时,可忽略地形效应.