三角形网格下求解二维浅水方程的和谐Godunov格式

为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。     

Refraction Effect in Geothermal Heat Flow Due to a 3-D Prism in a Two-Layered Earth

A mathematical model of the disturbance of the stationary geothermal field due to a three-dimensional perturbing body embedded near the surface of a two layered earth. The theoretical analysis is based on the generalized theory of the double-layer potential, similar to the boundary integral method likely to direct current geoelectricity problems. Special attention is paid to cases in which a pan of the body's surface is tangential to the planar boundaries of the substratum or of the surface of the earth The numerical calculations are performed for a three-dimensional prism using methods similar to the boundary element method (BEM). Numerous graphs are shown for the disturbance of the temperature field and of the heat flow on the surface of the earth and along vertical profiles (boreholes).     

基于探地雷达的混凝土无损检测反演成像方法

本文针对水利水电工程中的混凝土构件无损探测问题,提出了一种新的反演成像方法.传统方法主要依赖人工经验观察偏移探测数据,精度无法保证.这种方法以Maxwell方程的TM问题为数学模型,在混凝土的表面放置发射源与接收器,根据探地雷达的接收数据,采用同伦优化方法结合收敛速度较快的阻尼高斯牛顿方法作为反演算法,其成像结果更具定量化和可视性.该成像方法克服了传统雷达剖面图只能近似地反映混凝土中埋藏物及缺陷的深度、大小等缺点,不仅能够探测混凝土内部异物的位置,而且能够比较准确地确定异物的大小及属性.通过计算机模拟以及对实际资料的处理,验证了该反演成像结果直观、可靠、具有一定的应用价值.     

一类测井响应方程存在的必要条件

本文将测井分析中所用的绝大多数响应方程归并为一类,给出了这类方程的一般形式。同时,根据一个基本事实导出了这类方程存在的必要条件,从而建立了检验和评价它们的理论标准,并已得到了实验验证。     

Seismic rays in an ellipsoidal earth model: theoretical analysis and estimation of deviations due to the ellipticity

Summary. A first-order form of the Euler's equations for rays in an ellipsoidal model of the Earth is obtained. The conditions affecting the velocity law for a monotonic increase, with respect to the arc length, in the angular distance to the epicentre, and in the angle of incidence, are the same in the ellipsoidal and spherical models. It is therefore possible to trace rays and to compute travel times directly in an ellipsoidal earth as in the spherical model. Thus comparison with the rays of the same coordinates in a spherical earth provides an estimate of the various deviations of these rays due to the Earth's flattening, and the corresponding travel-time differences, for mantle P -waves and for shallow earthquakes. All these deviations are functions both of the latitude and of the epicentral distance. The difference in the distance to the Earth's centre at points with the same geocentric latitude on rays in the ellipsoidal and in the spherical model may reach several kilometres. Directly related to the deformation of the isovelocity surfaces, this difference is the only cause of significant perturbation in travel times. Other differences, such as that corresponding to the ray torsion, are of the first order in ellipticity, and may exceed 1 km. They induce only small differences in travel time, less than 0.01s. Thus, we show that the ellipticity correction obtained by Jeffreys (1935) and Bullen (1937) by a perturbational method can be recovered by a direct evaluation of the travel times in an ellipsoidal model of the Earth. Moreover, as stated by Dziewonski & Gilbert (1976), we verify the non-dependence of this correction on the choice of the velocity law.     

基于非稳态匹配的角度域叠前道集去调谐方法

地层调谐效应改变了不同偏移距处反射波的干涉模式,使得远偏移距数据振幅和频率信息发生畸变,将会降低AVO/AVA分析及反演的置信度和分辨率.文中首先给出了调谐作用下反射波合成记录,通过正演模拟说明了调谐对叠前道集的影响.考虑到不同地层厚度下调谐效应的差异性,借助局域Lamoureux窗实现地震数据的自适应分解,利用角度数据之间的差异,构建了角度域地震数据的非平稳匹配目标函数,形成了非平稳匹配去调谐方法,实现了叠前道集的振幅和波形拉伸校正.本文方法在实际应用中取得了较好的效果,能够有效地改善大角度地震数据品质,为储层预测与流体识别奠定了数据基础.     

引入PLSA模型的光学遥感图像舰船检测

提出一种基于概率潜在语义分析(Probabilistic Latent Semantic Analysis,PLSA)的检测算法,首先通过PLSA将目标表述为潜在成分的概率组合,然后利用统计模式识别方法对获取的潜在成分概率进行判别,从而完成最终的检测。其中,生成的潜在成分反映了目标与特征之间相互出现的频率关系,并以潜在成分在目标中概率差异的形式对上述不对应现象给出了直观描述。实验结果表明,所提出的算法对多种复杂情况下的光学图像舰船检测具有很好的准确性和鲁棒性。     

One‐dimensional dynamics of saturated incompressible porous media: analytical solutions and influence of inertia terms

The complexity of formulations for the hydromechanical coupled mechanics of porous media is typically minimised by simplifying assumptions such as neglecting the effect of inertia terms. For example, three formulations commonly employed to model practical problems are classified as fully dynamic, simplified dynamic and quasi‐static. Thus, depending on the porous media conditions, each formulation will have advantages and limitations. This paper presents a comprehensive analysis of these limitations when solving one‐dimensional fully saturated porous media problems in addition to a new solution that considers a more general loading situation. A phase diagram is developed to assist on the selection of which formulation is more appropriate and convenient regarding particular cases of porosity and hydraulic conductivity values. Non‐dimensional formulations are proposed to achieve this goal. Results using the analytical solutions are compared against numerical values obtained with the finite element method, and the effect of porosity is investigated. Copyright © 2016 John Wiley & Sons, Ltd.     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

地心坐标反解大地坐标的非线性方程组牛顿迭代法分析

目前的一些地心坐标向大地坐标的转换模型在计算速度,稳定性或精度方面存在一定的局限性。本文探讨了非线性方程组数值迭代的求解方法及其在MATLAB 7.0中的编程实现,并将结果与基于一元三次方程求解的严密方法做了比较分析,结果证明该方法是可以在实际中参考应用的。