径向分层TI孔隙介质井孔中激发的模式波的数值研究

横向各向同性介质是地层中普遍存在的一种各向异性介质.本文对径向分层TI孔隙介质包围井孔中激发的斯通利波和弯曲波的传播特性进行了理论计算,发现模式波在低频时更多的是反应原状地层的信息,而随着频率的增加侵入带参数逐渐起控制作用;Biot理论描述的地层衰减比速度更容易受井壁附近地层参数的影响.利用灵敏度曲线定量研究了不同频率下地层各个参数对相速度和衰减系数的贡献大小,主要结果显示模式波的衰减受水平渗透率影响明显,而垂直渗透率的变化对模式波几乎无影响;斯通利波对水平向传播的横波速度比弯曲波的灵敏度高.从单极子和偶极子声源在井孔中激发的全波波形也可发现,声波测井仪器较宽的声源频带和合适的源距设置有利于对不同径向深度上的地层声学参数进行成像.     

轴对称二维任意非均匀介质中偏轴点源激励的位场数值解

对具有任意多平面分层，且每一平面分层中又具有任意多径向柱面分层的复杂非均匀介质中偏轴点源激励的位场分布进行了高效数值分析．对于这种三维非均匀位场分布，应用数值本征模展开描述其径向的基模和高阶模，用集合反射矩阵和集合透射矩阵匹配各平面界面的边界条件，各平面分层中的位场则用递推算法获得．偏轴点源位场的求解，为电法测井中高分辨薄五层探测提供了理论分析基础．     

Dynamics of rigid strip foundations embedded in orthotropic elastic soils

This study is concerned with the dynamic response of an arbitrary shaped rigid strip foundation embedded in an orthotropic elastic soil. The foundation is subjected to time-harmonic vertical, horizontal and moment loadings. The boundary-value problem related to an embedded foundation is analysed by using the indirect boundary integral equation method. The kernel functions of the integral equations are displacement and traction Green's functions of an anisotropic elastic half plane. Exact analytical solutions are used for the Green's functions. The boundary integral equation is solved by using numerical techniques. Selected numerical results are presented for the impedances of rectangular and semi-circular rigid strip foundations embedded in four types of anisotropic soils. A discussion on the influence of soil anisotropy and frequency of excitation on the impedances is presented. The versatility of the analysis is demonstrated by considering the through soil interaction between two semi-circular strip foundations.     

Assessing radial transmissivity variation in heterogeneous aquifers using analytical techniques

Pumping tests are one of the most commonly used in situ testing techniques for assessing aquifer hydraulic properties. Numerous researches have been conducted to predict the effects of aquifer heterogeneity on the groundwater levels during pumping tests. The objectives of the present work were as follows: (1) to predict drawdown conditions and to estimate aquifer properties during pumping tests undertaken in radially symmetric heterogeneous aquifers, and (2) to identify a method for assessing the transmissivity field along the radial coordinate in radially symmetric and fully heterogeneous transmissivity fields. The first objective was achieved by expanding an existing analytical drawdown formulation that was valid for a radially symmetric confined aquifer with two concentric zones around the pumping well to an N concentric zone confined aquifer having a constant transmissivity value within each zone. The formulation was evaluated for aquifers with three and four concentric zones to assess the effects of the transmissivity field on the drawdown conditions. The specific conditions under which aquifer properties could be identified using traditional methods of analysis were also evaluated. The second objective was achieved by implementing the inverse solution algorithm (ISA), which was developed for petroleum reservoirs to groundwater aquifer settings. The results showed that the drawdown values are influenced by a volumetric integral of a weighting function and the transmissivity field within the cone of depression. The weighting function migrates in tandem with the expanding cone of depression. The ability of the ISA to predict radially symmetric and log‐normally distributed transmissivity fields was assessed against analytical and numerical benchmarks. The results of this investigation indicated that the ISA method is a viable technique for evaluating the radial transmissivity variations of heterogeneous aquifer settings. Copyright © 2013 John Wiley & Sons, Ltd.     

Effectiveness of vertical drains in mitigation of liquefaction

One method of mitigating the damaging effects of earthquake-induced liquefaction is to provide rapid dissipation of excess pore pressures by the use of vertical drains through the liquefiable material. Drain systems are currently designed using a chart-based approach. Field experience suggests that the performance of these installations cannot yet be accurately predicted.

In this paper, high quality centrifuge testing is used to help clarify drain behaviour. It will be established, supported by centrifuge test data, that the pore water from a radially expanding zone of soil contributing to drainage through the drains is developed. Naturally, the geometry of this expanding zone changes with time. It will be shown that fluid from deeper strata is drained first, reducing the effectiveness of the drain for near-surface soil layers. It is concluded that these zones are useful in analysing more complicated drain system geometries.     

