平面点集凸壳的一种快速算法

提出一种计算平面点集凸壳的快速算法———八方向极值快速凸壳算法。该算法首先对平面点集进行一次扫描,从而快速查找到东、南、西、北、东南、西南、东北、西北8个方向上的极值点,构造出一个更接近凸壳的初始凸壳,从而在后续的点集扫描中可以排除更多的内点,使该算法计算效率更高。该算法的空间复杂度为O(N);其时间复杂度虽然无法突破最坏情况下O(NlogN)的理论下限,但其期望时间复杂度已达到线性水平,并且可以容易地扩展到三维和高维空间。     

Implementation of the semi-Lagrangian advection scheme on a quasi-uniform overset grid on a sphere

The semi-Lagrangian advection scheme is implemented on a new quasi-uniform overset （Yin-Yang） grid on the sphere. The Yin-Yang grid is a newly developed grid system in spherical geometry with two perpendicularly-oriented latitude-longitude grid components （called Yin and Yang respectively） that overlapp each other, and this effectively avoids the coordinate singularity and the grid convergence near the poles. In this overset grid, the way of transferring data between the Yin and Yang components is the key to maintaining the accuracy and robustness in numerical solutions. A numerical interpolation for boundary data exchange, which maintains the accuracy of the original advection scheme and is computationally efficient, is given in this paper. A standard test of the solid-body advection proposed by Williamson is carried out on the Yin-Yang grid. Numerical results show that the quasi-uniform Yin-Yang grid can get around the problems near the poles, and the numerical accuracy in the original semi-Lagrangian scheme is effectively maintained in the Yin-Yang grid.     

利用毛管压力曲线分形分维方法研究流动单元

利用取心井铸体薄片获得的图像资料和毛管压力曲线,通过图像分形几何学方法以分维数的形式定量地表征出了复杂的微观孔隙喉道结构特征,发现能够很好地划分和评价孔隙岩石中油、气、水的渗流差异,可以用于储层微观流动单元表征。文中阐述了岩石微观孔隙喉道结构分形的理论基础、计算方法和应用于表征流动单元的依据。建立了中国西部砾岩低渗透油藏微观孔隙喉道分维数与孔隙度、渗透率之间计算图版,据此在油藏中利用常规测井资料获得的孔隙度、渗透率参数计算微观孔隙喉道分维数,开展全油藏流动单元划分与评价,取得了良好的效果。研究结果表明,利用毛管压力曲线分形分维方法研究储层微观流动单元是一种很有效的途径。     

长江三峡巴东复杂斜坡系统成因研究

三峡库区巴东县新城区所在的扇形大斜坡是一个复杂斜坡系统。斜坡后缘边界的挤压变形带是斜坡整体沿其底界面（T1j^3-T2b^1）之间岩层界面向长江发生重力蠕滑作用的结果,不存在“巴东断裂”问题。斜坡系统内部未发现成规模的褶皱,也不存在数量众多的内动力成因的断裂带。巴东斜坡区浅表生变形破坏类型可概括为侧向拉裂、“雪崩式”垮塌、表层风化、侧向滑移张剪破裂、岩层折断变位、软岩膝折剪破裂和压扭破裂又重胶结7种重力卸荷破坏力学模式。斜坡系统内部发育分期分区（块）滑动（滑坡）现象,但不存在构造地质意义上的“断层”形迹。在地貌形态上,巴东大斜坡被4条冲沟分割成5个斜坡单元。在斜坡地质结构、成因类型与空间分布上,巴东大斜坡可划分为3个层次,即表层崩塌滑坡成因为主的堆积层（第一层次）,冲沟分割且浅表生地质形迹发育的层状顺倾的中间基岩层（第二层次）和整体连续顺倾的深层基岩层（第三层次）。分形几何计算证明,巴东大斜坡的地貌形态尚处于侵蚀发育的青（幼）年期。FLAC^3D数值模拟发现,在长江侵蚀下切的不同阶段,巴东斜坡体前缘和后缘接近底界面位置塑性变形区分布集中,但不具备沿深层界面发生整体滑动剪出的可能性。基本认识是,“巴东复杂斜坡系统”是在官渡口-东瀼口向斜南翼（单斜山）的地质背景下,持续经受长江快速侵蚀下切外动力作用,河谷岸坡快速临空导致其自身重力产生强烈的侧向卸荷与滑移等浅表生地质改造作用过程而形成的,可简单地概括为“重力成因论”。     

