单面边坡的两种动力反应形式及其临界高度

把单面边坡的动力反应分为高边坡动力反应和低边坡动力反应两类. 研究了两种不同动力反应形式临界边坡高度的影响因素. 量纲分析表明弹性条件下边坡动力反应临界高度Hthre与构成边坡材料的动弹性模量Ed的平方根成正比,与边坡材料的密度ρ成反比,与边坡动力输入的特征周期T成正比. 结合弹性波动理论,利用大量数值模拟和统计分析得到了Hthre的表达式. 当边坡的高度大于Hthre时,边坡的动力反应呈现高边坡动力反应规律,当边坡的高度小于Hthre时,边坡的动力反应呈现低边坡动力反应规律.     

A quantitative study of some aftershock sequences in greece

Temporal features of the aftershock activities following twelve moderate shallow earthquakes in Greece have been studied quantitatively, by making use of the modified Omori's formula and Akaike's Information Criterion (AIC).Although the earthquake sequences occurred in different seismotectonic regions in Greece, similar seismic patterns before the occurrence of large aftershocks have been observed. Aftershock activity shows an appreciable decrease from the level expected from the modified Omori formula before the occurrence of the large aftershock. Also in some cases, the aftershock activity recovers to a normal level or increases beyond it prior to the occurrence of the large aftershock, if the aftershock activity is monitored immediately after the occurrence of the main shock. In such case, the observed pattern may be useful in predicting a large aftershock.     

Apparent resistivity of azimuthal anisotropy layered media

Introduction The method of apparent resistivity with direct current (DC) measurement, adopted from geophysical exploration, is one of the most important precursory approaches in earthquake prediction, in which the geo-electric constitutive model is illustrated by the scalar theory Archies law and correspondingly a set of calculation methods for apparent resistivity of layered media has been established (QIAN, et al, 1985; YAO, 1989). But this is not consistent with the fact that in earthqu…     

层状方位各向异性介质的视电阻率计算

从电性各向异性的欧姆定律出发，推导了直流电法层状方位各向异性介质中的电位分布、边界条件及视电阻率计算公式．以四极对称装置系统为例，对具有相同各向异性系数的4层模型采用核函数递推法作了理论数值模拟，得到了不同方向的电阻率测深曲线及其等值线形态．结果表明理论公式是正确的，测深曲线既反映了分层介质的电阻率差异，又反映了各层中电阻率的各向异性特征．      

High frequency S wave envelope synthesis using a multiple nonisotropic scattering model: application to aftershocks from the 2008 Wenchuan earthquake

Based on the formulation of a multiple non-isotropic scattering process, a characteristic source time is introduced to define the initial impulse width of energy density at the source. An analytical expression of the initial intensity spectral density of a seismic wave is incorporated into the integral equation of seismic wave energy density. And, a recursive formula of Green’s function is derived to obtain the higher order Green’s function, which is included to describe the stronger non-isotropic scattering process. Then, the effect of the scattering pattern on the energy density envelope is investigated by the modified scattering theory. Significant differences are found in the decay of the energy density envelopes with distances using different scattering patterns. The envelope synthesized by the forward dominated scattering pattern is larger than the results obtained by the isotropic and backward dominated scattering pattern. Different scattering patterns are also used to fit the observation data from the aftershocks of the 2008 Wenchuan earthquake. It is concluded that the envelopes synthesized by the forward scattering pattern can match the data better than the isotropic and backward dominated scattering cases, and a new interpretation of the coda wave is given. Finally, using the forward dominated scattering pattern, the envelope broadening of the observed data is reproduced.     

框架-阻尼框筒结构附加阻尼比简化计算方法

针对框架-阻尼框筒新型减震结构体系,考虑钢板阻尼墙屈服后强化和框架弹塑性,推导了结构附加阻尼比简化计算公式。以某框架-阻尼框筒结构为例,采用试验和有限元分析方法确定公式的参数取值,以此得到结构附加阻尼比的公式计算值,通过与结构动力弹塑性分析值进行对比,论证了简化计算公式的可用性。进一步基于参数分析给出了不同因素对公式计算值的影响。该简化计算方法为框架-阻尼框筒结构方案的快速设计提供了有效途径。     

General velocity formula of boundary layer above mobile sediment bed induced by asymmetric waves

