Normally incident wave propagation across a joint set with the virtual wave source method

This paper presents an analytical study on normally incident wave propagation across a joint set. Combined with the displacement discontinuity method (DDM), the virtual wave source method (VWSM) is adopted to take account of multiple wave reflections among joints. After theoretical formulation, a comparison between results from VWSM and those from the method of characteristics (MC) is carried out. These results agree well with each other, and the applicability of VWSM into wave propagation study is verified. Subsequently, an interesting phenomenon, i.e., the magnitude of transmission coefficient increases with joint number on some circumstances, is pointed out and studied. By analyzing differently arriving transmitted waves and their superposition, such phenomenon is explained.     

Seismic response of an embedded pile in a transversely isotropic half-space under incident P-wave excitations

A rigorous mathematical formulation is presented for the analysis of a thin cylindrical shell embedded in a transversely isotropic half-space under vertically incident P-wave excitation. By virtue of a set of ring-loads Green's functions for the shell and a group of dynamic fundamental solutions for the half-space under arbitrary interfacial dynamic loads, the problem is shown to be reducible to a pair of Fredholm integral equations. By utilizing an adaptive-gradient family capable of capturing regular-to-singular solution transitions smoothly, an accurate numerical procedure is developed. To assess the effect of material anisotropy on the dynamic load-transfer process, a set of comprehensive numerical results presented for various material and geometrical conditions. The accuracy of the proposed numerical scheme is confirmed by its comparison with a benchmark solution for the corresponding isotropic problem.     

Influence of water saturation on horizontal and vertical motion at a porous soil interface induced by incident SV wave

Recent field observations have indicated that water saturation of soils may strongly affect the vertical ground motion. A study is therefore carried out to investigate the effect of saturation on horizontal and vertical motion at an interface of porous soils with potential contributions directed to site evaluation based on field observations of both the horizontal and vertical motion. The problem described in this paper corresponds to an SV wave incident at the interface between the overlying soil and the underlying rock formation. The soils are modeled as partially water-saturated porous material with a small amount of air inclusions, while the rock are approximately regarded as ordinary one-phase solid. Theoretical formulation is developed for the computation of amplitudes of horizontal and vertical interface motion, which are expressed as functions of the degree of saturation, the angle of incidence as well as the frequency. Numerical results are given for a typical sand to illustrate the influence of saturation on the interface motion in two directions and their ratios. The present study demonstrates that the effect of water saturation may be substantial on both the horizontal and vertical motion as well as on their ratios, implying the importance of such effects in the interpretation of field observations.     

Influence of water saturation on horizontal and vertical motion at a porous soil interface induced by incident P wave

This paper describes a study on the influence of water saturation on horizontal and vertical motion at the interface between porous soil and rock formation due to an incident P wave. The soil is modeled as a partially water-saturated porous material with a small amount of air inclusions, while the underlying rock is approximately regarded as an ordinary one-phase solid. Theoretical formulation is developed for the computation of motion amplitudes in both horizontal and vertical components, which are considered as functions of the degree of saturation, the angle of incidence as well as the frequency. Numerical results are provided to illustrate the effect of saturation on displacement amplitudes in both components and their ratios. The results show that even a slight decrease of complete saturation may lead to a substantial change in the motion amplitudes in both horizontal and vertical components. In general, partial saturation may cause lower horizontal-to-vertical ratios over the entire range of incident angles except the normal incidence, implying the potential importance of saturation effects in the interpretation of field observations.     

Forced wave motion in an unbouded space surrounding a lined,spherical cavity

Summary This investigation considers the point-symmetric response of an unbounded, elastic medium. A spherical cavity in the unbounded medium is lined by an elastic, spherical shell attached to the surrounding medium. The shell liner is subjected to specified time-varying surface stress acting upon its internal surface. A detailed evaluation of the shell liner shielding characteristics is presented.     

P波入射时含圆柱形空腔板弹性波散射与动应力集中

采用波函数展开法给出了板的横截面处孔洞在平面P波入射下动应力集中问题的解。本文使用大圆弧假定法来满足板中的自由表面的边界条件。数值结果表明,上表面和孔洞间的距离与下表面和孔洞间的距离都对结果有影响。如果孔洞半径相对板的厚度很小,解近似等于孔洞全空间解。所给出的解是半解析解。因此能被用来验证诸如BEM,FEM,FDM等数值法。     

倾斜地层中的井孔声场研究

