Numerical Simulation of Non-Linear Wave Propagation in Waters of Mildly Varying Topography with Complicated Boundary

In this paper, the characteristics of different forms of mild slope equations for non-linear wave are analyzed, and new non-linear theoretic models for wave propagation are presented, with non-linear terms added to the mild slope equations for non-stationary linear waves and dissipative effects considered. Numerical simulation models are developed of non-linear wave propagation for waters of mildly varying topography with complicated boundary, and the effects are studied of different non-linear corrections on calculation results of extended mild slope equations. Systematical numerical simulation tests show that the present models can effectively reflect non-linear effects.     

中国西部地区基于经验模型的土壤水分反演

地表土壤水分含量的时空分布信息是十分重要的,常常作为水文模型、气候模型、生态模型的输入参数,同时,也是干旱预报、农作物估产等工作的重要指标。被动微波遥感是监测土壤含水量最有效的手段之一。相比红外与可见光,它具有波长长,穿透能力强的优势。相比主动微波雷达,被动微波辐射计具有监测面积大、周期短,受粗糙度影响小,对土壤水分更为敏感,算法更为成熟的优势。目前,已研究出许多反演土壤水分的方法.本课题的主要内容是借助AMSR-E土壤水分影像数据、MODIS归一化植被指数(NDVI)影像数据和MODIS分类影像数据,利用ENVI软件进行遥感图像数据处理,运用统计分析方法建立NDVI与土壤水分的经验模型,研究中国西部地区稀疏植被覆盖区土壤水分的反演。     

静态精密单点定位精度分析

阐述了精密单点定位技术的应用价值,分析了观测数据和精密钟差采样间隔对静态精密单点定位的影响,提供了工程实际使用有益的参考.     

南海北部深水盆地沉积-构造的差异性及其油气意义

南海北部深水区自西向东依次分布着琼东南盆地、珠江口盆地、台西南盆地等新生代被动陆缘盆地,这些盆地经历了大致相当的从裂陷到坳陷的构造演化史,但在张裂活动过程中存在着明显的沉积-构造的差异性。构造沉降特征分析显示:在同一构造带上自西向东有盆地主要构造沉降发生的时段逐步变晚的趋势;在不同构造带上自北向南有盆地主要构造沉降发生的时段逐步变晚的趋势。这种沉积-构造的差异性对烃源岩的发育类型、分布及生储盖组合等方面有明显的控制作用,表现为:裂谷期构造沉降幅度大的盆地,陆相烃源岩发育,以陆生陆储陆盖型成藏组合为主;裂后期构造沉降幅度大的盆地,海相烃源岩发育规模较大,海生海储海盖型成藏组合及混生海储海盖型生储盖组合所占分量逐渐增多。推测渐新统湖相-湖沼相及海陆过渡相源岩和中新统海相烃源岩应是南海北部深水区油气的主要来源,陆生海储海盖型、海生海储海盖型及混生海储海盖型生储盖组合应是深水区基本生储盖组合类型。     

GPS精密单点定位的数据处理方法综述

在GPS精密单点定位数据处理中,通常采用基于非差观测值的处理方法。然而精密单点定位也可通过观测值不同差分组合以减少某些未知参数,且在理论上与非差模型等价。文中介绍基于非差和差分模型的精密单点定位的几种数据处理方法,并对各种模型进行比较。     

A continuous/discontinuous Galerkin framework for modeling coupled subsurface and surface water flow

We consider conjunctive surface-subsurface flow modeling, where surface water flow is described by the shallow water equations and ground water flow by Richards’ equation for the vadose zone. Coupling between the models is based on the continuity of flux and water pressure. Numerical approximation of the coupled model using the framework of discontinuous Galerkin (DG) methods is formulated. In the subsurface, the local discontinuous Galerkin (LDG) method is used to approximate ground water velocity and hydraulic head; a DG method is also used to approximate surface water velocity and elevation. This approach allows for a weak coupling of the models and the use of different approximating spaces and/or meshes within each regime. A simplified LDG method based on continuous approximations to water head is also described. Numerical results that investigate physical and numerical aspects of surface–subsurface flow modeling are presented. This work was supported by National Science Foundation grant DMS-0411413.     

Analysis of a rigid footing lying on three-layered soil using the finite difference method

A numerical procedure is described for the analysis of the vertical deformation and the stress distribution of the strip footings on layered soil media. Three layers of soil with different stiffness are considered with the middle soil layer the thinnest and most stiff layer. The soil media is discretized and using the theory of elasticity, the governing differential equations are obtained in terms of vertical and horizontal displacements. These equations along with appropriate boundary and continuity conditions are solved by using the finite difference method. The vertical and horizontal displacements, strains and stresses are found at various nodes in the soil media. Parametric studies are carried out to study the effect of the placement depth of the middle soil layer, the relative ratios of the moduli of deformation of the soil layers on the vertical displacement of the footing and the vertical stress distribution. These studies reveal that the middle thin but very stiff layer acts like a plate and redistributes the stresses on the lower soft soil layer uniformly. The displacement on the top and bottom of the middle soil layer is almost the same showing that the compression of the middle layer is negligible as it is very stiff.     

Traction image method for irregular free surface boundaries in finite difference seismic wave simulation

In this study, we propose a new numerical method, named as Traction Image method, to accurately and efficiently implement the traction-free boundary conditions in finite difference simulation in the presence of surface topography. In this algorithm, the computational domain is discretized by boundary-conforming grids, in which the irregular surface is transformed into a 'flat' surface in computational space. Thus, the artefact of staircase approximation to arbitrarily irregular surface can be avoided. Such boundary-conforming gridding is equivalent to a curvilinear coordinate system, in which the first-order partial differential velocity-stress equations are numerically updated by an optimized high-order non-staggered finite difference scheme, that is, DRP/opt MacCormack scheme. To satisfy the free surface boundary conditions, we extend the Stress Image method for planar surface to Traction Image method for arbitrarily irregular surface by antisymmetrically setting the values of normal traction on the grid points above the free surface. This Traction Image method can be efficiently implemented. To validate this new method, we perform numerical tests to several complex models by comparing our results with those computed by other independent accurate methods. Although some of the testing examples have extremely sloped topography, all tested results show an excellent agreement between our results and those from the reference solutions, confirming the validity of our method for modelling seismic waves in the heterogeneous media with arbitrary shape topography. Numerical tests also demonstrate the efficiency of this method. We find about 10 grid points per shortest wavelength is enough to maintain the global accuracy of the simulation. Although the current study is for 2-D P-SV problem, it can be easily extended to 3-D problem.