基于二级近似离散复镜像法的低频电磁格林函数计算

本文讨论了利用二级近似离散复镜像法实现低频格林函数快速、精确的计算.通过数值分析,给出了在低频电磁场计算中该方法近似参数的选取原则.(1)二级近似的积分区间里谱格林函数采样个数N为十倍的近似多项式的个数M;(2)总的积分区间L2大小约为40/r,r为收发距离;(3)两个积分区间的分界值L1为总的积分区间大小L2与第二个...     

National-scale data-driven rainfall induced landslide susceptibility mapping for China by accounting for incomplete landslide data

China is one of the countries where landslides caused the most fatalities in the last decades.The threat that landslide disasters pose to people might even be greater in the future,due to climate change and the increasing urbanization of mountainous areas.A reliable national-scale rainfall induced landslide suscep-tibility model is therefore of great relevance in order to identify regions more and less prone to landslid-ing as well as to develop suitable risk mitigating strategies.However,relying on imperfect landslide data is inevitable when modelling landslide susceptibility for such a large research area.The purpose of this study is to investigate the influence of incomplete landslide data on national scale statistical landslide susceptibility modeling for China.In this context,it is aimed to explore the benefit of mixed effects mod-elling to counterbalance associated bias propagations.Six influencing factors including lithology,slope,soil moisture index,mean annual precipitation,land use and geological environment regions were selected based on an initial exploratory data analysis.Three sets of influencing variables were designed to represent different solutions to deal with spatially incomplete landslide information:Set 1(disregards the presence of incomplete landslide information),Set 2(excludes factors related to the incompleteness of landslide data),Set 3(accounts for factors related to the incompleteness via random effects).The vari-able sets were then introduced in a generalized additive model(GAM:Set 1 and Set 2)and a generalized additive mixed effect model(GAMM:Set 3)to establish three national-scale statistical landslide suscep-tibility models:models 1,2 and 3.The models were evaluated using the area under the receiver operating characteristics curve(AUROC)given by spatially explicit and non-spatial cross-validation.The spatial pre-diction pattern produced by the models were also investigated.The results show that the landslide inven-tory incompleteness had a substantial impact on the outcomes of the statistical landslide susceptibility models.The cross-validation results provided evidence that the three established models performed well to predict model-independent landslide information with median AUROCs ranging from 0.8 to 0.9.However,although Model 1 reached the highest AUROCs within non-spatial cross-validation(median of 0.9),it was not associated with the most plausible representation of landslide susceptibility.The Model 1 modelling results were inconsistent with geomorphological process knowledge and reflected a large extent the underlying data bias.The Model 2 susceptibility maps provided a less biased picture of landslide susceptibility.However,a lower predicted likelihood of landslide occurrence still existed in areas known to be underrepresented in terms of landslide data(e.g.,the Kuenlun Mountains in the northern Tibetan Plateau).The non-linear mixed-effects model(Model 3)reduced the impact of these biases best by introducing bias-describing variables as random effects.Among the three models,Model 3 was selected as the best national-scale susceptibility model for China as it produced the most plausible portray of rainfall induced landslide susceptibility and the highest spatially explicit predictive perfor-mance(median AUROC of spatial cross validation 0.84)compared to the other two models(median AUROCs of 0.81 and 0.79,respectively).We conclude that ignoring landslide inventory-based incomplete-ness can entail misleading modelling results and that the application of non-linear mixed-effect models can reduce the propagation of such biases into the final results for very large areas.     

Genetic waveform modeling for the crustal structure in Northeast Japan