研究声波在倾斜充液井孔中的传播对于声波测井数据处理和解释具有重要意义.应用三维交错网格有限差分方法模拟了处于倾斜各向同性分层地层中的井孔声场.首先,针对均匀地层中单极子声源在裸眼井中激发的声场,将有限差分的结果和实轴积分法的结果进行对比验证.然后,采用单极子和偶极子两种声源,针对地层分界面和井轴间的不同倾角,计算了相应的声场分布和井轴上的接收波形.数值计算的结果表明,当声源处于倾斜分界面以下,即处于快速（下方）地层,接收器处于倾斜分界面以上（慢速）地层时,随着地层倾斜角度的加大,测得的慢度值从接近上方慢速地层值逐渐减小直至接近下方快速地层的值.任何源距情况下测得的首波慢度均小于上方地层实际的纵波慢度.并且,慢度与源距的关系曲线随源距的加大逐渐平缓.用偶极子声源激发得到的横波慢度与纵波结果相同,并表现得比纵波对倾角的改变更敏感.上述结论在本文中用声场快照和利用合成接收波列的慢度计算得以清楚显示,并且用射线声学理论验证.     

Harmonic wave diffraction by two circular cavities in a poroelastic formation

This study investigates the dynamic interaction of time harmonic plane waves with a pair of parallel circular cylindrical cavities of infinite length buried in a boundless porous elastic fluid-saturated medium. The novel features of Biot dynamic theory of poroelasticity along with the appropriate wave field expansions, the pertinent boundary conditions, and the translational addition theorems for cylindrical wave functions are employed to develop a closed-form solution in the form of infinite series. The analytical results are illustrated with numerical examples in which two empty cavities are insonified by a fast compressional or a shear wave at end-on incidence. The basic dynamic field quantities such as the hoop stress amplitude and the radial displacement of the elastic frame are evaluated and discussed for representative values of the parameters characterizing the system. The effects of the proximity of the two cavities, the incident wave frequency and type are examined. Particular attention has been focused on multiple scattering interactions in addition to the slow wave coupling effects which is known to be the primary distinction of the scattering phenomenon in poroelasticity from the classical elastic case. Limiting case involving two empty cylindrical cavities in an elastic solid is considered and excellent agreement with a well-known solution is established.     

On the electron dynamics in the field of a whistler wave propagating along a magnetic field in a weakly inhomogeneous plasma

The relativistic motion of electrons in the field of a finite amplitude circularly polarized whistler propagating along a constant magnetic field in a plasma with longitudinal inhomogeneity is considered. It is shown that the equations of the particle motion exhibit a constant of the motion despite the spatial dependence of the wave parameters. The existence of the constant allows one to reduce the equations of motion to a canonical form describing one-dimensional oscillations of a particle with the Hamiltonian slowly varying in the process of the oscillations.     

Heat transfer between a thermochemical plume channel and the surrounding mantle in the presence of horizontal mantle flow

Heat and mass transfer processes in the conduit of a thermochemical plume located beneath an oceanic plate far from a mid-ocean ridge (MOR) proceed under conditions of horizontal convective flows penetrating the plume conduit. In the region of a mantle flow approaching the plume conduit (in the frontal part of the conduit), the mantle material heats and melts. The melt moves through the plume conduit at the average velocity of flow v and is crystallized on the opposite side of the conduit (in the frontal part of the conduit). The heat and the chemical dope transferred by the conduit to the mantle flow are carried away by crystallized mantle material at the velocity v.The local coefficients of heat transfer at the boundary of the plume conduit are theoretically determined and the balance of heat fluxes through the side of the plume conduit per linear meter of the conduit height. The total heat generation rate, transmitted by the Hawaiian plume into the upper and lower mantle, is evaluated. With the use of regular patterns of heat transfer in the lower mantle, which is modeled on the horizontal layer, heated from below and cooled from above, the diameter of the plume source, the kinematic viscosity of the melt in the plume conduit, and the velocity of horizontal lower-mantle flows are evaluated and the dependences of the temperature drop, viscosity and Rayleigh number for the lower mantle on the diameter of the plume source are presented.     

ELF wave generation in the ionosphere using pulse modulated HF heating: initial tests of a technique for increasing ELF wave generation efficiency

This paper describes the results of a preliminary study to determine the effective heating and cooling time constants of ionospheric currents in a simulated modulated HF heating, ‘beam painting’ configuration. It has been found that even and odd harmonics of the fundamental ELF wave used to amplitude modulate the HF heater are sourced from different regions of the ionosphere which support significantly different heating and cooling time constants. The fundamental frequency and its odd harmonics are sourced in a region of the ionosphere where the heating and cooling time constants are about equal. The even harmonics on the other hand are sourced from regions of the ionosphere characterised by ratios of cooling to heating time constant greater than ten. It is thought that the even harmonics are sourced in the lower ionosphere (around 65 km) where the currents are much smaller than at the higher altitudes around 78 km where the currents at the fundamental frequency and odd harmonics maximise.     

含有旅行时间的水平反射界面上P-SV转换波转换点坐标的唯一解析解

