Unsteady source flow in weakly heterogeneous porous media

Average nonuniform flows in heterogeneous formations are modeled with the aid of the nonlocal effective Darcy's law. The mean head for flow toward source of instantaneous discharge in a heterogeneous medium of given statistics represents the fundamental solution of the average flow equation and is called the Mean Green Function (MGF). The general representation of the MGF is obtained for weakly heterogeneous formations as a functional of the logconductivity correlation function. For Gaussian logconductivity correlation, the MGF is derived in terms of one quadrature in time t and it is analyzed for isotropic media of any dimensionality d and for 3D axisymmetric formations. The MGF is further applied to determining the mean head distribution for flow driven by a continuous source of constant discharge. The large time asymptotic of the mean head is analyzed in details.     

新型边界积分方程的无单元解及其应用

目前隧道及大型地下工程往往在裂隙岩体中开挖,而裂隙与地下空间的距离及裂隙的扩展条件,制约着隧道及地下工程的稳定性。应用能考虑孔洞和裂纹问题的新型边界积分方程与无网格加辽金法结合,建立一种新型的边界无单元法。在该方法中基本的未知量是由边界上的面力和边界上位移密度函数构成的复变量边界函数 。文中应用的边界积分公式和Muskhelishvili的积分公式直接相关。将无网格构造方法引入新型的边界积分方程,建立了新型的边界无单元法。应用该方法详细分析了含隧道和裂纹间相互关系等问题,其数值结果与解析结果吻合很好,说明该方法的正确性和可行性。     

正交各向异性岩体孔径变形法的二维地应力分析方法改进

根据正交各项异性材料的弹性理论,推导出偏轴坐标系下的三维、二维本构方程;利用应力变分法对钻孔周围的应力分布作出修正,进一步得出孔壁径向位移的计算公式。结合孔径变形法测量地应力的特点,建立正交各向异性岩体地应力的分析理论,通过一典型工程与传统的各向同性理论进行了对比。     

TTI介质弹性波频率-空间域有限差分数值模拟

由周期性薄互层引起的VTI介质是研究比较广泛的一类各向异性介质。当VTI介质对称轴偏离垂向,本构坐标系与观测坐标系不重合时,会形成观测坐标系下的TTI介质。引入25点优化差分算子,推导出二维TTI介质频率域弹性波动方程;为压制边界反射,采用完全匹配层法吸收边界条件,并计算出优化差分系数;最后采用集中力源,模拟了弹性波在TTI介质中的传播过程。从波场快照和地面共炮记录可以看出,笔者采用的数值模拟算法能有效压制数值频散。TTI介质中的波场传播比较复杂,纵波传播相对稳定,横波波前的三分叉现象比较明显,并存在振幅奇异性。当VTI介质的对称轴偏转后,还会增加地面地震记录的复杂性。     

结合抛物型缓坡方程计算波浪辐射应力

将波浪辐射应力与抛物型缓坡方程中的待求变量联系起来,提出了一种计算辐射应力的新方法,并用有限差分法对控制方程进行了数值求解。数值结果表明这种方法精度高、编程简单、求解快速,可用于实际大区域波浪辐射应力的计算。     

An Extended Mild-Slope Equation

On the assumption that the vortex and the vertical velocity component of the current aresmall,a mild-slope equation for wave propagation on non-uniform flows is deduced from the basichydrodynamic equations,with the terms of (V_hh)~2 and (V_h~2)h included in the equation.The terms of bot-tom friction,wind energy input and wave nonlinearity are also introduced into the equation.The wind en-ergy input functions for wind waves and swells are separately considered by adopting Wen′s(1989)empiri-cal formula for wind waves and Snyder′s observation results for swells.Thus,an extended mild-slope equa-tion is obtained,in which the effects of refraction,diffraction,reflection,current,bottom friction,wind en-ergy input and wave nonlinearity are considered synthetically.     

海浪场非平衡信息论熵的研究初探

鉴于传统熵与时间无关 ,不能反映海 -气耦合机制下海浪场的基本属性———不可逆性 ,文章将非平衡信息论熵的概念引入海浪场的研究。依据海浪场相空间的Fokker Planck方程 ,首次推导出了非平衡信息论熵的熵平衡方程、熵密度、熵流密度、熵产生密度、广义力和广义力密度的具体表达式。     

Effects of Topography and Current on Horizontal Irrotational Waves in Shallow Water

Based on the Boussinesq assumption,derived are couple equations of free surface elevationand horizontal velocities for horizontal irrotational flow,and analytical expressions of the correspondingpressure and vertical velocity.After the free surface elevation and horizontal velocity at a certain depth areobtained by numerical method,the pressure and vertical velocity distributions can be obtained by simplecalculation.The dispersion at different depths is the same at the O(ε)approximation.The waveamplitude will decrease with increasing time due to viscosity,but it will increase due to the matching ofviscosity and the bed slope.thus,flow is unstable.Numerical or analytical results show that the waveamplitude.velocity and length will increase as the current increases along the wave direction.but theamplitude will increase.and the wave velocity and length will decrease as the water depth decreases.     

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation is briefly des-cribed.Several currently used wave numerical models for shallow water based on different wave theoriesare discussed.Wave energy conservation models for the simulation of shallow water waves are introduced,with emphasis placed on the SWAN model,which takes use of the most advanced wave research achieve-ments and has been applied to several theoretical and field conditions.The characteristics and applicabilityof the model,the finite difference numerical scheme of the action balance equation and its source termscomputing methods are described in detail.The model has been verified with the propagation refractionnumerical experiments for waves propagating in following and opposing currents;finally.the model is ap-plied to the Haian Gulf area to simulate the wave height and wave period field there,and the results arecompared with observed data.     

二维各向异性介质中三分量波动方程有限元法模拟

利用六面体单元和三线性插值函数 ,推导出在任意弹性各向异性介质中三维三分量波动方程所满足的有限元方程 ,同时给出也适用于矩形单元和双线性插值的二维三分量各向异性介质中的波动方程模拟的一般性公式 ,最后讨论两个数值模拟的结果。