Seismic response of flexible cantilever retaining walls with dry backfill

Failure of retaining walls during earthquakes has occurred many times in the past. Although significant progress has been made in analysing the seismic response of rigid gravity type retaining walls, considerable difficulties still exist in the seismic-resistant design of the flexible cantilever type of retaining walls because of the complex nature of the dynamic soil–structure interaction. In this paper the seismic response of cantilever retaining walls with dry backfill is simulated using centrifuge modelling and numerical modelling. It is found that bending moments on the wall increased significantly during an earthquake. After the end of base shaking, the residual moment on the wall was significantly higher than the moment under static loading. The numerical simulation is able to model quite accurately the main characteristics of acceleration, bending moment, and displacement recorded in the centrifuge test.     

Development of spectral hazard maps for a proposed revision of the Indonesian Seismic Building Code

The need to revise the current Indonesian Seismic Hazard Map contained in Indonesian Earthquake Resistant Building Code SNI 03-1726-2002 which partially adopts the concept of UBC 1997, was driven among others by the desire to better reflect the potential larger earthquake disasters faced by the nation in the future. The much larger than maximum predicted Aceh Earthquake (M _w 9.0–9.3) of 2004, followed by the destruction observed during the 2005 Nias Earthquake (M _w 8.7) urgently underline to need to consider the new conceptual approach and technological shift shown in the transition of UBC 1997 to IBC 2006. This paper presents research works for developing spectral hazard maps for Indonesia. Some improvements in seismic hazard analysis were implemented using recent seismic records. Seismic sources were modeled by background, fault, and subduction zones by considering a truncated exponential model, pure characteristic model or both models. A logic tree method was performed to account for the epistemic uncertainty and several attenuation functions were selected. Maps of PGA and spectral accelerations for a short period (0.2 s) and for a 1-s period were then developed using a probabilistic approach. The maps will be proposed as a revision for the current seismic hazard map in the Indonesian Seismic Building Code.     

Characterisation of carbonate minerals from hyperspectral TIR scanning using features at 14 000 and 11 300 nm

Rapid characterisation of carbonate phases in hyperspectral reflectance spectra acquired from drill core material has important implications for mineral exploration and resource modelling. Major infrared active features of carbonates lie in the thermal region around 6500 nm, 11 300 nm and 14 000 nm, with the latter two features being most useful for differentiating mineral species. A scatter diagram of the wavelength of the 14 000 nm feature vs that of the 11 300 nm feature, powerfully differentiates carbonates. Although the wavelength of the 11 300 nm peak is easily measured, the 14 000 nm trough and peak are commonly weak and their wavelengths can confidently be used only after filtering the spectra, e.g. selecting only those with the trough and peak separated by 175–230 nm, typical of common carbonates. The method is demonstrated with drillhole 120R from the Rosebery polymetallic VHMS deposit in western Tasmania, which has been scanned with the HyLogger-3 system. A 14 000–11 300 plot shows a high degree of clustering of the drillhole 120R data close to the library spectra of calcite, dolomite, Fe-dolomite, ankerite, kutnohorite, rhodochrosite, Fe-rhodochrosite and siderite. The interpreted compositions of the carbonate spectral populations strongly correlate with the chemical populations of 144 analysed carbonates and provide a highly resolved spatial framework for interpreting carbonate alteration.     

Kappa (k) derived from accelerograms recorded in the 2008 Wenchuan mainshock,Sichuan, China

High-frequency spectral decay factor, kappa (k), in the accelerograms of the Wenchuan mainshock was measured using strong motion data from 52 stations within 311 km of the epicenter. The derived k range from 0.0034 s to 0.0468 s. The correlation of k versus fault distance was given, which is k = 0.01288 + 5.9068 × 10^–5 R for the N-S component, k = 0.01881 + 1.4219 × 10^–5 R for the E-W component, and k = 0.00855 + 5.6086 × 10^–5 R for the U-D component. The analysis on the spatial variation of k demonstrates that k relates to source effect and propagation effect besides local site effect. Ground motions for the 52 stations were simulated using derived k and compared to actual recordings in terms of waveforms, amplitude spectra and response spectra. The results show agreement at shorter periods (<1 s), but a slight overestimation at longer periods (1–7 s).     

判断水淹程度的碳氧比能谱测井方法研究

碳氧比能谱测井方法是在套管井中用于划分油水层、确定含油水饱和度及岩性的一种有效方法。根据碳氧比能谱测井的原始谱线资料,提出了一种判断油层水淹程度的指数——水淹指数。通过实际资料验证,该指数对水淹程度有很好的判断作用,具有良好的应用前景。     

The Fourier Integral Method: An efficient spectral method for simulation of random fields