We propose a genetic algorithm (GA) search procedure for waveform modeling of local crustal earthquakes for optimal one-dimensional (1-D) crustal velocity model. Both waveforms and travel-time data are used for the structure determination. The use of travel times in model evaluation improves the waveform modeling performance in the sense of computation speed and accuracy. We applied this method to broadband waveforms of a local crustal earthquake (M 4.2) in Northeast Japan. P-wave velocities of the crustal model are found to be 4.95 ± 0.30, 5.9 ± 0.02, and 6.51 ± 0.20 km/s for a surface layer, upper crust and lower crust, respectively. The surface layer thickness and the Conrad and Moho depths are found to be 3.01 ± 0.8, 17.77 ± 0.4 and 34.59 ± 1.0 km, respectively. For epicentral distances <200 km, our synthetic waveforms match the observed ones generally well. Early arrivals are mainly observed at stations near the Pacific coast in the forearc area having a thinner crust. In contrast, delayed arrivals appear at stations near the volcanic front and back-arc areas where low-velocity anomalies exist due to the effect of the Pacific slab dehydration and the hot upwelling flows in the mantle wedge. In general, our results agree well with the main tectonic setting of the study area, which confirms the reliability of the proposed approach. Despite a 1-D velocity model is too simple to represent the complex crustal structure, it is still required for the conventional routine analysis of seismology, such as earthquake location and source parameter studies. The current approach is considered as a step toward the genetic full waveform modeling for the 3-D velocity model estimation.     

标准化降水指数与广义极值分布干旱指数在西北地区应用的对比分析

从气象干旱的角度出发,对比分析了基于广义极值分布理论构建的新干旱指数GEVI(Generalized Extreme Value Index)和我国广泛应用的SPI (Standardized Precipitation Index)指数,详细讨论了SPI和GEVI指数在西北地区干旱监测中的适用性.结果表明,SPI和GEVI指数在计算原理上具有相似性,在一定程度上衡量某一地区干旱程度和反映干旱范围方面具有一致性和可替代性.但因两者采取不同的分布形态,对干旱等级的衡量效果不尽相同.GEVI指数在拟合降水量分布函数时更为详尽、客观.在西北5省(区)的干旱监测中,对地处干旱半干旱区的甘肃来说,GEVI指数在于旱强度判定方面优于SPI指数;对年降水量最少的新疆地区来说,GEVI指数比SPI指数在干旱监测方面更为适用.     

Numerical assessment of liquefaction mitigation effects on residential houses: Case histories of the 2007 Niigata Chuetsu-offshore earthquake

In this study, the effectiveness of liquefaction countermeasures for residential houses was explored using a fully coupled dynamic effective-stress finite element procedure. Numerical analyses were conducted on two wooden houses that were damaged to different degrees due to dune liquefaction during the 2007 Niigata Chuetsu-offshore earthquake. House A, which was only improved by horizontal drainage pipes to lower the ground water level, was completely destroyed; however, house B, which was improved by a horizontal drainage system, soil–cement mixtures, and steel-pipe piles, was slightly deformed. Numerical results show that the effects of the sand dune slope on the damage to the two houses were somewhat different. For house B, it was found that the steel-pipe piles were more effective. Two countermeasures generally led to a greater degree of reduction in both lateral and vertical displacements of house B than only a single countermeasure employed. In addition, the combined implementation of steel-pipe piles and soil–cement mixtures was the most effective among the cases with two countermeasures.     

Effects of sample size on the accuracy of geomorphological models

Commonly, the most costly part of geomorphological distribution modelling studies is gathering the data. Thus, guidance for researchers concerning the quantity of field data needed would be extremely practical. This paper scrutinises the relationship between the sample size (the number of observations varied from 20 to 600) and the predictive ability of the generalized linear model (GLM), generalized additive model (GAM), generalized boosting method (GBM) and artificial neural network (ANN) in two data settings, i.e., independent and split-sample approaches. The study was performed using empirical data of periglacial processes from an area of 600 km² in northernmost Finland at grid resolutions of 1 ha (100 × 100 m) and 25 ha (500 × 500 m). A rather sharp increase in the predictive ability of the models was observed when the number of observations increased from 20 to 100, and the level of robust predictions was reached with 200 observations. The result indicates that no more than a few hundred observations are needed in geomorphological distribution modelling at a medium scale resolution (ca. 0.01–1 km²).     

时频域油气储层低频阴影检测

