Estimation of seismic spectral acceleration in Peninsular India

Peninsular India (PI), which lies south of 24°N latitude, has experienced several devastating earthquakes in the past. However, very few strong motion records are available for developing attenuation relations for ground acceleration, required by engineers to arrive at rational design response spectra for construction sites and cities in PI. Based on a well-known seismological model, the present paper statistically simulates ground motion in PI to arrive at an empirical relation for estimating 5% damped response spectra, as a function of magnitude and source to site distance, covering bedrock and soil conditions. The standard error in the proposed relationship is reported as a function of the frequency, for further use of the results in probabilistic seismic hazard analysis.     

利用神经网络技术在储层非均质性上的研究

储层非均质性对油气田的勘探和开发效果影响很大。结合测井资料二次解释的孔隙度(φ)、泥质含量(Vsh)、粒度中值(Md)、束缚水饱和度(Swi)等资料,借助BP网络算法,对渗透率进行预测。BP网络算法是人工神经网络算法中的一种,采用的是误差逆传算法,即把网络输出的误差归纳为各连接权的过程,通过把输出层的误差逐层向输入层逆向传播以分摊给各层单元,从而获取各层单元的参考误差,以便调整相应的联接权。如此反复,直至网络输出与期望输出之间的误差减少到允许误差为止。通过这样学习过程,消除了系统误差,提高了渗透率预测的精度,并为定性定量描述储层非均质性打下了坚实的基础。通过对萨尔图油田中部葡萄花油层组的非均质性进行了定量半定量的研究,认为萨尔图油田中部葡萄花油层组PI2小层层内非均质性最严重,与沉积微相密切相关。单砂层垂向上渗透率的变化以正韵律和复合式韵律为主,局部发育反韵律模式。经对葡萄花油层组PI1~PI4系统研究,层内非均质性强弱依次为PI2a、PI3、PI2b、PI1、PI4;层间非均质性在PI3和PI4间表现得最强,其他相对要弱;各小层平面非均质性相差无几,整体表现为较强的非均质性。     

Application of Wavelet Decomposition to Removing Barometric and Borehole Water Tidal Response in Level

Wavelet decomposition is used to analyze barometric fluctuation and earth tidal response in borehole water level changes. We apply wavelet analysis method to the decomposition of barometric fluctuation and earth tidal response into several temporal series in different frequency ranges. Barometric and tidal coefficients in different frequency ranges are computed with least squares method to remove barometric and tidal response. Comparing this method with general linear regression analysis method, we find wavelet analysis method can efficiently remove barometric and earth tidal response in borehole water level. Wavelet analysis method is based on wave theory and vibration theories. It not only considers the frequency characteristic of the observed data but also the temporal characteristic, and it can get barometric and tidal coefficients in different frequency ranges. This method has definite physical meaning.     

How chalcophile is rhenium? An experimental study of the solubility of Re in sulphide mattes

Although the Re/Os isotopic system has proved to be a valuable aid in understanding the evolution of the Earth's mantle, interpretation of the data is currently hindered by limited knowledge of the high-temperature geochemical behaviour of Re. In particular, the extent to which Re in the mantle is hosted by sulphide or silicate phases is poorly known. We report the results of an experimental study of the solubility of Re in sulphide melts (called here “mattes”) coexisting with a Re-rich Re–Fe alloy in the system Fe–Re–S–O over a range of fO₂, fS₂ and temperatures, which allow extrapolation to conditions pertinent to the Earth's mantle. The solubility of Re in mattes increases with increasing fS₂, with Re dissolving as Re⁴⁺ at high fS₂ and Re⁰ at low fS₂. The effect of fO₂ is negligible except at high fO₂ where O in the matte becomes important. At constant fS₂, an increase in temperature leads to an increase in the solubility of Re⁰ in the matte, but a decrease in the solubility of Re⁴⁺. These results, coupled with data for the Re solubility in silicate melts taken from the literature, allow the calculation of Re matte/silicate–melt partition coefficients (D_Re^matte/sil) for a range of conditions. The calculated D_Re^matte/sil show a large dependence on fO₂, and a lesser dependency on fS₂, with Re behaving as a chalcophile element for relatively reduced MORB-type mantle and as a lithophile element for the oxidised sources of island-arc basalts. The sensitivity of Re to fS₂ and fO₂ reconciles the apparent discrepancies between previous estimates of this parameter, which can vary by more than five orders of magnitude within the range of fS₂ and fO₂ covered by terrestrial basaltic magmas.     

Sphere-boundary edge and sphere-boundary corner contacts model in DDA for simulating particulate media in 3-D

In modelling particulate media, such as soils, using 3-D DDA, there are four types of contacts: sphere to sphere, sphere to boundary face, sphere to boundary edge, and sphere to boundary corner contacts. The first two were studied by the authors in a previous work (Beyabanaki and Bagtzoglou 2012). In this paper, we present a new contact model for sphere-boundary edge and sphere-boundary corner contacts in sphere-based 3-D DDA. The model includes a new algorithm to search for contacts, detect the contact types and calculate contact points. Moreover, formulas for contact sub-matrices are derived. The proposed contact model has been implemented into a sphere-based 3-D DDA program and three test cases are studied in order to verify the workability of the new contact model. The numerical results obtained demonstrate the capability of the model to deal with sphere-boundary interaction in particulate media.     

