地下管线三向地震动一致激励与非一致激励数值分析

基于某工业地下管道建立有限元模型,考虑管土相互作用和行波效应,综合分析在三向地震动一致激励与非一致激励下的管道动力响应结果,并且简单分析部分相关因素对管道动力响应结果的影响。结论如下:对于地下管道来说,非一致激励与一致激励作用下的位移响应曲线在峰值和形状方面存在较为明显的不同,具体表现为:非一致作用水平向位移远大于一致作用,而竖直向位移稍大于一致作用;对于应力响应,不管一致还是非一致激励,同一截面各处的应力响应有明显不同,并且总体来说非一致作用下较大;位于管道走向变化段、土层变化处和管道弯曲段的截面的位移响应峰值和应力响应一般会产生突变,说明这些因素对管道的动力响应具有较为明显的影响。     

Induction in laterally non-uniform conductors: Theory and numerical models

The study of electromagnetic induction in laterally non-uniform conductors is briefly reviewed. The two-dimensional perturbation problem is considered and the two polarization cases which arise from Maxwell's equations are discussed. Techniques for the solution of the equations for laterally non-uniform conductors are discussed with emphasis on numerical methods.     

2011年日本东北大地震特低频地磁信号的分形标度特征研究

本文选择2011年3月11日M9.0日本东北大地震震中附近三个地磁台站16个月(2010/01/01—2011/04/30)的特低频地磁观测资料,采用去倾扰动分析方法,得到了这三个台站地磁记录的非均匀标度特征随时间的变化,提出了一种能反映地磁三分量非均匀标度特征同步变化的指标,并据此探讨了特低频地磁信号分形标度特征变化与日本东北大地震之间的可能关联性,初步结果表明,这三个台站的特低频地磁信号分形标度特征指标在日本东北大地震前25~50天左右呈现出同步的异常增加,基于磁静日观测资料的随机合成地磁数据的统计检验结果可知前述异常并非随机异常,而是一种具有统计显著性的异常,可能反映了日本东北大地震对周边地磁信号内在的非线性系统特征产生了有统计意义的影响.     

Enhanced bed load sediment transport by unsteady flows in a degrading channel

Laboratory flume experiments were done to investigate bed load sediment transport by both steady and unsteady flows in a degrading channel. The bed, respectively composed of uniform sand, uniform gravel, or sand-gravel mixtures, always undergoes bulk degradation. It is found that both uniform and non-uniform bed load transport is enhanced greatly by unsteady flows as compared to their volume-equivalent steady flows. This enhancement effect is evaluated by means of an enhancement factor, which is shown to be larger with a coarser bed and lower discharges. Also, the fractional transport rates of gravel and sand in non-uniform sand-gravel mixtures are compared with their uniform counterparts under both steady and unsteady flows. The sand is found to be able to greatly promote the transport of gravel, whilst the gravel considerably hinders the transport of sand. Particularly, the promoting and hindering impacts are more pronounced at lower discharges and tend to be weakened by flow unsteadiness.     

Non-uniform sediment incipient velocity

Based on the mechanism of non-uniform sediment incipient motion, the dragging force and uplift force coefficient expressions are put forward for the non-uniform bed material exposure and close alignment state. The incipient velocity formula of the non-uniform sediment is then established. It is shown that the formula structure is reasonable, and fine particles of the non-uniform sediment are more difficult to set into motion than the same sized uniform particles, whereas coarser particles are easier to set into motion than their uniform equivalents. The validity of the formula is verified by field and experiment data.     

Response of an arch dam to non-uniform excitation generated by a seismic wave scattering model

Non-uniform ground motions are generated based on a single record available at a site and seismic wave scattering analysis. The Chino Hills 2008 earthquake records at the Pacoima Dam site are used to indicate the accuracy of the method. Dynamic analysis of the Pacoima dam-reservoir-foundation under uniform and non-uniform ground motions is carried out using the EACD-3D2008 software, and the results are compared to recorded responses at different locations on the dam. There is good agreement between computed and recorded displacements of the dam for non-uniform excitation. For uniform excitation, the displacements are underestimated in comparison with those obtained from recorded excitation. Significant intensification of stresses, especially near the foundation, and different patterns of stress distribution are observed for non-uniform excitation in comparison with uniform excitation. For uniform excitation maximum stresses occur in the crown cantilever near the crest, but for non-uniform excitation the maximum stresses occur along the sides and near the foundation.     