为了准确刻画地震信号的局部层次结构,实现高效率的三维地震资料瞬时谱分解,检测油气储层的低频阴影,构造了广义S变换.广义S变换通过引入两个参数,改造S变换的小波函数,使其小波能根据信号处理的具体应用需要而调整.模型信号的仿真发现,广义S变换具有更加优越和灵活可调的时频聚集性能.文中分析了低频阴影的机理,并利用广义S变换对三维实际地震资料进行了瞬时谱分解,它不仅能检测油气储层的低频阴影,而且可以刻画油气储层的岩性边界和空间展布,减小油气储层检测的多解性.     

三维层状孔隙介质中弹性波的一种积分表达式I：理论

孔隙介质弹性波传播理论在地球物理勘探、地震工程和岩土动力学等领域有着广泛的应用.而孔隙介质中的弹性波受孔隙度、渗透率、流体黏滞系数等参数的影响,因此研究波场的传播特征将有助于分析和提取这些信息.本文在Biot理论的基础上,针对三维层状孔隙介质模型,利用在合成理论地震图的研究中已经被证实具有稳定、高效且适用范围较广的Luco-Apsel-Chen（LAC）广义反透射方法,给出了弹性波场的一种积分形式的半解析解,可通过数值方法高效、准确地计算层状孔隙介质中的理论波场,所以该积分形式的半解析解可为三维层状孔隙介质波场传播特征的理论数值模拟研究提供一种新的途径和手段.     

Generalized cyclic p–y curve modeling for analysis of laterally loaded piles

Owing to their simplicity and reasonable accuracy, Beam on Nonlinear Winkler Foundation (BNWF) models are widely used for the analysis of laterally loaded piles. Their main drawback is idealizing the soil continuum with discrete uncoupled springs representing the soil reactions to pile movement. Static p–y curves, obtained from limited full-scaled field tests, are generally used as a backbone curve of the model. However, these empirically derived p–y curves could not incorporate the effects of various pile properties and soil continuity. The strain wedge method (SWM) has been improved to assess the nonlinear p–y curve response of laterally loaded piles based on a three-dimensional soil–pile interaction through a passive wedge developed in front of the pile. In this paper, the SWM based p–y curve is implemented as the backbone curves of developed BNWF model to study the nonlinear response of single pile under cyclic lateral loading. The developed nonlinear model is capable of accounting for various important soil–pile interaction response features such as soil and pile yielding, cyclic degradation of soil stiffness and strength under generalized loading, soil–pile gap formation with soil cave-in and recompression, and energy dissipation. Some experimental tests are studied to verify the BNWF model and examine the effect of each factor on the response of laterally loaded pile embedded in sand and clay. The experimental data and computed results agree well, confirming the model ability to predict the response of piles under one-way and two-way cyclic loading. The results show that the developed model can satisfactorily simulate the pile stiffness hardening due to soil cave in and sand densification as observed in the experiment. It is also concluded from the results that the gap formation and soil degradation have significant effects on the increase of lateral pile-head deflection and maximum bending moment of the pile in cohesive soils.     

广义非线性动态最小二乘问题的一个直接解算方法

构建“数字地球”、“数字国家”等数字化科学工程的基础是数据[1 ] ,其数据具有多源、多维、多种类型、多种时态、多种精度并具有非线性特征等特点[2 ] ,首先要进行数据处理并应采用全新的广义非线性动态最小二乘法[3] [4] ,数据处理方法的核心是广义非线性动态最小二乘问题参数估计的函数模型及其解算方法 ,迄今国内外对这方面的研究尚不多。本文在作者前期研究、提出的广义非线性动态最小二乘函数模型参数估计迭代法求解[5] 的基础上 ,进一步研究、提出了一种广义非线性动态最小二乘模型参数估计的直接解算方法 ,将问题分离 ,把待求参数减半 ,直接求解。从而大大降低求解问题的维数 ,大大减少计算难度和计算工作量 ,这是国内外首次研究提出的一种比迭代法更快速、更有效、更科学的解算方法。为多源、多类、多时态数据处理开辟了一新途径 ,也大大扩大了广义非线性动态最小二乘法的应用面