流体饱和多孔隙垂向非均匀固体的弹性波传播理论

建立含有可压缩粘滞流体的多孔隙垂向非均匀弹性固体内之应力波传播理论。论述了其中的一个问题。即当Ｐｏｉｓｅｕｉｌｌｅ流动假设成立的低频率范围且流体与固体的质量密度可相比拟（例如含水饱和岩石之类）时，流体与固体间虽有相对运动但无摩擦的情形。在一定的近似条件下，可以存在两类膨胀波和一类旋转波，它们都存在频散     

分层密度流体的边界波动在不透水刚性底床的运动行为——Ⅰ造波条件对波型演化的关系

在一座二维内波实验水槽(12 m×0.7 m×0.5 m)、于总水深H等于40 cm的条件下,以淡、盐水依密度分层布置的两层流体,从事孤立内波生成的研究。内波水槽内的上下层水深比H₁/H₂、造波区位能差(η₀)、上下层的密度比皆为实验主要的控制变因。在各种不同实验条件配合之下,试验及比较以上变因对内波振幅能量及相关传递系数的影响,进而产生几种不同的波动特性:上举型孤立内波、下沉形孤立内波、非周期性的重力内波等,且不同形式的重力内波需要配合固定的造波条件使能达成。     

饱和多孔介质中的非均匀平面波

讨论了两相饱和多孔介质中的非均匀平面波,在所采用的多孔介质模型中,利用了每种组分的微观密度保持不变而与体积分数相关的宏观密度是变化的假设。针对理想无粘流体充填的弹性多孔固体,深入讨论了两种非均匀平面波的传播特性以及质点运动的轨迹线。     

行星构造:寻求地球演化的踪迹

地质构造是记录地球内、外动力地质作用过程的标志.和地球相似,太阳系其他天体上也发育丰富的地质构造.以研究天体表面的地质构造及其动力学机制为目的的"行星构造学"是建立在构造地质学、遥感地质学和地球物理学等学科基础上的一门新兴前沿学科.由于天体的大小、组分和轨道位置不同,表面构造特征及其形成机制各异.对比研究地球和其他天体...     

A semi-analytical method for the wave-induced seabed response

A general semi-analytical method is presented for the analysis of seabed response under wave action during a storm. The seabed is idealized as a poro-elastic medium filled with a single compressible fluid with anisotropic flow. The coupled process of fluid flow and deformation of soil skeleton is formulated in the framework of Biot's theory. The analysis for the response of homogeneous seabed of finite thickness under a plane progressive wave is developed first, followed by an extension for the case of a layered seabed. A generalization for three-dimensional response of seabed is also developed for a general wave field which provides the analyses for seabed response under short-crested or standing waves in the vicinity of a structure. Some numerical examples illustrating the proposed analyses are also presented.     

在沿岸流作用下的线性短峰波解析解

由于目前三维短峰波的研究甚少涉及普遍的波-流相互作用机制,为此,针对波浪倾斜入射岸堤、反射、并伴随强烈沿岸流作用的典型近海波况,构造出有限水深短峰波一般线性系统,并得到其一般解析解和特定的深水、浅水波解析解。从而,进一步完善了现有的线性水波理论体系,为非线性短峰波理论在波-流相互作用领域里的探讨提供了一个可靠基础。     

Large transverse bedforms and the character of boundary-layers in shallow-water environments

Theoretical work, laboratory studies, and field observations indicate that the oscillatory boundary layers generated by the tidal wave differ fundamentally in dynamics and kinematic structure from the unidirectional boundary layers of rivers. Unique to the former are mass-transport currents attributable to: (1) the wave motion itself, and (2) bed curvature in the presence of the oscillatory flow. The implication of this difference for bed-material transport is that the larger flow-transverse bedforms of shallow-water environments are divisible hydraulically between two major classes: (A) those related to tidal conditions, under which the fluid reverses in direction of flow with each reversal of the tide, permitting the initiation and maintenance of bed features by the spatially reversing, curvature-related mass transport, and (B) those related strictly to rivers and river-like flows, in which the fluid motion is unidirectional, and therefore the only mechanisms available for bedform initiation and maintenance are those creating a finite spatial lag between the transport rate and the bed waviness. Forms of Class B are best called dunes and bars, and only those attributable to Class A should be termed sand waves. The latter, restricted to oscillatory boundary layers of tidal origin, apparently correspond to the very much smaller; but also commonly symmetrical, ripple marks produced in wind-wave oscillatory boundary layers.     

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

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

Lava worlds: From early earth to exoplanets

