Three‐phase medium model for filled rock joint and interaction with stress waves

A three‐phase medium model is proposed in describing the dynamic property of filled rock joints and an analytical study on longitudinal wave transmission normally across a three‐phase rock joint is presented. Parameters in the three‐phase medium model were determined by a series of modified split Hopkinson pressure bar (SHPB) tests, where a sand or clay layer was used to represent an artificially filled rock joint. The effect of the unloading path on the transmitted wave was discussed by comparing the analytical and SHPB test results. The derived wave transmission coefficients across the filled joint agreed very well with those from the test results. Both the analytical and the test results showed that the wave transmission coefficients were affected by the mechanical properties of the fillings. Parametric studies with respect to the volume ratios of water and air in the three‐phase medium and the type of filling material have also been performed. Copyright © 2010 John Wiley & Sons, Ltd.     

Exact solution for granular flows

In this paper, we present the exact solution of the Riemann problem for the nonlinear one‐dimensional so‐called shallow‐water or Saint‐Venant equations with friction proposed by SAVAGE and HUTTER to describe debris avalanches. This model is based on the depth‐averaged thin layer approximation of granular flows over sloping beds and takes into account a Coulomb type friction law with a constant friction coefficient. A particular configuration of the Riemann problem corresponds to a dam of infinite length in one direction from which granular material is released from rest at a given time over an inclined rigid or erodible bed. We solve analytically and numerically the depth‐averaged long‐wave equations derived in a topography‐linked coordinate system for all the possible Riemann problems. The detailed mathematical proof of the derivation of the analytical solutions and the analysis of their structure and properties is intended, first of all, for geophysicists, mathematicians, and physicists because of the possible extension of this study to more complex problems (geometries, friction laws, …). The numerical solution of the first‐order finite‐volume method based on a Godunov‐type scheme is compared with the proposed exact Riemann problem solution. This solution is used to solve the dam‐break problem and analyze the influence of the thickness of the erodible bed on the speed of the granular front. Comparison with existing experimental results shows that, for an erodible bed, the equations lack fundamental physical significance to reproduce the observed dynamics of erosive granular flows. Copyright © 2012 John Wiley & Sons, Ltd.     

Propagation of Love wave in fiber‐reinforced medium lying over an initially stressed orthotropic half‐space

In this paper, the propagation of Love waves in a fiber‐reinforced layered medium lying over an elastic orthotropic half‐space under initial stress has been investigated. We have obtained the velocity equations for Love wave in this media. It is observed that propagation of Love wave is influenced by reinforced parameters and initial stressed parameter. The velocity of Love wave has been computed for three different cases. Our computed equation of Love wave coincides with the standard equation of Love wave for the case of homogeneous layer and homogeneous half‐space (AEH Love, 1911). To study the effect of reinforced and initial stressed parameters, we have computed the numerical values for phase velocity and plotted in several figures. It is observed that the phase velocity decreases with the increases of reinforced parameters and initial stressed parameter. Using MATLAB software, GUI has been developed to generalize the effect of various parameters. Copyright © 2014 John Wiley & Sons, Ltd.     

Exact analytical solutions of nonlinear problems of tsunami wave run-up on slopes with different profiles

A review of papers investigating tsunami wave run-up on a beach is given and the control parameters of the problem are revealed. There are two such parameters in the case of ideal fluid: the bottom sloping angle and the breaking parameter. A stage-by-stage approach for finding run-up characteristics is formulated: the linear calculation of shoreline oscillations and the subsequent non-linear transformation of the solution according to the Riemann method. Solution of the nononedimensional problems of wave run-up on a beach in the linear formulation is obtained.     

Source signature and elastic waves in a half‐space under a momentary shear line impulse

The transient deformation of an elastic half‐space under a line‐concentrated impulsive vector shear load applied momentarily is disclosed in this paper. While in an earlier work, the author gave an analytical–numerical method for the solution to this transient boundary‐value problem, here, the resultant response of the half‐space is presented and interpreted. In particular, a probe is set up for the kinematics of the source signature and wave fronts, both explicitly revealed in the strained half‐space by the solution method. The source signature is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Fourteen wave fronts exist behind the precursor shear wave S: four concentric cylindrical, eight plane, and two relativistic cylindrical initiated at propagating centres that are located on the stationary boundaries of the solution domain. A snapshot of the stressed half‐space reveals that none of the 14 wave fronts fully extend laterally. Instead, each is enclosed within point bounds. These wave arresting points and the two propagating centres of the relativistic waves constitute the source signature. The obtained 14 wave fronts are further combined into 11 disparate wave fronts that are grouped into four categories: an axis of symmetry wave—so named here by reason of being a wave front that is contiguous to the axis of symmetry, three body waves, five surface waves and two inhibitor waves—so named here by reason that beyond them the material motion dies out. Of the three body waves, the first is an unloading shear wave, the second is a diffracted wave and the third is a reflected longitudinal two‐branch wave. Of the two inhibitor waves, the first is a two‐joint relativistic wave, while the second is a two‐branch wave. The wave system, however, is not the same for all the dependent variables; a wave front that appears in the behaviour of one dependent variable may not exist in the behaviour of another. It is evident from this work that Saint–Venant's principle for wave propagation problems cannot be formulated. Therefore, the above results are valid for the particular proposed model for the momentary line‐concentrated shear load. The formulation of the source signature, the wave system, and their role in the half‐space transient deformation are presented here. Copyright © 2003 John Wiley & Sons, Ltd.     