青藏高原东北缘上地幔多尺度层析成像

利用布设在青藏高原东北缘地区的甘肃宽频带地震台阵记录到的远震P波走时数据,采用小波域参数化和基于L1范数的稀疏约束反演算法的多尺度层析成像方法,得到了该地区400km深度范围内上地幔的P波速度结构.本文采用的多尺度层析成像方法可以自适应数据非均匀采样的情况,有效降低谱泄漏效应和反问题的多解性,明显提高解的分辨率和可靠性.层析成像结果表明青藏高原东北缘上地幔整体上显示为低速特征,扬子地块上地幔则显示为高速特征,两者之间上地幔存在清晰的块体边界带,该边界带位于东经104°—105°之间并且随深度的增加逐渐东移.该特征暗示了青藏高原上地幔物质向东扩张的机制,但在西秦岭上地幔顶部不存在物质运移的通道.青藏高原东北缘内部也具有明显的分区特征,松潘—甘孜地块上地幔P波速度整体呈低速特征,而柴达木地块的上地幔顶部具有相对高速特征,而在上地幔200km以下这两个地块间的差别逐渐减小.1654年天水地震和1879年武都地震都发生在扬子地块与青藏高原的碰撞交汇区,其震中下方上地幔显示为高低速转换结构.     

Static inelastic analysis of RC shear walls

A macro-model of a reinforced concrete (RC) shear wall is developed for static inelastic analysis. The model is composed of RC column elements and RC membrane elements. The column elements are used to model the boundary zone and the membrane elements are used to model the wall panel. Various types of constitutive relationships of concrete could be adopted for the two kinds of elements. To perform analysis, the wall is divided into layers along its height. Two adjacent layers are connected with a rigid beam. There are only three unknown displacement components for each layer. A method called single degree of freedom compensation is adopted to solve the peak value of the capacity curve. The post-peak stage analysis is performed using a forced iteration approach. The macro-model developed in the study and the complete process analysis methodology are verified by the experimental and static inelastic analytical results of four RC shear wall specimens. Supported by: National Natural Science Foundation of China, Grant number 59895410     

Theory and application of equivalent transformation relationships between plane wave and spherical wave

Based on the governing equations and the equivalent models,we propose an equivalent transformation relationships between a plane wave in a one-dimensional medium and a spherical wave in globular geometry with radially inhomogeneous properties.These equivalent relationships can help us to obtain the analytical solutions of the elastodynamic issues in an inhomogeneous medium.The physical essence of the presented equivalent transformations is the equivalent relationships between the geometry and the material properties.It indicates that the spherical wave problem in globular geometry can be transformed into the plane wave problem in the bar with variable property fields,and its inverse transformation is valid as well.Four different examples of wave motion problems in the inhomogeneous media are solved based on the presented equivalent relationships.We obtain two basic analytical solution forms in Examples Ⅰ and Ⅱ,investigate the reflection behavior of inhomogeneous half-space in Example Ⅲ,and exhibit a special inhomogeneity in Example Ⅳ,which can keep the traveling spherical wave in constant amplitude.This study implies that our idea makes solving the associated problem easier.     

Effect of non-linear soil–structure interaction on seismic response of tall slender structures

Some structures may be very massive and may have to be located on relatively soft soil. In such cases, the soil adjacent to the structure behaves in a non-linear fashion and affects the response of the structure to the dynamic loading. An approximate hybrid approach to analyse soil–structure systems accounting for soil non-linearities has been developed in this paper. The approach combines the consistent infinitesimal finite-element cell method (CIFECM) and the finite-element method (FEM). The CIFECM is employed to model the non-linear (near-field) zone of the soil supporting the structure as a series of bounded media. The material properties of the bounded media are selected so that they are compatible with the average effective strains over the whole bounded medium during the excitation. The linear zone of soil away from the foundation, the far-field, is modelled as an unbounded medium using the CIFECM for unbounded media. The structure itself is represented by the FEM. The proposed method is used to model the dynamic response of a one-mass structure and a TV-tower supported on a homogenous stratum and excited by an earthquake. It was found that the secondary soil non-linearity might increase or decrease the base forces of tall slender structures depending on the type of structure, frequency content of the input motion and the dynamic properties of the near-field soil.     

Scattering of shear waves by an elliptical cavity in a radially inhomogeneous isotropic medium

Complex function and general conformal mapping methods are used to investigate the scattering of elastic shear waves by an elliptical cylindrical cavity in a radially inhomogeneous medium. The conformal mappings are introduced to solve scattering by an arbitrary cavity for the Helmholtz equation with variable coefficient through the transformed standard Helmholtz equation with a circular cavity. The medium density depends on the distance from the origin with a power-law variation and the shear elastic modulus is constant. The complex-value displacements and stresses of the inhomogeneous medium are explicitly obtained and the distributions of the dynamic stress for the case of an elliptical cavity are discussed. The accuracy of the present approach is verified by comparing the present solution results with the available published data. Numerical results demonstrate that the wave number, inhomogeneous parameters and different values of aspect ratio have significant influence on the dynamic stress concentration factors around the elliptical cavity.     