The magma ocean concept was first conceived to explain the geology of the Moon, but hemispherical or global oceans of silicate melt could be a widespread “lava world” phase of rocky planet accretion, and could persist on planets on short-period orbits around other stars. The formation and crystallization of magma oceans could be a defining stage in the assembly of a core, origin of a crust, initiation of tectonics, and formation of an atmosphere. The last decade has seen significant advances in our understanding of this phenomenon through analysis of terrestrial and extraterrestrial samples, planetary missions, and astronomical observations of exoplanets. This review describes the energetic basis of magma oceans and lava worlds and the lava lake analogs available for study on Earth and Io. It provides an overview of evidence for magma oceans throughout the Solar System and considers the factors that control the rocks these magma oceans leave behind. It describes research on theoretical and observed exoplanets that could host extant magma oceans and summarizes efforts to detect and characterize them. It reviews modeling of the evolution of magma oceans as a result of crystallization and evaporation, the interaction with the underlying solid mantle, and the effects of planetary rotation. The review also considers theoretical investigations on the formation of an atmosphere in concert with the magma ocean and in response to irradiation from the host star, and possible end-states. Finally, it describes needs and gaps in our knowledge and points to future opportunities with new planetary missions and space telescopes to identify and better characterize lava worlds around nearby stars.     

Wave propagation in a general anisotropic poroelastic medium: Biot’s theories and homogenisation theory

Anisotropic wave propagation is studied in a fluid-saturated porous medium, using two different approaches. One is the dynamic approach of Biot’s theories. The other approach known as homogenisation theory, is based on the averaging process to derive macroscopic equations from the microscopic equations of motion. The medium considered is a general anisotropic poroelastic (APE) solid with a viscous fluid saturating its pores of anisotropic permeability. The wave propagation phenomenon in a saturated porous medium is explained through two relations. One defines modified Christoffel equations for the propagation of plane harmonic waves in the medium. The other defines a matrix to relate the relative displacement of fluid particles to the displacement of solid particles. The modified Christoffel equations are solved further to get a quartic equation whose roots represent complex velocities of the four attenuating quasi-waves in the medium. These complex velocities define the phase velocities of propagation and quality factors for attenuation of all the quasi-waves propagating along a given phase direction in three-dimensional space. The derivations in the mathematical models from different theories are compared in order to work out the equivalence between them. The variations of phase velocities and attenuation factors with the direction of phase propagation are computed, for a realistic numerical model. Differences between the velocities and attenuations of quasi-waves from the two approaches are exhibited numerically.     

关于深水环境下内波、内潮汐沉积分类的探讨

深水环境下内波、内潮汐沉积从在地层记录中被发现已过去20年,在这20多年里虽然积累了不少资料,但至今未对内波、内潮汐沉积进行细分。本文以地层记录中已发现的内波、内潮汐沉积为基础,结合内波、内潮汐的破碎过程和海洋物理学中内波的研究现状,对深水内波、内潮汐沉积的分类进行了探索。将其分为3个层次下的9种类型,分别为：正压内潮汐沉积、斜压内潮汐沉积、正压短周期内波沉积、斜压短周期内波沉积、等深流叠加内波沉积、低密度浊流叠加内波沉积、长周期内波叠加沉积、驻波沉积和其他内波叠加沉积。地层记录中已发现的内波、内潮汐沉积分别归为：正压内潮汐沉积、正压短周期内波沉积和长周期内波叠加沉积3类,其余类型尚未被发现。这样的分类较合理地解释了为什么在现代海洋中内波、内潮汐无处不在,但地层记录中的内波、内潮汐沉积却少得可怜的现象,同时将内波、内潮汐沉积研究和海洋物理学中对内波的研究结合起来,并能将内波、内潮汐沉积研究置于更为广阔的研究背景之中,使之可以和大洋突发事件、天文旋回及大洋环流等联系起来。     

Mudbank regime off the Kerala coast during monsoon and non-monsoon seasons

Field experiments conducted in the nearshore ocean to understand the dynamics of mudbank off Kerala, south-west coast of India, are highlighted. Real time monitoring of the nearshore ocean off Purakkad, Kerala was accomplished using pressure transducers for nearshore surface wave measurements, and current sensors for nearshore velocity measurements. Comprehensive information on the spatial structure of mudbank was obtained from aerial surveys. Extensive data collected on surface waves and currents in the nearshore ocean, indicate that the infra-gravity (IG) waves (leaky modes and trapped edge wave modes), and far infra-gravity (FIG) waves coupled with strong shoreline reflections and undertow play an important role in the dynamics associated with the mudbanks off Kerala during the monsoon season. During the non-monsoon season evidence for progressive edge waves in the infragravity frequency band, an energetic gravity wave band and a strong undertow with weak reflections was observed.     