Sediment movement in the wave boundary layer above a mobile sediment bed is complex.A velocity formula for the boundary layer is proposed for sheet flow induced by asymmetric waves above a mobile sediment bed.The formula consists of a free stream velocity and a defect function which contains a phase-lead,boundary layer thickness and mobile sediment bed.Phase-lag of sediment movement is considered in the formula for the mobile sediment bed.The formula needs six dependent variables about asymmetric wave and sediment characteristics.Asymmetry effects on parameters(orbital amplitude,roughness height,bed shear stress,and boundary layer thickness)are properly considered such that the formula can yield velocity differences among onshore,offshore,acceleration,and deceleration stages.The formula estimates the net boundary layer velocity resulting from the mobile sediment bed and asymmetric boundary layer thickness.In addition,a non-constant phase-lead also contributes to the net boundary layer velocity in asymmetric oscillatory sheet flow.Results of the formula are as good as that of a two-phase numerical model.Sheet flow transport induced by asymmetric waves,and the offshore net sediment transport rate with a large phase-lag under velocity-skewed waves,can be adequately estimated by the formula with a power sediment concentration function.     

LOCAL SCOUR DEPTH AROUND SPURS IN THE LOWER YELLOW RIVER

1 INTRODUCTION The Lower Yellow River is a suspended river channel bed much higher than the adjacent ground behind the dikes, well known all over the world. The floods flow east into the sea under the constraints of dikes. The safety of the dikes or embankments mainly relies on the protection provided by the river training works that have been constructed and maintained for several thousand years. There are three types of river training works and they are spurs, guide banks, and reve…     

Seismic safety of high concrete dams

Peak ground acceleration(PGA) estimation is an important task in earthquake engineering practice.One of the most well-known models is the Boore-Joyner-Fumal formula,which estimates the PGA using the moment magnitude,the site-to-fault distance and the site foundation properties.In the present study,the complexity for this formula and the homogeneity assumption for the prediction-error variance are investigated and an effi ciency-robustness balanced formula is proposed.For this purpose,a reduced-order Monte Carlo simulation algorithm for Bayesian model class selection is presented to obtain the most suitable predictive formula and prediction-error model for the seismic attenuation relationship.In this approach,each model class(a predictive formula with a prediction-error model) is evaluated according to its plausibility given the data.The one with the highest plausibility is robust since it possesses the optimal balance between the data fi tting capability and the sensitivity to noise.A database of strong ground motion records in the Tangshan region of China is obtained from the China Earthquake Data Center for the analysis.The optimal predictive formula is proposed based on this database.It is shown that the proposed formula with heterogeneous prediction-error variance is much simpler than the attenuation model suggested by Boore,Joyner and Fumal(1993).     

Deriving a practical form of IDF formula using transformed rainfall intensities

A general form of formula is presented for the rainfall Intensity–Duration–Frequency (IDF) relationship. This formula is derived from the nearly normal probability distribution function of transformed intensities. In order to transform the raw intensities, a correcting non‐constant spread technique, the Kruskal–Wallis statistic, and the Box–Cox transformation are adopted. These transformations enable to express a simpler model for the IDF formula that agrees well with traditional IDF relationships. Since the proposed method allows the estimation of any percentile value of intensities with a single equation, the intensity percentile at arbitrary duration can be generated easily. The validity of the formula derived by means of the proposed method is assessed using data from major weather stations in Korea. The results show that the percentile intensities produced using the proposed method are in good agreement with those of traditional frequency analysis. Copyright © 2013 John Wiley & Sons, Ltd.     

Vertical concentration profile of nonuniform sediment

A new formula for the concentration profile of nonuniform sediment is derived using the diffusion differential equation with equilibrium bed concentration.The interaction coefficient for nonuniform sediment particles as a function of both relative diameter and geometric standard deviation of nonuniform sediment is taken into account in the settling velocity.The diffusion coefficient is obtained from the logarithmic velocity profile.This new formula possesses several advantages as follows:(1)avoids theoretical defect of the Rouse formula,which states that the sediment concentration is infinite at the bottom and zero at the water surface.(2)suitable for an arbitrary fraction of nonuniform sediment,and(3)easy to apply with a simple form.The formula yields less difference among concentration profiles for various particles than that for uniform sediment and reveals the lawfor fractional concentration distribution of nonuniform sediment.The calculated concentration agrees well with the measured data from the Yangtze River.     

An aerodynamic formula to compute evaporation from open water surfaces

A formula to compute evaporation from open water surfaces is derived using wind profile characteristics. The mass transfer coefficient is given as a function of the shape and size of the water body in addition to the turbulent boundary layer parameters. The formula has been applied successfully to lakes and reservoirs of areas ranging from a fraction of an acre to thousands of acres. It is also given in a simplified form for application to cases where the wind profile measurements are not made.     

A probabilistic investigation of infiltration in the vadose zone: proposal for a new formula of infiltration rate

The infiltration rate in the unsaturated zone is analysed from a probabilistic point of view. It is shown that the empirical formulas of Horton and Kostiakov, without apparent physical basis, are explained in a probabilistic approach. Horton's formula reflects a Markovian process contrary to Kostiakov's formula. This approach made it possible to explain why Kostiakov's formula is more powerful than that of Horton. A new equation of infiltration is proposed. The three formulas were compared, for four types of soil, with the model of van Genuchten based on the Richards equation. Copyright © 2006 John Wiley & Sons, Ltd.     