Spectral methods for bulk and guided waves in weakly inhomogeneous media

Plane wave spectra localized around interference maxima provide the foundation for charting the propagation of source excited wavefields through weakly inhomogeneous bulk and layered environments. Attention is focused here on recent developments in spectral analysis, uniformization, restoration, and inversion, centered around the ‘bare bones’ spectral objects represented by ray fields in a bulk medium, by adiabatic modes in a layered medium, and by beams.     

THEORETICAL RESISTIVITY SOUNDING RESULTS OVER A TRANSITION LAYER MODEL *

An earth model with a transition layer (anisotropic inhomogeneous) is considered. The inhomogeneity in σ_v (vertical conductivity) of the transition layer is represented by a power law variation. Expressions for potential distribution in the upper layer, transition layer and bottom layer are obtained by solving appropriate differential equation for each layer. By utilizing the boundary conditions, expressions of apparent resistivity for Wenner and Schlumberger configurations are derived. Numerical analysis is performed for linear and quadratic variation of σ_v. The results are presented in the form of theoretical apparent resistivity curves for both configurations. Negative apparent resistivities are the interesting feature of this analysis.     

Meridional circulation in the sun and stars

Mean-field hydrodynamics advanced to clear explanations for the origin and properties of the global meridional flow in stellar convection zones. Qualitative arguments and analysis of basic equations both show that the meridional circulation is driven by non-conservative centrifugal and buoyancy forces and results from a slight disbalance between these two drivers. The deviations from the thermal wind balance are relatively large near the boundaries of convection zones. Accordingly, the meridional flow attains its largest velocities in the boundary layers and decreases inside the convection zone. This picture, however, is neither supported nor dismissed by the conflicting results of recent helioseismic soundings or 3D numerical experiments. The relevant physics of the differential temperature and its possible relation to the solar oblateness are briefly discussed.     

Vertical and horizontal vibrations of a rigid disc on a multilayered transversely isotropic half-space

A half-space containing horizontally multilayered regions of different transversely isotropic elastic materials as well as a homogeneous half-space as the lowest layer is considered such that the axes of material symmetries of different layers and the lowest half-space to be as depth-wise. A rigid circular disc rested on the free surface of the whole half-space is considered to be under a forced either vertical or horizontal vibration of constant amplitudes. Because of the involved integral transforms, the mixed boundary value problems due to mixed condition at the surface of the half-space are changed to some dual integral equations, which are reduced to Fredholm integral equations of second kind. With the help of contour integration, the governing Fredholm integral equations are numerically solved. Some numerical evaluations are given for different combinations of transversely isotropic layers to show the effect of degree of anisotropy of different layers on the response of the inhomogeneous half-space.     

The structure of the Vladikavkaz Fault Zone based on the study utilizing a complex of geological-geophysical methods

The Vladikavkaz Fault Zone is a regional boundary between the Alpine folded-block system of the Greater Caucasus and Pre-Caucasus foredeep (Ossetian Depression). A complex of geological and geomorphological techniques was used for its study. The data of recent activity of discontinuous and folded dislocations were obtained as a result. Additional application of the microseismic sounding method allowed us to trace the behavior at depth and along the direction of the Vladikavkaz Fault Zone’s branches known on the Earth’s surface from geological and geophysical surveys. An idea of the three-dimensional structure of this system of tectonic faults was made up. It was shown that the structure of the upper section of the fault zone has a form of listric foldthrust, and in the deep layers of the Earth’s crust it is represented by a well-formed contrast zone with subvertical plunging. Expressive deformations of young sediments and a clear manifestation of the fault zone at a depth make it possible to consider the Vladikavkaz Fault as a major seismogenerating structure.     

Effects of soil–structure interaction on the dynamic properties and seismic response of piled structures

This paper presents a simple and stable procedure for the estimation of periods and dampings of piled shear buildings taking soil–structure interaction into account. A substructuring methodology that includes the three-dimensional character of the foundations is used. The structure is analyzed as founded on an elastic homogeneous half-space and excited by vertically incident S waves. The strategies proposed in the literature to estimate the period and damping are revised, and a modified strategy is proposed including crossed impedances and all damping terms. Ready-to-use graphs are presented for the estimation of flexible-base period and damping in terms of their fixed-base values and the system configuration. Maximum shear forces together with base displacement and rocking peak response are also provided. It is shown that cross-coupled impedances and kinematic interaction factors need to be taken into account to obtain accurate results for piled buildings.     