Investigating origin of the inadequate medium range predictability of the lower tropospheric ultra-long waves in tropics

The predictability of planetary or ultra-long scale waves is limited by the large growth of errors in these scales in almost all the medium range forecast models. Understanding the cause for the enormous build up of error is, therefore, a necessary task for improving the prediction of planetary waves. A diagnostic analysis of the systematic error energetics has been performed in the Global Forecast System model to investigate the reasons for poor predictability of the lower tropospheric ultra-long waves (wavenumber bands 1–4) in tropics using the analysis-forecast system of horizontal wind field at 850 hPa level during the boreal summer period. For this purpose, systematic error energy is computed in spatial as well as in wavenumber domain. Non-linear inter-scale transfer of error has been formulated and evaluated through energy exchanges among participating triads. The study reveals that the error is generated in the prognostic model initially with a small magnitude at the different locations around tropical convergence zone (TCZ) attributed to the inaccuracy in representing different physical processes like cumulus convection applied in the model. At subsequent evolution of forecasts, error increases and spreads along the TCZ due to its non-linear advection to the higher scales and eventually to the ultra-long scales attributed to the inherent dynamics of the model evaluated through the process of wave-wave exchange of error energy in terms of the triad interactions. The continuous generation and then, non-linear propagation of error up to the planetary scales in the course of prediction increase the uncertainty in ultra-long scales which actually inhibit to predict accurately the planetary scale waves in tropics during medium range forecasts. This work suggests caveats to the modeler’s community in the predictability study of tropical ultra-long waves.     

Wave propagation in thermoelastic saturated porous medium

Biot’s theory for wave propagation in saturated porous solid is modified to study the propagation of thermoelastic waves in poroelastic medium. Propagation of plane harmonic waves is considered in isotropic poroelastic medium. Relations are derived among the wave-induced temperature in the medium and the displacements of fluid and solid particles. Christoffel equations obtained are modified with the thermal as well as thermoelastic coupling parameters. These equations explain the existence and propagation of four waves in the medium. Three of the waves are attenuating longitudinal waves and one is a non-attenuating transverse wave. Thermal properties of the medium have no effect on the transverse wave. The velocities and attenuation of the longitudinal waves are computed for a numerical model of liquid-saturated sandstone. Their variations with thermal as well as poroelastic parameters are exhibited through numerical examples.     

A non-hydrostatic model for the numerical study of landslide-generated waves

This paper proposes and demonstrates a two-layer depth-averaged model with non-hydrostatic pressure correction to simulate landslide-generated waves. Landslide (lower layer) and water (upper layer) motions are governed by the general shallow water equations derived from mass and momentum conservation laws. The landslide motion and wave generation/propagation are separately formulated, but they form a coupled system. Our model combines some features of the landslide analysis model DAN3D and the tsunami analysis model COMCOT and adds a non-hydrostatic pressure correction. We use the new model to simulate a 2007 rock avalanche-generated wave event at Chehalis Lake, British Columbia, Canada. The model results match both the observed distribution of the rock avalanche deposit in the lake and the wave run-up trimline along the shoreline. Sensitivity analyses demonstrate the importance of accounting for the non-hydrostatic dynamic pressure at the landslide-water interface, as well as the influence of the internal strength of the landslide on the size of the generated waves. Finally, we compare the numerical results of landslide-generated waves simulated with frictional and Voellmy rheologies. Similar maximum wave run-ups can be obtained using the two different rheologies, but the frictional model better reproduces the known limit of the rock avalanche deposit and is thus considered to yield the best overall results in this particular case.     

Numerical and analytical investigation of compressional wave propagation in saturated soils

In geotechnical earthquake engineering, wave propagation plays a fundamental role in engineering applications related to the dynamic response of geotechnical structures and to site response analysis. However, current engineering practice is primarily concentrated on the investigation of shear wave propagation and the corresponding site response only to the horizontal components of the ground motion. Due to the repeated recent observations of strong vertical ground motions and compressional damage of engineering structures, there is an increasing need to carry out a comprehensive investigation of vertical site response and the associated compressional wave propagation, particularly when performing the seismic design for critical structures (e.g. nuclear power plants and high dams). Therefore, in this paper, the compressional wave propagation mechanism in saturated soils is investigated by employing hydro-mechanically (HM) coupled analytical and numerical methods. A HM analytical solution for compressional wave propagation is first studied based on Biot’s theory, which shows the existence of two types of compressional waves (fast and slow waves) and indicates that their characteristics (i.e. wave dispersion and attenuation) are highly dependent on some key geotechnical and seismic parameters (i.e. the permeability, soil stiffness and loading frequency). The subsequent HM Finite Element (FE) study reproduces the duality of compressional waves and identifies the dominant permeability ranges for the existence of the two waves. In particular the existence of the slow compression wave is observed for a range of permeability and loading frequency that is relevant for geotechnical earthquake engineering applications. In order to account for the effects of soil permeability on compressional dynamic soil behaviour and soil properties (i.e. P-wave velocities and damping ratios), the coupled consolidation analysis is therefore recommended as the only tool capable of accurately simulating the dynamic response of geotechnical structures to vertical ground motion at intermediate transient states between undrained and drained conditions.     

