A novel coupled FDM-DEM modelling method for flexible membrane boundary in laboratory tests

Difficulties are involved in discrete element method (DEM) modelling of the flexible boundary, that is, the membranes covering the soil sample, which can be commonly found in contemporary laboratory soil tests. In this paper, a novel method is proposed wherein the finite difference method (FDM) and DEM are coupled to simulate the rubber membrane and soil body, respectively. Numerical plane strain and triaxial tests, served by the flexible membrane, are implemented and analysed later. The effect of the membrane modulus on the measurement accuracy is considered, with analytical formulae derived to judge the significance of this effect. Based on an analysis of stress-strain responses and the grain rotation field, the mechanical performances produced by the flexible and rigid lateral boundaries are compared for the plane strain test. The results show that (1) the effect of the membrane on the test result becomes more significant at larger strain level because the membrane applies additional lateral confining pressure to the soil body; (2) the tested models reproduce typical stress and volumetric paths for specimens with shear bands; (3) for the plane strain test, the rigid lateral boundary derives a much higher peak strength and larger bulk dilatation, but a similar residual strength, compared with the flexible boundary. The latter produces a more uniform (or ‘diffuse') rotation field and more mobilised local kinematics than does the former. All simulations show that the proposed FDM-DEM coupling method is able to simulate laboratory tests with a flexible boundary membrane.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

地下洞库勘察中的深斜孔钻探和压水试验技术

青岛市黄岛区拟建设海底地下油气储藏库,前期地质勘察工作要求施工倾角45°、深220m的取心钻孔4个,并在孔内分段进行压水试验.本文介绍了如何利用XY-4型金刚石小口径钻机施工基岩深斜孔及压水试验技术.     

变形分析的Bayes估计与检验

由于观测条件的变化，观测值方差将产生随机波动。本文根据Ｂａｙｅｓ估计理论，推导了线性模型待估参数的分布函数及观测值方差估计公式，并将其用于位移检验。计算结果表明，在同样的置信度下，本方法更符合实际。文中还对参数估值的统计性质进行了讨论。     

塔中4号油藏渗透率分布的非均质性及渗流方向研究

研究油藏内部油水渗流方向时,主要考虑渗透率影响,引入两个水平方向渗透率比、水平与垂向渗透率比参数来分析。前者反映样品所在储层平面上两个水平方向渗流能力的差异,后者反映样品所在储层水平与垂向渗流能力的差异,并经小柱样测试和全岩心测试两种方法确定两个参数。通过对塔中4号油藏CⅠ油组9口井2306个样品两种渗透率比值分布直方图分析,研究了CⅠ油组渗透率的非均质性及其渗流方向:CI油组储层在水平方向的渗流能力有很大差异,具有明显的方向性。储层垂向渗流能力明显小于水平方向的渗流能力。     

密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

Predicting ground-motion variations at the Turkey-flat test site,California

Results of a single group participating in an international experiment are analyzed. The experiment served to verify computational predictions of the ground-motion variations due to near-surface geological effects at a site established for that purpose by the California Department of Conservation. Based on an acceleration record at a rock location, and geotechnical model of medium, records at the other locations of a nearby sedimentary deposit were predicted. A 2-D finite-difference sensitivity analysis suggested that the lateral wave-propagation effects are negligibly small, and locally 1-D computations are sufficient for the present site. Those computations are compared with observations not available to the authors during the blind prediction. Peak accelerations, peak velocities and RMS accelerations were predicted with errors less than 159%, 114% and 62%, respectively. Maxima of the response spectra were fitted within a factor of 2. The predicted and observed Husid's plots (i.e., the normalized cumulative plots of the acceleration squared) have the correlation coefficients 0.98. The detected misfits do not show any simple relation to the instrument location, component, frequency, or time.     

宁夏香山-天景山弧形断裂带新活动特征及1709年中卫南71/2级地震形变带

本文根据野外地质填图和水系位移测量结果,论述了香山-天景山弧形断裂带新生代有两个不同活动性质的阶段,即早期阶段的强烈挤压和晚期阶段的左旋走滑兼挤压。分析、讨论了不同活动阶段的时间界限和转变原因。指出了1709年中卫南71/_2级地震形变带的表现形式、延伸范围     

Crustal structure at the easternmost termination of the Variscan belt based on CELEBRATION 2000 and ALP 2002 data

The eastern margin of the Variscan belt in Europe comprises plate boundaries between continental blocks and terranes formed during different tectonic events. The crustal structure of that complicated area was studied using the data of the international refraction experiments CELEBRATION 2000 and ALP 2002. The seismic data were acquired along SW–NE oriented refraction and wide-angle reflection profiles CEL10 and ALP04 starting in the Eastern Alps, passing through the Moravo-Silesian zone of the Bohemian Massif and the Fore-Sudetic Monocline, and terminating in the TESZ in Poland. The data were interpreted by seismic tomographic inversion and by 2-D trial-and-error forward modelling of the P waves. Velocity models determine different types of the crust–mantle transition, reflecting variable crustal thickness and delimiting contacts of tectonic units in depth. In the Alpine area, few km thick LVZ with the Vp of 5.1 km s^− 1 dipping to the SW and outcropping at the surface represents the Molasse and Helvetic Flysch sediments overthrust by the Northern Calcareous Alps with higher velocities. In the Bohemian Massif, lower velocities in the range of 5.0–5.6 km s^− 1 down to a depth of 5 km might represent the SE termination of the Elbe Fault Zone. The Fore-Sudetic Monocline and the TESZ are covered by sediments with the velocities in the range of 3.6–5.5 km s^− 1 to the maximum depth of 15 km beneath the Mid-Polish Trough. The Moho in the Eastern Alps is dipping to the SW reaching the depth of 43–45 km. The lower crust at the eastern margin of the Bohemian Massif is characterized by elevated velocities and high Vp gradient, which seems to be a characteristic feature of the Moravo-Silesian. Slightly different properties in the Moravian and Silesian units might be attributed to varying distances of the profile from the Moldanubian Thrust front as well as a different type of contact of the Brunia with the Moldanubian and its northern root sector. The Moho beneath the Fore-Sudetic Monocline is the most pronounced and is interpreted as the first-order discontinuity at a depth of 30 km.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.