基于小波域高斯尺度混合模型的天气雷达图像高分辨率插值

采用可精确刻画雷达回波强度数据统计特征的小波域高斯尺度混合（GSM）模型作为雷达图像先验模型,进行天气雷达图像插值,在提高图像分辨率的同时有效重建降水回波中局部强回波值、小尺度变化细节等一些重要空间分布统计特征。分析和总结雷达回波强度数据小波频率域统计特点,建立小波域GSM模型;匹配天气雷达图像小波系数和GSM模型,利用贝叶斯理论估计更小尺度的小波系数,进行小波逆变换,完成高分辨率天气雷达图像插值。试验表明,该算法能从低分辨率图像中估计出高分辨率高频系数,且所利用的先验模型充分考虑降水数据本身的特点,可有效捕获降水回波结构的非高斯特征和局部相关特性,重建雷达图像中的局部变化细节。      

Bilge keel loads and hull pressures created by bilge keels fitted to a rotating cylinder

This paper presents bilge keel loads and hull pressure measurements carried out on a rotating cylinder in a free surface water basin. A flat plate bilge keel and one more complex shaped bilge keel were studied to investigate the geometry effect. The draft of the cylinder was varied to study the effect of the vicinity of the free surface on the bilge keel loads and hull pressures. The rotation axis of the cylinder was fixed to define a pure roll experiment (one degree of freedom).The cylinder was subject to forced oscillations of varying amplitude leading to a KC range of 0.3–16. Using Fourier analysis the first three harmonic coefficients representing the normal bilge keel load were derived. The first harmonic drag and inertia coefficients are in good agreement to existing experimental data obtained for wall bounded flat plates fitted in a U-shaped water tunnel as reported by Sarpkaya and O’Keefe (1996). New insight is gained by the fact that the addition of higher harmonic contributions is essential to capture the time varying bilge keel normal force.The pressure measurements next to the bilge keel are compared to measurements reported by Ikeda et al. (1979). Similar findings are obtained, showing that the pressure on the hull in front of the moving bilge keel is KC independent while the vortex system in the wake of the bilge keel leads to KC dependent hull pressure distributions. The hull pressure jump over the bilge keel correlates well to the force coefficient on the bilge keel. The complex nature of the vortex induced hull pressures is manifested. The empirically derived hull pressure distribution by Ikeda et al. (1979) for the time instant of maximum velocity is shown to correlate reasonably well to the measured data with some conservatism in the absolute value.Although a cylinder is very different from a ship-shaped section, the experiments provide essential insight into the physics associated with roll damping and into the factors that should be included in a roll damping prediction method.     

Wavelet‐based non‐stationary response analysis of a friction base‐isolated structure

A wavelet‐based stochastic formulation has been presented in this paper for the seismic analysis of a base‐isolated structural system which is modelled as a two‐degree‐of‐freedom (2‐DOF) system. The ground motion has been modelled as a non‐stationary process (both in amplitude and frequency) by using modified Littlewood–Paley basis wavelets. The proposed formulation is based on replacing the non‐linear system by an equivalent linear system with time‐dependent damping properties. The expressions of the instantaneous damping and the power spectral density function (PSDF) of the superstructure response have been obtained in terms of the functionals of input wavelet coefficients. The proposed formulation has been validated by simulating a ground motion process. The effect of the frequency non‐stationarity on the non‐linear response has also been studied in detail, and it has been clearly shown how ignoring the frequency non‐stationarity in the ground motion leads to inaccurate non‐linear response calculations. Copyright © 2000 John Wiley & Sons, Ltd.     

A new method to calculate the resonant response of a sunspot model atmosphere to magneto-atmospheric waves

In order to understand the observed oscillations in sunspots we present a new method for calculating the resonant response of a realistic semi-empirical model of the sunspot umbral atmosphere and subphotosphere to magneto-atmospheric waves in a vertical magnetic field. The depth dependence of both the adiabatic coefficient and the turbulent pressure is taken into account. This requires an extension of the wave equations by Ferraro & Plumpton (1958). We compare the coefficients of wave transmission, re flection, and conversion between fast mode and slow mode waves for different assumptions, compare the results with those from earlier modelling efforts, and point out possible sources of mistakes. The depth dependence of the adiabatic coefficient strongly influences the resulting spectrum of resonance frequencies. The condition of a conservation of wave flux is violated if the depth dependence of the turbulent pressure is not properly considered.     

计算S型海洋重力仪交叉耦合改正的测线系数修正法

提出了测线系数修正法,即在交叉点处采用平差的思想对每条测线的交叉耦合(CC)改正监视项的系数进行修正,然后采用新系数重新计算CC改正。实验结果表明,采用修正后的系数能补偿外界的动态测量条件与仪器生产厂家计算CC改正系数时设计的假设条件之间的差异对CC改正造成的影响,提高海洋重力数据处理的精度。