Influence of Outcrop Scale Fractures on the Effective Stiffness of Fault Damage Zone Rocks

We combine detailed mapping and microstructural analyses of small fault zones in granodiorite with numerical mechanical models to estimate the effect of mesoscopic (outcrop-scale) damage zone fractures on the effective stiffness of the fault zone rocks. The Bear Creek fault zones were active at depths between 4 and 15 km and localize mesoscopic off-fault damage into tabular zones between two subparallel boundary faults, producing a fracture-induced material contrast across the boundary faults with softer rocks between the boundary faults and intact granodiorite outside the boundary faults. Using digitized fault zone fracture maps as the modeled fault geometries, we conduct nonlinear uniaxial compression tests using a novel finite-element method code as the experimental “laboratory” apparatus. Map measurements show that the fault zones have high nondimensional facture densities (>1), and damage zone fractures anastamose and intersect, making existing analytical effective medium models inadequate for estimation of the effective elastic properties. Numerical experiments show that the damage zone is strongly anisotropic and the bulk response of the fault zone is strain-weakening. Normal strains as small as 2% can induce a reduction of the overall stiffness of up to 75%. Fracture-induced effective stiffness changes are large enough to locally be greater than intact modulus changes across the fault due to juxtaposition of rocks of different lithologies; therefore mesoscopic fracturing is as important as rock type when considering material or bimaterial effects on earthquake mechanics. These results have important implications for earthquake rupture mechanics models, because mesoscopic damage zone fractures can cause a material contrast across the faults as large as any lithology-based material contrast at seismogenic depths, and the effective moduli can be highly variable during a single rupture event.     

High-resolution seismic stratigraphy of the Yamato Basin, Japan Sea and its geological application

Abstract High-resolution seismic stratigraphy of the Yamato Basin, Japan Sea, was successfully established using core-log-seismic data integration. The construction of synthetic seismograms by the combination of physical properties and well-log data from the Ocean Drilling Program (ODP) Site 797 was the key to accomplishing the high-resolution seismic stratigraphy. To achieve resolution comparable with well-log data and core lithology, single channel seismic reflection data taken from ODP underway geophysics were reprocessed, and then carefully compared with synthetic seismogram, core and well log profiles to identify seismic units. Ten seismic stratigraphic units were identified at the site, and seismic stratigraphic interpretation was successfully extended from the site to the nearby area along the Yamato Basin margin. The opal-A/opal-CT (biogenic silica/metastable diagenetic silica) boundary has different appearances at places from strong to weak, and mostly discontinuous. One of the significant results achieved from this study is clear distinction of the opal-A/CT boundary from a very strong reflector, which appears at 22 m below the opal-A/CT boundary. Through well-log and physical properties characterization of the different units, resistivity was found to be the best indicator of diatom content and with gamma-ray it also is an indicator of chert layers in the opal-CT zone. Velocity is not greatly effected by diatom ooze in the opal-A zone, however, it shows strong peaks and has an indirect relationship with gamma-ray in the opal-CT zone. Finally, successful correlation of Gamma-ray Attenuation Porosity Evaluator density and resistivity peaks with strong seismic reflectors from upper and lower stratified layers may provide new information on the late Neogene paleoceanography of the Japan Sea in high-resolution scale.     

点源两层大地三维地电模型视电阻率边界元解

本文用两层均匀大地人工点源作用下的空间电位分布作为基本解，导出了下伏基岩中陷含三维地质体情况下，关于电位的边界积分方程组，用边界单元法解积分方程求得地表电位，从而求得了一定装置下的视电阻率。本文方法的独特优点是积分只须在异常体边界面上进行了测线可以沿任意方向。用联剖法计算了均匀半空间下含球形异常体的ρｓ曲线，所得结果与其他学者用解析方法所求得的结果一致。用中间梯度法计算下层基岩中含有一长方体，测线     

Love wave dispersion: A comparison of results for a semi-infinite medium with inhomogeneous layers and for its approximation by homogeneous layers

Love wave dispersion in various semi-infinite media consisting of inhomogeneous layers is discussed. The phase and group velocities are computed when shear wave velocity and density in each inhomogeneous layer are varying exponentially with depth. At the beginning one or two inhomogeneous layers over a homogeneous semi-infinite medium are considered. The dispersion results for these structures are compared with those for their approximations with homogeneous layers. Comparisons show that differences of phase and group velocities for the original models from those for their approximated models (i) increase with the increase of wave number and (ii) are larger for group velocity than for phase velocity. The difference is approximately proportional to the rate of change of parameters in the layers. Finally, dispersion curves are obtained for model IP3MC, which consists of many inhomogeneous and homogeneous layers over a homogeneous semi-infinite medium. The results are compared with the observed group velocity data across the Indian Peninsula.