EFFECT OF COARSENING OF SURFACE BED MATERIAL ON NON-EQUILIBRIUM SEDIMENT TRANSPORT PROCESS DURING RIVER DEGRADATION

1INTRODUCTIONNon-equilibriumsedimenttransportina"at'Uralstreamwithnon-uniformbedmaterialisasubjectofilltensiveresearch.ThesedimentconcentrationmaybenotequaltothetransportcapacityofflowinanalluvialrivedItgraduallyapproachesequilibriumbydegradationoraggradationinalongdistance.usuallyover100kilometers.Thisproblemisofgreatimportancetothepredictionofthedistributionoferosionordepositionalongachannel,especiallyforlargeriversinChina.Scientistsstudiednonequilibriumsedimenttransportandthemainresul…     

Shaking table test of utility tunnel under non-uniform earthquake wave excitation

A series of shaking table tests were conducted on scaled utility tunnel models with and without construction joints under non-uniform input earthquake wave excitation. Details of experimental setup are first presented with particular focuses on: design and fabrication of double-axis laminar shear box with a rectangular hole opened on its side walls; design of two devices for measuring the slippage between the interface of test soil and the structure, and the relative deformation and rotation between joints of the structure model; and procedure for construction of input earthquake wave. The experiments were conducted in three phases. Phase 1 is free-field test. A 2-norm index is suggested to quantify the boundary effect and it is found that the designed laminar box does not impose significant boundary effect. Phases 2 and 3 are model tests in longitudinal and transversal directions, respectively. Test results are discussed in items of shear force–slip relationship at the soil–model structure interaction surface, movement and rotation of the construction joint, and effect of non-uniform earthquake input. The comparison shows that structural response under non-uniform earthquake excitation is larger than that under uniform excitation. The effect of spatial distribution of earthquake excitation should be considered in the seismic design of utility tunnel.     

Using finite element and Okada models to invert coseismic slip of the 2008 Mw 7.2 Yutian earthquake, China, from InSAR data

We invert InSAR from ascending and descending passes of ENVISAT/ASAR data to estimate the distribution of coseismic slip for the 2008 Mw 7.2 Yutian earthquake by separately using uniform finite element models (FEMs), non-uniform FEMs, and analytical elastic half-space (Okada) models. We use the Steepest Descent method and Laplacian smoothing to regularize and estimate the slip distribution for each of these different models. Fault surface ruptures interpolated from Quickbird optical images constrain the strike of the fault. The uniform FEM and Okada models assume a fault embedded in a homogeneous and elastic finite domain or half-space, respectively, while the non-uniform FEM describes a domain with spatially variable elastic properties according to the geological structure of Tibet plateau and Tarim. We find that the estimated slip distribution of the non-uniform FEM is similar to that of the uniform FEM but has a slightly larger moment (~2 %). However, more noticeable discrepancies occur between the slip distributions of the uniform FEM and Okada model, where the slip pattern estimated with the FEM looks more scattered and locates at a lower depth. We conclude that, for the case of the Yutian earthquake, the improvement brought by simulating the distribution of geologic material properties is insignificant. This suggests that the uniform geologic structure of Okada models may be sufficient for simulating large intra-continental earthquakes.     

Numerical and analogue modelling of induction effects in laterally non-uniform conductors

Modelling of electromagnetic fields in laterally non-uniform geoelectric structures is briefly reviewed. Attention is paid to progress achieved specifically during the last three years in numerical and analogue techniques of solving two- and three-dimensional problems. Some new important results are also mentioned.     

Modelling cyclic elastic behaviour of sands

The cyclic stress-strain behaviour of cohesionless soils is considered to be composed of elastic and inelastic components, and the cyclic elastic stress-strain relationships for one type of uniform sand are investigated by conducting isotropic and laterally constrained compression tests under low frequency cyclic loads. The effects of relative density, stress history, and stress state on the cyclic elastic behaviour are studied. Simple and practically applicable elastic stress-strain models are developed based on Hertzian contact theory. It was observed in both types of tests that elastic strains can accurately be modelled by a power function of stresses.     