可压缩流体饱和孔隙介质中孔隙压力波传播数值分析

首次将高精度时空守恒元/解元方法推广到可压缩流体饱和孔隙介质中孔隙压力波传播的数值计算中。将孔隙度梯度从源（汇）项中分离,直接引入流通量,改进了理论模型。通过对孔隙介质激波问题的数值模拟,验证了方法的精度和有效性。在此基础上,提出了孔隙介质中二维黎曼问题,并揭示了孔隙压力波存在接触间断、激波、膨胀波、压缩波等复杂的结构特征。该成果对二氧化碳地质封存、二氧化碳提高石油采收率、页岩气压裂开采以及地震破裂过程的研究具有重要的理论与应用意义。     

Source signature and elastic waves in a half‐space under a sustainable line‐concentrated impulsive normal force

First, the response of an ideal elastic half‐space to a line‐concentrated impulsive normal load applied to its surface is obtained by a computational method based on the theory of characteristics in conjunction with kinematical relations derived across surfaces of strong discontinuities. Then, the geometry is determined of the obtained waves and the source signature—the latter is the imprint of the spatiotemporal configuration of the excitation source in the resultant response. Behind the dilatational precursor wave, there exists a pencil of three plane waves extending from the vertex at the impingement point of the precursor wave on the stress‐free surface of the half‐space to three points located on the other two boundaries of the solution domain. These four wave‐arresting points (end points) of the three plane waves constitute the source signature. One wave is an inhibitor front in the behaviour of the normal stress components and the particle velocity, while in the behaviour of the shear stress component, it is a surface‐axis wave. The second is a surface wave in the behaviour of the horizontal components of the dependent variables, while the third is an inhibitor wave in the behaviour of the shear stress component. An inhibitor wave is so named, since beyond it, the material motion is dying or becomes uniform. A surface‐axis wave is so named, since upon its arrival, like a surface wave, the dependent variable in question features an extreme value, but unlike a surface wave, it exists in the entire depth of the solution domain. It is evident from this work that Saint‐Venant's principle for wave propagation problems cannot be formulated; therefore, the above results are a consequence of the particular model proposed here for the line‐concentrated normal load. Copyright © 2002 John Wiley & Sons, Ltd.     

Cambrian wave‐dominated tidal‐flat deposits,central Wisconsin,USA

In central Wisconsin, Cambrian strata of the Elk Mound Group record deposition on open‐coast, wave‐dominated tidal flats. Mature, medium‐grained quartz arenite is dominated by parallel‐bedding with upper‐flow regime parallel‐lamination, deposited during high‐energy storms that also produced three‐dimensional bedforms on the flats. Abundant wave ripples were produced as storms waned or during fair weather, in water depths ≤2 m. Indicators of variably shallow water (washout structures and stranded cnidarian medusae) and subaerial exposure (adhesion marks, rain‐drop impressions and desiccation cracks, including cracked medusae) are abundant. Parallel‐bedded facies preserve a Cruziana ichnofacies, similar to other Cambrian tidal‐flat deposits. Flats were dissected by small, mainly straight channels, the floors of which were grazed intensely by molluscs. Most channels were ephemeral but some developed low levées, point bars and cut‐banks, probably reflecting stabilization by abundant microbial mats and biofilms. Channels were filled with trough cross‐bedding that is interpreted to have been produced mainly during storm runoff. The strata resemble deposits of open‐coast, wave‐dominated tidal flats on the east coast of India and west coast of Korea. Ancient wave‐dominated and open‐coast tidal flats documented to date appear to have been limited to mud‐rich strata with ‘classic’ tidal indicators such as flaser bedding and tidal bundles. The Cambrian (Miaolingian to early Furongian) Elk Mound Group demonstrates that sandy, wave‐dominated tidal flats also can be recognized in the stratigraphic record.     

Spatial and Temporal Features and Trend of Different Level Heat Waves Over China