Simulating cross-sectional geometry of the main channel in response to changes in water and sediment in Lower Yellow River

Journal of Geographical Sciences - To understand the non-equilibrium morphological adjustment of a river in response to environmental changes, it is essential to (i) accurately identify how past...     

URANS investigation of the interaction between the free surface and a shallowly submerged underwater vehicle at steady drift

An axisymmetric underwater vehicle (UV) at a steady drift angle experiences the complex three-dimensional crossflow separation. This separation arises from the unfavorable circumferential pressure gradient developed from the windward side toward the leeward side. As is well known, the separated flow in the leeward side gives rise to the formation of a pair of vortices, which affects considerably the forces and moments acting on the UV. In this regard, the main purpose of the present study is to evaluate the role of the leeward vortical flow structure in the hydrodynamic behavior of a shallowly submerged UV at a moderate drift angle traveling beneath the free surface. Accordingly, the static drift tests are performed on the SUBOFF UV model using URANS equations coupled with a Reynolds stress turbulence model. The simulations are carried out in the commercial code STARCCM+ at a constant advance velocity based on Froude number equal to F_n = 0.512 over submergence depths and drift angles ranging from h = 1.1D to h = ∞ and from β = 0 to β = 18.11°, respectively. The validation of the numerical model is partially conducted by using the existing experimental data of the forces and moment acting on the totally submerged bare hull model. Significant interaction between the low-pressure region created by the leeward vortical flow structure and the free surface is observed. As a result of this interaction, the leeward vortical flow structure appears to be largely responsible for the behavior of the forces and moments exerted on a shallowly submerged UV at steady drift.     

Temporal variation of modulated-wave-train geometries and their influence on vertical bending moments of a container ship

A towing experiment was conducted using a modulated wave train to investigate the vertical bending responses of a hydro-structural container ship model. In the experiment, a spatially periodic modulated wave train, as a model of a freak wave in successive high waves mimicking the so-called three sisters, was generated by the recently established higher-order spectral method wave generation (HOSM-WG) method. HOSM-WG enables us to control the location and timing of the maximum crest height in a wave tank. With precise control of the towing carriage, an experiment was conducted in which the timing of the encounters between the ship model and the modulated wave train was accurately determined. The maximum sagging moment (SM) was found to increase in proportion with the encounter wave height. However, because of differences in the relative depth of the fore and aft troughs, the maximum SM is highly variable for a given wave height. The temporal wave-geometry evolution caused the relative trough-depth to vary significantly within a wave period in the vicinity of the maximum crest height. As a result, depending on the encounter timing, the SM varied considerably for a given wave height. The temporal variation of the wave geometry is a robust feature of a modulated wave train and is common between the spatially periodic and temporally periodic modulated wave trains.     

基于几何变换的图形采集与编辑算法研究

分析了解析几何法的符号绘制算法,采用了几何变换实现图形绘制的算法,实现了椭圆弧的采集、编辑及Bezier曲线和椭圆弧外接矩形的算法。     

FY-3微波成像仪遥感图像地理定位方法研究

MWRI (MicroWave Radiation Imager) is one of the payloads on our next generation polar meteorological satellite FY-3. MWRI conically scans with a fixed incident angle on the earth surface. It is the first time for Chinese remote sensor to use this scan mode. In this work, we present a geolocation method for FY-3 MWRI’s remote sensing image based on its special scan geometry. The integrated coordinate systems and the specific relationships with these coordinate systems are defined. A spatial relationship model between the remote sensing data and the earth-based coordinate system is established. This method also includes an algorithm of satellite orbit computation, which is used to get the satellite’s instantaneous velocity vector from its position. This method has been applied to MWRI’s remote sensing image geolocation. The results show that the accuracy of this method can achieve 1 pixel. The 33 GCPs (Ground Control Points) which are in the regiones of FY-3 MWRI’s observation have been collected and used to analyze the precision of the geolocation. By statistical analysis, the error along-track is about 1.5km, and the error along-scan is about 3.0km. It is obvious that this method fulfills the requirement of precision for FY-3 MWRI whose space resolution exceeds 5km.     

基于特征提取的目标变化检测方法研究

针对不同时相遥感影像中特定目标的变化,研究了基于特征提取的目标变化检测方法。该方法首先根据目标的几何特征进行检测与提取,然后对提取出的两目标影像进行比较,得到特定目标的变化。实验结果表明该方法对特定目标的变化检测有较好的效果。