The Fourier Integral Method (FIM) of spectral simulation, adapted to generate realizations of a random function in one, two, or three dimensions, is shown to be an efficient technique of non-conditional geostatistical simulation. The main contribution is the use of the fast Fourier transform for both numerical calculus of the density spectral function and as generator of random finite multidimensional sequences with imposed covariance. Results obtained with the FIM are compared with those obtained by other classic methods: Shinozuka and Jan Method in 1D and Turning Bands Method in 2D and 3D, the points for and against different methodologies are discussed. Moreover, with the FIM the simulation of nested structures, one of which can be a nugget effect and the simulation of both zonal and geometric anisotropy is straightforward. All steps taken to implement the FIM methodology are discussed.     

Intra-seasonal variations of kinetic energy of lower tropospheric zonal waves during northern summer monsoon

Space spectral analysis of zonal (u) and meridional (v) components of wind and time spectral analysis of kinetic energy of zonal waves at 850 hPa during monsoon 1991 (1st June 1991 to 31st August 1991) for the global belt between equator and 40°N are investigated. Space spectral analysis shows that long waves (wavenumbers 1 and 2) dominate the energetics of Region 1 (equator to 20°N) while over Region 2 (20°N to 40°N) the kinetic energy of short waves (wavenumbers 3 to 10) is more than kinetic energy of long waves. It has been found that kinetic energy of long waves is dominated by zonal component while both (zonal and meridional) the components of wind have almost equal contribution in the kinetic energy of short waves. Temporal variations of kinetic energy of wavenumber 2 over Region 1 and Region 2 are almost identical. The correlation matrix of different time series shows that (i) wavenumber 2 over Regions 1 and 2 might have the same energy source and (ii) there is a possibility of an exchange of kinetic energy between wavenumber 1 over Region 1 and short waves over Region 2. Wave to wave interactions indicate that short waves over Region 2 are the common source of kinetic energy to wavenumber 2 over Regions 1 and 2 and wavenumber 1 over Region 1. Time spectral analysis of kinetic energy of zonal waves indicates that wavenumber 1 is dominated by 30–45 day and bi-weekly oscillations while short waves are dominated by weekly and bi-weekly oscillations. The correlation matrix, wave to wave interaction and time spectral analysis together suggest that short period oscillations of kinetic energy of wavenumber 1 might be one of the factors causing dominant weekly (5–9 day) and bi-weekly (10–18 day) oscillations in the kinetic energy of short waves.     

Scale Matching with Factorial Kriging for Improved Porosity Estimation from Seismic Data

When seismic data and porosity well logs contain information at different spatial scales, it is important to do a scale-matching of the datasets. Combining different data types with scale mismatch can lead to suboptimal results. A good correlation between seismic velocity and rock properties provides a basis for integrating seismic data in the estimation of petrophysical properties. Three-dimensional seismic data provides an unique exhaustive coverage of the interwell reservoir region not available from well data. However, because of the limitations of measurement frequency bandwidth and view angles, the seismic image can not capture the true seismic velocity at all spatial scales present in the earth. The small-scale spatial structure of heterogeneities may be absent in the measured seismic data. In order to take best advantage of the seismic data, factorial kriging is applied to separate the small and large-scale structures of both porosity and seismic data. Then the spatial structures in seismic data which are poorly correlated with porosity are filtered out prior to integrating seismic data into porosity estimation.     

基于遥感视反射率图像的矿化信息识别及其应用

成矿作用过程中产生的矿化蚀变,在遥感图像上往往是一种非直接的弱相关信息。文章以川西西支沟地区为例,通过典型矿床地面波谱分析,并在视反射率图像模拟及其与典型矿床地面波谱反射率的相关分析基础上,对矿化蚀变信息进行反演识别研究。西支沟地区的研究表明,基于视反射率图像的矿化信息反演识别分析可增强遥感图像矿化信息,在多金属矿化信息识别方面有较理想的找矿预测效果。     

频谱分解技术在泥页岩储层含气性预测中的应用——以柴达木盆地鱼卡凹陷为例

针对柴达木盆地北缘鱼卡凹陷中侏罗统大煤沟组泥页岩,利用地震频谱曲线在不同流体介质储层中的形态、斜率、主频与中心频率的差值等属性参数特征,识别陆相泥页岩层系的地震频谱属性特征,预测含气性。研究结果表明,中侏罗统大煤沟组上部能量系数较低,平均为0.6719,下部能量系数较高,平均为0.7521,结合现场含气量和全烃特征认为,可以以能量系数0.7为界线,作为识别下部和上部含气性的标志,下部明显比上部含气性好。平面上进一步圈定了2个有利目标区,得到柴页1井勘探验证。因此,频谱分解技术可应用于中国陆相页岩气的勘查开发。