Full 3-D numerical modeling of borehole electric image logging and the evaluation model of fracture

A full 3-D finite element method numerical modeling program is written based on the principle and technical specification of borehole electric image well logging tool. The response of well logging is computed in the formation media model with a single fracture. The effect of changing fracture aperture and resistivity ratio to the logging response is discussed. The identification ability for two parallel fractures is also present. A quantitative evaluation formula of fracture aperture from borehole electric image logging data is set up. A case study of the model well is done to verify the accuracy of the for-mula. The result indicates that the formula is more accurate than the foreign one.     

Relative Seismic Quiescence Before Large Aftershocks

--In terms of the modified Omori's formula and Akaike's Information Criterion, thirteen aftershock sequences have been analysed quantitatively, to check if any anomalous change of seismic activity took place before the occurrence of the largest aftershock. The analysis follows the method proposed by . The epicentre distribution of the main shocks covers the Aegean Sea (Greece) area and the surrounding regions. The earthquakes occurred from 1988 to 1998.¶Although difficulties existed for a quantitative treatment of seismic quiescence, the results of the temporal analysis of the aftershock activities seem to be promising. In eight of thirteen cases the observed temporal pattern is similar. The aftershock activity decreases significantly before the occurrence of the largest aftershock, with respect to the level expected according to the modified Omori's formula. Thereafter, the activity recovers to its previous level just before or immediately after the occurrence of the main aftershock. Moreover, an estimation of the reliability of the method is made, based on the results of this study as well as of previous similar investigations. From a total of 41 aftershock sequences which have been analysed following Matsu'ura's method, 29 of them (71%) show a remarkable relative quiescence before the occurrence of the main shock.¶Although it is difficult to construct a model which explains the precursory pattern, the real time monitoring and check for seismic quiescence could at least contribute to a 'qualitative' prediction of the largest aftershock which could at times be as large and disastrous as the main shock.     

Improved First-Order Approximation of Group Velocities in Weakly Anisotropic Media

A first-order approximation of the group velocity is derived for qP and qS waves in weakly anisotropic media. The formula gives an explicit expression of the group velocity in terms of elastic parameters and wave normal and is independent of any reference isotropic media. The approximated group velocity differs from the first order phase velocity in direction and in magnitude, the difference being of the first order in direction and the second order in magnitude. The accuracy of the approximate formula is tested on two examples of TI media. The formula well approximates the qS-waves group velocity surface even in the presence of triplications.     

Estimating inelastic response spectra from elastic spectra

Inelastic displacement response spectra are determined for a broad class of single-degree-of-freedom hysteretic structures. Based on these spectra, effective linear period and damping parameters are defined as a function of ductility. A simple empirical formula is derived which may be used to estimate the mid-period range inelastic response spectrum of a general hysteretic structure given the linear response spectrum of the excitation. The estimates obtained from this formula are compared with those obtained by the Newmark-Hall method, the substitute–structure method and the ATC–3 tentative procedure. It is found that the empirical formula not only gives good estimates of the average behaviour of the inelastic spectrum, but also reproduces some of the details of the spectrum.     

多点输入下大跨结构地震反应的频域精细传递矩阵法

地震输入问题一直是工程结构抗震研究所关注的焦点.在对大跨度结构进行抗震设计时,需要考虑多点地震输入的影响已成为国内外学术界和工程界的共识.本文根据傅氏变换,运用精细传递矩阵法推导出了大跨结构在多点地震输入下的频域精细传递矩阵,并在频域内对大跨度结构进行了动力分析.这种方法公式简单,能够快速、高精度地进行结构的地震反应分析.算例显示了本文方法的有效性.     

A new drain spacing formula

Abstract

A drain spacing formula is derived considering the variation in radial flux and the area above the drain level in the radial flow zone. The extent of the radial flow zone is ascertained by applying a mass balance and differentiability criterion of the water surface profile at the interface of radial and Dupuit-Forchheimer flow zones. The radial flow zone extends from the centre of the tile drain a distance of 2/π times the depth to impervious layer below the drain. For a normal ratio of recharge rate to hydraulic conductivity (R/K ≤ 0.0025), the water surface profile in the radial flow zone computed using Hooghoudt's formula is very different from the profile obtained by the new drain spacing formula; however, Hooghoudt's formula computes the maximum water table height which marginally differs from that found by the present method. For a ratio of high recharge rate to hydraulic conductivity (R/K = 0.1) and close drain spacing (L/D = 2), the difference in the maximum heights is 21%. Hooghoudt's formula overestimates the maximum water table position for L/D < 40. Unlike Hooghoudt's equivalent depth, the equivalent depth obtained using the present method is a function of the ratio of recharge rate to hydraulic conductivity.