变差函数参数的计算

为了求出地质统计学中变差函数理论模型的参数，首先把相应模型化成线性函数的形式，多面手应用线性方程组非负解的理论计算出它的参数。     

结构内设备的摩擦摆隔震响应分析

基于系统耦和振动微分方程，分析了结构内摩擦摆隔震设备的响应规律和隔震效果。计算结果表明，经过合理设计的摩擦摆系统能够显著控制设备的地震绝对加速度响应从而有效提高设备的抗震可靠度。     

AN APPARENT TEMPERATURE MODEL BASED ON THERMAL EQUILIBRIUM

In this paper,a comprehensive analysis and review of the advances in the research of human body comfort degree is made.Based on the studies of human thermal equilibrium theories contributed by previous researchers,different apparent temperature models corresponding to different environments are set up in accordance with the current situation and the characteristics of metabolism in different groups of people in China.After case studies and comparison with present dominant human body comfort degree statistical models,it is proved that the apparent temperature models have high rationality and wide adaptability.Furthermore,the comfort scaling standard according to apparent temperature is suggested,which is suitable for the middle latitude regions in China.     

基于矩形网格的有限差分走时计算方法

对于大多数速度场，地震波沿射线传播的初至波走时，可以用有限差分外推的方法在二维或三维数值网格上计算出来. 在保证精度的条件下，为提高计算效率和适应性，本文推导了基于任意矩形网格和局部平面波前近似的有限差分初至波走时计算方法. 另外，该方法对首波和散射波做了合适的处理，而且不会碰到传统射线法存在的阴影区和焦散区等问题. 简单模型和复杂的Marmousi模型试算的结果表明，该方法精度较高并适用于强纵、横向变速的复杂介质. 基于该方法的Kirchhoff叠前深度偏移, 在主要构造和目的层位置的成像效果上基本达到了波动方程法叠前深度偏移的位置成像效果. 由于未考虑续至波等有效能量，在成像的保幅性上不如波动方程法叠前深度偏移的效果，但其计算效率则明显高于全格林函数法和波动方程法.      

变降深抽水附近地下水非稳定运动的积分方程解

根据完整井变流量抽水附近地下水运动的解析解,利用积分方程的方法,推导出了水井水位降深以任意一种函数关系随时间变化时含水层中地下水非稳定运动积分方程形式的解,并提出了该积分解相应的近似计算方法.最后通过一个实例验证了该方法的实用性.     

SH Wave Propagation in Viscoelastic Media

When comparing solutions for the propagation of SH waves in plane parallel layered elastic and viscoelastic (anelastic) media, one of the first things that becomes apparent is that in the elastic case the location of the saddle points required to obtain a high frequency approximation are located on the real p axis. This is true of the branch points also. In a viscoelastic medium this is not typical. The saddle point corresponding to an arrival lies in the first quadrant of the complex p-plane as do the branch points. Additionally, in the elastic case the saddle point and branch points lie on a straight line drawn through the origin (the positive real axis in the complex p-plane), while in the viscoelastic case this is generally not the case and the saddle point and branch points lie in such a manner as to indicate the degree of their complex values.In this paper simple SH reflected and transmitted particle displacement arrivals due to a point torque source at the surface in a viscoelastic medium composed of a layer over a half space will be considered. The path of steepest descent defining the saddle point in the first quadrant will be parameterized in terms of a real variable and the high frequency solutions and intermediate analytic results obtained will be used to formulate more specific constraints and observations regarding saddle point location relative to branch point locations in the complex p-plane.As saddle point determination for an arrival is, in general, the solution of a non-linear equation in two unknowns (the real and imaginary parts of the complex saddle point p ₀), which must be solved numerically, the use of analytical methods for investigating this problem type is somewhat limited.Numerical experimentation using well documented solution methods, such as Newton's method, was undertaken and some observations were made. Although fairly basic, they did provide for the design of algorithms for the computation of synthetic traces that displayed more efficient convergence and accuracy than those previously employed. This was the primary motivation for this work and the results from the SH problem may be used with minimal modifications to address the more complicated subject of coupled P-SV wave propagation in viscoelastic media.Another reason for revisiting a problem that has received some attention in the literature was to approach it in a fairly comprehensive manner so that a number of specific observations may be made regarding the location of the saddle point in the complex p-plane and to incorporate these into computer software. These have been found to result in more efficient algorithms for the SH wave propagation and a significant enhancement of the comparable software in the P-SV problem.     

时域瞬变场电磁场有限差分法

模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。     

非饱和土壤Richards方程入渗求解探讨

根据理查兹 (Richards)方程求解非饱和土壤垂直一维入渗率历时规律 ,因缺乏土壤水分特性函数K(θ)、D(θ)的一般有效公式 ,而采用与实际相差甚远的定值或线性函数假设 ,使成果隐含过多经验成分或需要反复迭代计算。笔者在文献 [1]中提出的非饱和土壤水分特性函数的一般有效基础公式 ,用分离变量法直接求解Richards方程 ,用土壤水动力学方法分别导出了含有渗前土湿因子的新入渗公式。它摆脱了同类入渗公式对水分函数过于简单的概化假设和经验成分。为放弃迄今水文学产流计算中惯用的下渗容量曲线“截首留尾”的经验做法 ,以直接使用含有渗前土湿因子的下渗公式提供了理论依据 ,从而将有利于产流计算方法的改进和计算精度的提高。     

降雨条件下岩土饱和-非饱和渗流分析

基于Buckley-Leverett两相渗流方程,提出了新的岩土饱和-非饱和渗流数学模型,利用有限差分方法给出了隐式压力显示饱和度数值求解方法,编制了饱和-非饱和渗流计算程序.结合工程实例,模拟了降雨入渗条件下边坡岩体渗流场孔隙压力变化和含水饱和度变化,模拟结果验证了所提出的模型对饱和-非饱和渗流的有效性.     

川西亚高山三个森林群落的湿林冠蒸发速率

于3个假说和林冠上方2m处的气象变量，采用Penman-Monteith组合模型估算了一个生长季节内川西亚高山林区分别以云杉(SF)、冷杉(FF)和白桦(BF)为优势树种的3个林分的湿林冠蒸发速率(Er)。研究结果表明，SF、FF和BF的湿林冠蒸发量(E)分别为44．51mm、88．51mm和57．8mm，分别占总降雨量的9．2％、16．6％和10．2％。与SF和BF相比，FF具有最高的月平均Er和蒸发比例。SF、FF和BF的平均Er分别为0．097mm／h(变化范围：0．028-0．487mm／h)、0．242mm／h(变化范围：0．068～0．711mm／h)和0．149mm／h(0．060～0．576mm／h)。最高和最低的月平均Er分别在6月(SF、FF和BF分别为0．120mm／h、0．317mm／h和0．169mm／h)和10月(SF、FF和BF分别为0．083mm／h、0．187mm／h和0．101mm／h)。8：00至16：00期间的平均点Er显著高于0：00至8：00以及16：00至0：00期间的平均Er。Er显著的日变化和月变化主要归因于林冠上方的太阳辐射、空气温度和相对湿度的变化。