地下洞室群对地面运动影响问题的级数解答:SH波入射

采用波函数展开法给出平面SH波入射下半空间中洞室群对地面运动影响问题的一个级数解答，并分析了入射波频率和入射角以及两洞室之间距离等参数对地面运动的影响。数值结果表明，地下洞室群对沿线地面运动具有显著的放大作用：两个洞室情况地面运动可以达到单个洞室情况的1．6倍、无洞室情况的3．5倍以上。建议地铁等地下工程在设计和规划时应考虑工程建设后对沿线设计地震动的影响。     

海洋生物地球化学模式研究进展

海洋生物地球化学模式是定量认识物质的海洋生物地球化学循环、理解其控制机制以及预测体系变动的重要手段。20世纪90年代以来,该研究领域的进展主要体现在海洋生物地球化学循环的物理输送和生态动力学过程以及年际、年代际变动的模拟3个方面。物理过程模拟方面的进展,集中在寡营养海区上层海水营养盐的供应机制问题上,在经典的上升流、垂直扩散之外,提出涡旋可能构成一种重要的物理输入过程。而生态动力学过程的模拟方面,90年代前期考虑食物网基本结构,由浮游植物、浮游动物和细菌三大类群构成生物状态变量,氮和磷营养盐以及颗粒碎屑构成其他状态变量;90年代后期,开始引入铁和硅的限制问题,考虑不同浮游植物和浮游动物群落结构的影响,特别是浮游植物粒级结构变化的预测可能是未来该领域力图解决的一个技术问题。年际变化的模拟,多围绕ENSO事件对初级生产的影响及其机制问题展开;年代际和地质年代尺度的体系变动问题仍存在争论,相对缺乏有效的数值模拟研究。该研究领域未来应加强生物—化学过程的函数表达、物理模式、中尺度过程、边界交换以及资料获取技术等方面的研究,以应对目前面临的诸多问题与挑战。     

Analytical Study of Underground Explosion-Induced Ground Motion

A theoretical approach is presented to study the ground motion induced by an underground tunnel explosion. The ground motion is caused by two coupled stress waves, i.e., the reflected body wave and the secondary surface wave or Raleigh wave. Based on the principle of conservation of momentum at the wavefronts, the reflected body waves along the ground surface are derived. The interaction of the body wavefront and the ground surface induces the secondary surface wave which transfers outwards on the ground. The particle velocity and particle acceleration on the ground surface are subsequently derived. The analytical results are compared with results from numerical simulations and empirical formulae with different material damping ratios. The effects of the loading density and the material damping on the ground motion are investigated. Finally, the limitations of the proposed theoretical approach for ground motion prediction are discussed.     

半空间饱和土中圆形衬砌对弹性压缩波的散射

借助于Biot 波动理论和弹性波的传播理论,采用复变函数和多级坐标法,对半空间饱和土中圆形衬砌结构对弹性稳态压缩波的散射问题进行求解和分析。利用一个半径很大的圆弧来逼近半空间直边界,将待解问题转化为稳态弹性压缩波在一个大圆孔和一个弹性衬砌结构的散射问题。通过引入势函数,将饱和土的Biot波动方程和衬砌的弹性波动方程解耦成Helmholtz 方程,借助复变函数级数展开便可以预先写出该组Helmholtz方程的通解。然后,通过引用复变量,把饱和土和衬砌结构中的应力、位移及孔压用设定的势函数表示出来,再利用半空间饱和土和衬砌结构的连续性条件和近似直边界的圆弧边界和衬砌内边界的边界条件求解出该组势函数的特解。最后,利用势函数的特解,得到饱和土中的位移,应力和孔压及衬砌结构的位移和应力;变换不同的参数求解衬砌结构内外边界的动应力和孔压的集中系数,通过对算例结果的分析得出一系列有益的结论。