转换波转换点位置的求取一直是一个引人关注的问题.Schneider（2002）从一个关于转换点的三次方程中推导出了P-SV转换波转换点位置的一种解析解,这种解是偏移距、横纵波速度和P-SV转换波旅行时的函数.我们严格地论证了Schneider（2002）给出的解析解的唯一性.这个解析解可以被看作为水平反射层P-SV转换波叠前偏移的精确算子.     

斜入射平面SH波在层状半空间中沉积谷地周围的三维散射

采用直接刚度法计算自由场动力响应,以层状半空间中移动斜线均布荷载动力格林函数模拟散射波场,采用间接边界元方法求解了层状半空间中沉积谷地对斜入射平面SH波的三维散射问题.由于文中采用的层状场地三维动力刚度矩阵是精确的,且用于模拟散射波场的均布移动斜线荷载可以直接施加在沉积交界面处而不存在奇异性,所以本文方法具有很高精度.文中以均匀半空间和基岩上单一土层中沉积谷地对入射SH波的散射为例进行了数值计算,研究表明,沉积谷地对地震波的三维散射与二维散射之间存在本质差别;层状半空间中沉积谷地与均匀半空间中沉积谷地附近地表位移存在显著差异.     

瑞雷波斜入射下层状半空间中沉积谷地周围的三维散射研究

针对层状半空间中沉积谷地对斜入射瑞雷波的三维散射问题, 采用直接刚度法计算自由场波场, 以层状半空间中移动斜线均布荷载动力格林影响函数求解三维散射波场, 建立了求解该问题的间接边界元方法. 通过与已有结果的比较, 验证了该方法的正确性, 并以均匀半空间以及弹性基岩上单一土层场地中沉积谷地为例进行了计算分析. 研究结果表明： 层状半空间与均匀半空间中沉积谷地对瑞雷波的散射存在显著差别; 层状半空间中瑞雷波的振动模态对沉积附近位移幅值有着重要影响; 土层刚度和厚度等参数也对沉积附近位移幅值大小及空间分布有着显著的影响.     

On the similarity in the formation mechanism of the fracture structure of a rock massif

Local segments of the Earth’s crust reside at the mechanical nonequilibrium and continuously obtain and dissipate mechanical energy. The energy exchange between the structural elements of a geophysical medium determines its state, especially if the medium is fragmented into blocks. The stationary state of a rock differs from its static equilibrium by the fact that the mechanical energy is conserved because the energy input is equal to the energy dissipation. From this point of view, the cracks, the faults, and the block structure are not simply the manifestations of rock destruction, but rather the mode of existence of a medium with large irreversible deformations. Then, the fracture structure, whose formation is actually a response of a rock to large irreversible deformation, becomes, in terms of physics, the characteristic of the state of a geophysical medium; for example, it allows one to assess the parameters of the deformation processes during the period of formation of the fracture structure. The present paper addresses the identification of the features of the fracture structure in geological objects of different scales.     

横波分裂裂缝检测方法在川西深层气藏中的应用

由横波分裂原理可知,当横波(或转换横波)在裂缝介质中传播时,横波将分裂成平行裂缝走向的快横波和垂直裂缝走向的慢横波.基于此,求取裂缝走向的问题可以转化为求取快横波的偏振方位角度的问题;由于快、慢横波的传播时差与裂缝发育密度密切相关,因而裂缝发育密度的预测问题可以转化为快、慢横波的传播时差求取问题.根据快、慢横波在径向分...     

Spatial structure of the magnetic field induced by an infrasonic wave in the oceanic waveguide

The magnetic field generated by an acoustic wave propagating in the oceanic waveguide has been considered. The effect of the self-induction factor on the spatial structure of this field has been studied. It has been indicated that there exists a frequency range where it is necessary to take into consideration self-induction. It has been indicated that the induced field is most substantial at frequencies when only the first normal mode exists. The dependences of the induced field on the depth, frequency, and geomagnetic field direction have been obtained and analyzed for this frequency range.     

3D scattering of obliquely incident plane SV waves by an alluvial valley embedded in a fluid-saturated,poroelastic layered half-space

The indirect boundary element method is used to study the 3D dynamic response of an infinitely long alluvial valley embedded in a saturated layered half-space for obli- quely incident SV waves. A wave-number transform is first applied along the valley＇s axis to reduce a 3D problem to a 2D plane strain problem. The problem is then solved in the section perpendicular to the axis of the valley. Finally, the 3D dynamic responses of the valley are obtained by an inverse wave-number transform. The validity of the method is con- firmed by comparison with relevant results. The differences between the responses around the valley embedded in dry and in saturated poroelastic medium are studied, and the effects of drainage conditions, porosity, soil layer stiffness, and soil layer thickness on the dynamic response are dis- cussed in detail resulting in some conclusions.