山体遮挡对滇池风生流的影响初探

用二维风生流数值模型模拟滇池湖流运动。滇池在主导风向西南风作用下,假定湖面风场是均匀的,数值模拟的湖流流态与实测湖流结果相差很大。而考虑山体遮挡影响,根据实测湖流期间现有的风情资料,在湖面上构造一非均匀风场,数值模拟结果与实测值基本一致。山体遮挡对滇池风生流的影响是不容忽视的。建议进一步进行湖流和湖面风向、风速监测,并建立过山气流数学模型,深入研究山体遮挡对湖泊风生流的影响。     

A sample of controlled experiments in diff radion tomography1

Various applications of a new geophysical reconstruction method, generalized acoustical diffraction tomography (GADT), which is based on transmission data as input are considered. Conventional diffraction tomography methods normally require linearization with respect to a uniform reference medium and regular sampling along a straight line. Thus, these methods will not work well when the background is strongly non-uniform and/or the acquisition geometry is arbitrary. However, GADT can, in principle, handle both irregularly spaced data, curved acquisition lines, and non-uniform background models. A number of controlled model tank and field experiments, where the model and the test object(s) are known a priori, have been carried out. After acquiring the tomographic data in each experiment, these are used to compute a reconstruction of the model, which can then be compared with the actual, known model. The method's ability to yield high-quality images of the different targets is demonstrated.     

超长混凝土框架梁弹性内力行波效应研究

鉴于超长混凝土框架结构行波效应规律尚缺乏系统研究,借助建筑结构通用有限元软件MI-DAS/GEN,采用一致与非一致两种激励方式对165～1100m/s地震波速作用下长度为33～264m的3层钢筋混凝土框架结构进行弹性时程分析。通过对比一致与非一致激励下的梁弯矩、剪力和轴力,详细讨论了其随波速、长度以及不同部位的变化规律。研究表明:波速越慢,模型越长,行波效应越明显。框架梁弯矩的增大作用主要集中在264m模型第1层两端,减小作用表现在各模型第1层中部和第2、3层的梁上;剪力与弯矩有相似的变化规律;梁中存在拉力与压力,中部各梁的拉力和压力较两端的梁要大。     

非一致地震激励地下综合管廊振动台模型试验研究(Ⅱ)——试验结果

本文以纵向非一致地震激励试验为主,对地下综合管廊振动台模型试验结果进行了分析,包括模型场地的加速度响应、结构的应变响应、结构和场地加速度响应关系、结构应变响应的原因、结构接头位移响应、结构在纵向非一致激励和一致激励下应变响应的区别。结果分析表明,本试验模型场地及模型结构设计合理,从中得到了关于地下综合管廊非一致地震激励作用下响应规律的一些新内容。     

Tomographic determination of the upper crustal structure in the Jiashi strong earthquake swarm region