Heat wave has become a severe problem over China in recent year. Based on daily air temperature data from 719 meteorological stations in China in the period of 1961-2012, the heat wave events were classified into three levels by duration: weak HW, medium heat wave and strong heat wave. This paper describes the spatial and temporal characteristics among three levels of heat wave over China during the period 1961-2013. The results showed that the number of hot days displayed a weakening trend from the 1960s to the early 1990s, followed by a strengthening trend from the late 1990s up to now. Long-term linear trends in hot days had significant regional differences, for instance, the number of hot days in southwest and southern China growing sharply than other regions; three levels heat wave occurred in different regions that weak heat wave usually occurred in the northwest and middle part of China while strong heat wave mainly happened in the southern part of China. When it came to the first day of heat wave, weak heat wave occurred earliest in China, and it happened around July 3rd while medium heat wave and strong heat wave always started on July 13th and July 14th, respectively. In general, the earliest heat wave was always found in southwest China while the last ones usually occurred in southern part and southwest China. Increasing trend of three levels heat wave were found in the entire country, except for southwest China. The trend of the last day of heat waves showed significant difference between northern and southern part of China. The last day of weak and medium heat wave occurred later in the south China, while strong heat wave happened earlier in north China.     

Torsional wave propagation in non‐homogeneous layer between non‐homogeneous half‐spaces

The study of surface wave in a layered medium has their possible application in geophysical prospecting. In the present work, dispersion equation for torsional wave in an inhomogeneous isotropic layer between inhomogeneous isotropic half‐spaces has been derived. Two cases are discussed separately for torsional wave propagation in inhomogeneous layer between homogeneous and non‐homogeneous half‐spaces, respectively. Further, two possible modes for torsional wave propagation are obtained in case of inhomogeneous layer sandwiched between non‐homogeneous half‐spaces. Closed form solutions for displacement in the layer and half‐spaces are obtained in each case. The study reveals that the layer width, layer inhomogeneity, frequency of inhomogeneity, as well as inhomogeneity in the half‐space has significant effect on the propagation of torsional surface waves. Displacement and implicit dispersion equation for torsional wave velocities are expressed in terms of Heun functions and their derivatives. Effects of inhomogeneity on torsional wave velocity are also discussed graphically by plotting the dimensionless phase velocity against dimensionless and scaled wave number for different values of inhomogeneity parameter. Copyright © 2012 John Wiley & Sons, Ltd.     

Coupled simulation of wave propagation and water flow in soil induced by high‐speed trains

The purpose of this paper is to simulate the coupled dynamic deformation and water flow that occur in saturated soils when subjected to traffic loads, which is a problem with several practical applications. The wave propagation causes vibrations leading to discomfort for passengers and people in the surroundings and increase wear on both the vehicle and road structure. The water flow may cause internal erosion and material transport in the soil. Further, the increased pore water pressure could reduce the bearing capacity of embankments. The saturated soil is modelled as a water‐saturated porous medium. The traffic is modelled as a number of moving wheel contact loads. Dynamic effects are accounted for, which lead to a coupled problem with solid displacements, water velocity and pressure as primary unknowns. A finite element program has been developed to perform simulations. The simulations clearly demonstrate the induced wave propagation and water flow in the soil. The simulation technique is applicable to railway as well as road traffic. Copyright © 2007 John Wiley & Sons, Ltd.     

爆破荷载应力波在无限弹性介质中传播的特征线法解

双指数爆破荷载的应力时程关系是爆破工程的数值分析常用的动力输入。其在弹性岩石中传播衰减的特征线法解是工程爆破数值分析中弹性以及弹塑性分析的重要参照,也是动力学数值分析算法开发的重要高阶校准工具。应用统一弹性波的特征线法解求解源自圆形孔洞的双指数爆破荷载在理想介质中的传播。介绍了源自圆形孔洞的双指数爆破荷载在周围弹性介质传播过程中的径向应力、周向应力、位移以及速度波的特征线法解的具体过程,在MATLAB中进行了4种波的特征线法解的数值实施。对径向应力、周向应力、位移以及速度波的特征线法解的结果进行了多方面讨论。结果表明,特征线法是求解双曲型偏微分方程的一种有效方法,双指数爆破荷载在弹性岩石中传播衰减的特征线法解可为爆破工程数值分析以及动力学数值分析算法开发提供重要的支撑。     

转换波综合静校正技术在新场地区的应用

转换波静校正问题一直被认为是转换波处理中较难解决的问题,目前已研究的转换波静校正方法都存在一定的不足。这里根据新场地区转换波资料的特点,采用多种静校正方法相结合的方式,构建了切实可行的转换波综合静校正技术。首先,利用模型法求取转换波长波长静校正量;然后,采用全局寻优静校正方法解决残留的长波长和中长波长转换波静校正问题;最后,在共检波点叠加道集上相关方法计算短波长静校正量。实际数据处理表明,该转换波静校正技术能较好地解决新场地区的转换波静校正问题,使地震记录质量得到明显改善。     

Wave propagation and strain localization in a fully saturated softening porous medium under the non‐isothermal conditions

The (THM) coupling effects on the dynamic wave propagation and strain localization in a fully saturated softening porous medium are analyzed. The characteristic polynomial corresponding to the governing equations of the THM system is derived, and the stability analysis is conducted to determine the necessary conditions for stability in both non‐isothermal and adiabatic cases. The result from the dispersion analysis based on the Abel–Ruffini theorem reveals that the roots of the characteristic polynomial for the THM problem cannot be expressed algebraically. Meanwhile, the dispersion analysis on the adiabatic case leads to a new analytical expression of the internal length scale. Our limit analysis on the phase velocity for the non‐isothermal case indicates that the internal length scale for the non‐isothermal THM system may vanish at the short wavelength limit. This result leads to the conclusion that the rate‐dependence introduced by multiphysical coupling may not regularize the THM governing equations when softening occurs. Numerical experiments are used to verify the results from the stability and dispersion analyses. Copyright © 2016 John Wiley & Sons, Ltd.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

The method of fundamental solution for 3‐D wave scattering in a fluid‐saturated poroelastic infinite domain

By using a complete set of poroelastodynamic spherical wave potentials (SWPs) representing a fast compressional wave P_I, a slow compressional wave P_II, and a shear wave S with 3 vectorial potentials (not all are independent), a solution scheme based on the method of fundamental solution (MFS) is devised to solve 3‐D wave scattering and dynamic stress concentration problems due to inhomogeneous inclusions and cavities embedded in an infinite poroelastic domain. The method is verified by comparing the result with the elastic analytical solution, which is a degenerated case, as well as with poroelastic solution obtained using other numerical methods. The accuracy and stability of the SWP‐MFS are also demonstrated. The displacement, hoop stress, and fluid pore pressure around spherical cavity and poroelastic inclusion with permeable and impermeable boundary are investigated for incident plane P_I and SV waves. The scattering characteristics are examined for a range of material properties, such as porosity and shear modulus contrast, over a range of frequency. Compared with other boundary‐based numerical strategy, such as the boundary element method and the indirect boundary integral equation method, the current SWP‐MFS is a meshless method that does not need elements to approximate the geometry and is free from the treatment of singularities. The SWP‐MFS is a highly accurate and efficient solution methodology for wave scattering problems of arbitrary geometry, particularly when a part of the domain extends to infinity.     

用遗传算法进行瑞利波反演

瑞利波反演属于非线性最优化问题,已有的局部线性化方法(如阻尼最小二乘法)常使解估计陷入目标函数的局部极小值之中,且严重地依赖初始模型.另外,目前被广泛应用的一些简单的瑞利波反演方法也存在明显的缺陷.全局优化法-遗传算法大大放松了对初始模型选择的要求,且不易陷入局部最优解中.本文用已有的简单宜行的反演解释方法获取初始模型,从而确定模型参数的搜索范围,再用遗传算法反演得到最终的介质模型,效果非常理想.     

非饱和土中热弹性波的传播特性分析

基于非饱和多孔介质的研究成果,考虑热效应和孔隙流体迂曲度的影响,研究了非饱和土中热弹性波的传播特性。利用非饱和土中耦合热的固-液-气三相介质的质量平衡方程、渗流连续方程、动量平衡方程和广义非Fourier热传导定律,建立了问题的热弹性波动方程。通过引入势函数,经过理论推导给出了非饱和土中热弹性波的弥散特征方程。结合数值算例,分析了几类热弹性波的波速和衰减系数随迂曲度、热膨胀系数和介质温度等热物理参数的变化规律。结果表明：孔隙水迂曲度的增大将引起P1波、P3波和S波的波速增大,而孔隙气体迂曲度的增大仅使得P2波的波速增大;热膨胀系数的增大将造成P1波波速的增大和热（T）波波速的减小;介质温度的升高将引起各类热弹性波波速的增大;频率、热膨胀系数和介质温度的变化对各类热弹性波的衰减系数均有较大影响,不可忽视。     

非饱和土地基中P1波通过波阻板的传播特性研究

基于单相弹性介质与非饱和多孔介质中波的传播理论和Snell定律,研究了P1波从非饱和土介质中入射到波阻板（wave impeding block,WIB）时的传播特性。推导了P1波在非饱和土地基中通过波阻板后反射波和透射波振幅比的解析解,分析了反射波振幅比和透射波振幅比与波阻板的剪切模量、密度以及入射角、入射频率、饱和度等各种参数的关系曲线,从而得到这些参数对非饱和土地基中波阻板隔振效果的影响规律。数值结果表明：波阻板隔振屏障的剪切模量和密度对其隔振性能影响较为显著,当波阻板剪切模量和密度在一定范围内时波阻板隔振屏障可获得较好的隔振性能;入射角对波阻板隔振性能影响较大,而入射频率对其隔振效果影响很小。