IntroductionBetween January and April of 1997, 7 earthquakes with M(6.0 occurred successively in Jiashi, Xinjiang. The continual occurrence of strong earthquakes within such a small area and in such a short period of time is exceptional for intraplate earthquakes. The Jiashi earthquake swarm took place on the northeast side of the Pamirs, where the Tarim basin, South Tianshan and West Karakoram meet (HU, et al, 1989). This is also a place where a number of active faults develop, so it is…     

High time resolution monitoring of tropospheric temperature with a Radio Acoustic Sounding System (RASS)

We have observed the time-height variation of the temperature field in the upper troposphere using a Radio Acoustic Sounding System (RASS) which consists of the MU radar and a high-power acoustic transmitter. The fast beam steerability of the MU radar has made it possible to measure temperature profiles in a fairly wide height range in the upper troposphere (5–11 km), even under intense wind conditions. Observations were continued for about 32 hr on 24–26 December, 1986 with a time-height resolution of 30 min and 150 m. During the observation period, the tropospheric jet was so intense that the acoustic wavefronts were severely distorted. Using wind velocity profiles observed by the MU radar we have numerically estimated the propagation of acoustic wavefronts, and further determined favorable pointing directions for the MU radar to receive significant backscattering from refractive index fluctuations produced by the acoustic waves. Conventional radiosonde soundings were carried out every 6 hr, which showed a temperature decrease of 4 K/day in the upper troposphere during the observation period. Temperature profiles taken by RASS agree well with the radiosonde results.     

Oblique surface wave propagation over a small undulation on the bottom of an ocean

The problem of oblique water wave diffraction by small undulation of the bottom of a laterally unbounded ocean is considered using linear water wave theory. A perturbation analysis is employed to obtain the velocity potential, the reflection and the transmission coefficients up to the first order in terms of integrals involving the shape functions c(x) representing the bottom undulation. Finite cosine transform is used to find the first order potential, and this potential is utilised in obtaining the first order reflection and transmission coefficients. Some particular forms of the shape function representing an exponentially damped undulation, a single hump and a patch of sinusoidal ripples are considered and the integrals for the reflection and transmission coefficients are evaluated. For the exponentially damped undulation, it is observed that the reflection ceases much before transmission while for the single hump, reflection and transmission go hand in hand up to a certain value of the wavenumber, after which they vanish. For the patch of sinusoidal ripples having the same wavenumber, the reflection coefficient up to the first order is found to be an oscillatory function in the quotient of twice the component of the wavenumber along x-axis and the ripple wavenumber. When this quotient becomes one, the theory predicts a resonant interaction between the bed and free surface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of the incident wave energy occurs if this number is large. Also, when a patch of ripples having different wavenumbers is considered the same result follows. Known results for the normal incidence are recovered as special cases for the patch of sinusoidal ripples. The theoretical observations are shown computationally.     

Ambiguous reflection coefficients for anelastic media

The calculation of reflection and transmission coefficients of plane waves at a plane interface between two homogeneous anelastic media may become ambiguous because it is not always obvious how to determine the sign of the vertical component of the slowness vector of the scattered waves. For elastic media, the sign is determined by applying so-called radiation condition when the slowness vector is complex-valued, but it has long been known that this approach does not work satisfactorily for anelastic media. Other approaches have been suggested, e.g., by requiring that the reflection and transmission coefficients should vary continuously with increasing incident angles, or by relating the sign to the direction of the energy flux. In the present paper, it is shown that these approaches may give different results, and that the results can be inconsistent with the elastic case even for weak attenuation. Instead, it is demonstrated that the ambiguity in the reflection coefficient can be resolved by expressing the seismic response of a point source over an interface as a superposition of plane waves and their reflection coefficients, and solving the resulting integral by the saddle point approximation. Although the saddle point itself (point of stationary phase) does not provide new insight, the ambiguity is removed by considering the steepest descent path through the point. Ray synthetic seismograms computed by this method compare well with synthetics computed by the reflectivity method, which does not suffer from the above-mentioned ambiguity since the integration path is taken along the real axis. This paper concentrates on the isotropic case, but it is discussed how the result may be extended to layered transversely isotropic media. The suggested approach, derived for a point source and plane layers, does not directly apply to 2-D or 3-D laterally inhomogeneous media, or to media of general anisotropy. A generalization of the result found is that the sign of the vertical slowness components should be chosen according to the energy flux direction for subcritical incidence and according to the radiation condition for supercritical incidence, even if this creates a discontinuity in the coefficients at the critical incidence angle. Such a discontinuity is sometimes necessary to get results which are consistent with the elastic case. It is discussed how the generalized result can be obtained by applying certain continuity criteria for the sub-and supercritical angle intervals, but the validity of this approach for general models remains to be proved.     

Shear waves at a corrugated interface between anisotropic elastic and visco-elastic solid half-spaces

Rayleigh’s method of approximation is employed to find out the reflection and transmission coefficients due to an incident plane SH wave at a corrugated interface between a laterally and vertically inhomogeneous anisotropic elastic solid half-space and a laterally and vertically inhomogeneous isotropic visco-elastic solid half-space. The lateral and vertical inhomogeneities are described by the exponential variations of elastic parameters. The formulae of reflection and transmission coefficients are derived in closed form for the first-order approximation of the corrugation. The effects of the corrugation of the interface, the inhomogeneity, the anisotropy, the visco-elasticity and the frequency of the incident wave on these coefficients are studied analytically and numerically for a specific model containing a periodic interface. The results of earlier workers have been reduced as particular cases from the present formulation.