Comparison of methods to model the gravitational gradients from topographic data bases

A number of methods have been developed over the last few decades to model the gravitational gradients using digital elevation data. All methods are based on second-order derivatives of the Newtonian mass integral for the gravitational potential. Foremost are algorithms that divide the topographic masses into prisms or more general polyhedra and sum the corresponding gradient contributions. Other methods are designed for computational speed and make use of the fast Fourier transform (FFT), require a regular rectangular grid of data, and yield gradients on the entire grid, but only at constant altitude. We add to these the ordinary numerical integration (in horizontal coordinates) of the gradient integrals. In total we compare two prism, two FFT and two ordinary numerical integration methods using 1" elevation data in two topographic regimes (rough and moderate terrain). Prism methods depend on the type of finite elements that are generated with the elevation data; in particular, alternative triangulations can yield significant differences in the gradients (up to tens of Eötvös). The FFT methods depend on a series development of the topographic heights, requiring terms up to 14th order in rough terrain; and, one popular method has significant bias errors (e.g. 13 Eötvös in the vertical–vertical gradient) embedded in its practical realization. The straightforward numerical integrations, whether on a rectangular or triangulated grid, yield sub-Eötvös differences in the gradients when compared to the other methods (except near the edges of the integration area) and they are as efficient computationally as the finite element methods.     

Numerical simulation research on multi-electrodes resistivity imaging survey array

Multi-electrodes Resistivity Imaging Survey （MRIS） is an array method of electrical survey. In practice how to choose a reasonable array is the key to get reliable survey results. Based on four methods of MRIS such as Wenner, Schlumberger, Pole-pole and Dipole-dipole the authors established the model, by studying the result of the forward numerical simulation modeling and inverse modeling, and analyzed the differences among the different forms of detection devices.     

A Comparative Study of Genetic Variation betwccn Chinese Mitten Crab Eriocheir Sinensis and Hepu Mitten Crab E. Hepuensis

Horizontal starch gel electrophoresis was used to investigate the biochemical genetic structure of Chinese mittencrab Eriocheir sinensis and Hepu mitten crab E. hepuensis. Sixteen putative enzyme-coding loci were examined for forty-nineChinese mitten crabs and thirty-eight Hepu mitten crabs. Nine loci, AAT-1 *, AAT-2 *, G3PDH*, GPI*, IDHP-1 *, IDHP-2 *,MDH-1 *, MDH-2 * and PGM * , are polymorphic in Chinese mitten crab, and seven, AAT-1 *, AAT-2 *, GPI *, IDHP-1 *,MDH-1 *, MDH-2 * and PGM * , are polymorphic in Hepu mitten crab. The proportion of polymorphic loci and the expectedheterozygosity are 0.562 5 and 0.080 3 for Chinese mitten crab, and 0.437 5 and 0.075 4 for Hepu mitten crab. The Nei' s ge-netic distance between two species is 0.002 4.     

龙门山及其邻区的地壳厚度和泊松比

根据龙门山及其周边地区(26°～35°N,98°～109°E)的132个台站的宽频带远震记录,使用H-k叠加方法计算地壳厚度和波速比.结果表明该区域的地壳厚度总体变化是:从东向西增加,东部的最小厚度为37.8km,西部的最大厚度是68.1 km,其中横跨龙门山断裂带的地壳厚度变化最大,从东南的41.5km增加到西北的52.5km.根据Airy均衡理论,用台站的高程和观测地壳厚度数据求得最小二乘意义下的壳幔密度差为0.649g/cm3,平均地壳厚度为37.9km.龙门山及其邻近地区基本上处于均衡状态.松潘-甘孜地体北部和西秦岭造山带具有低泊松比(v<0.26),扬子地台的西南部具有低一中泊松比(v<0.27),松潘-甘孜地体南部、三江褶皱带和四川盆地具有中一高泊松比(0.26≤P≤0.29).该地区的泊松比空间分布不支持青藏高原东部广泛分布的下地壳流的假说.龙门山断裂带南段及其附近地区的高泊松比(v≥0.30)可以看成是地壳具有较高的铁镁质组分和/或存在部分熔融.该地区下地壳可能是处于富含流体和温度较高的部分熔融状态.松潘-甘孜块体南部的上地壳物质向东运动,受刚性强度较大的扬子地台的阻挡,导致沿龙门山断裂带产生应变积累.当断层被地壳流体弱化,积累的应变能量快速释放,产生汶川Ms8.0地震.     

库车坳陷克拉2气田异常流体压力演化史

根据克拉2气田现今实测地层压力分析了异常高压发育的特征,并通过流体包裹体测试求取包裹体形成时期的古压力。计算结果表明:下白垩统和古近系库姆格列木群砂岩储层在距今8～6Ma(康村组沉积期末)时,古压力系数大致为1.09～1.18,发育过剩压力为4.6～9.6MPa的异常高压。在此基础上,以古、今压力数据为约束条件,利用数值模拟技术定量恢复出下白垩统与侏罗系异常流体压力的演化历史。研究发现:下白垩统储层流体压力经历了沉积型异常流体压力形成发展(白垩纪-库车组沉积期末)和沉积型异常流体压力萎缩—构造挤压型异常流体压力孕育(库车组沉积期末至今)两个阶段。康村组沉积早期沉积型超压开始出现,库车组沉积期末达到顶峰,过剩压力达40～50MPa;此后,在沉积型超压萎缩—构造挤压型超压孕育过程的耦合下,过剩压力曾降至25～30MPa,更新世起超压略有升高。侏罗系流体压力演化较为简单,吉迪克组沉积期开始出现超压,至库车组沉积期末达到高峰,此后异常高压逐渐萎缩。在克拉2气田的主要成藏时期(距今3～1Ma),烃源岩较储层具有更高的异常压力,源—储压力差为10～30MPa。     

地下结构对工程场地地震动场的影响

本文以日本神户市地铁大开站为工程背景,应用大型通用有限元ANSYS软件,分别建立了自然条件下的工程场地土层和建有地下车站结构的土层的二维有限元分析模型。计算了基岩一致地震动输入模式下2个土层模型的动力反应。通过将2个土层模型的动力反应进行比较分析,讨论了地下结构对工程场地地震动场的影响。分析结果表明,日本神户地铁大开站的存在对其邻近区域土层的地震反应有影响,但影响有限。     

GeoPW: Laying Blocks for the Geospatial Processing Web

Recent advances in Web‐related technologies have significantly promoted the wide sharing and integrated analysis of distributed geospatial data. Geospatial applications often involve diverse sources of data and complex geoprocessing functions. Existing Web‐based GIS focuses more on access to distributed geospatial data. In scientific problem solving, the ability to carry out geospatial analysis is essential to geoscientific discovery. This article presents the design and implementation of GeoPW, a set of services providing geoprocessing functions over the Web. The concept of the Geospatial Processing Web is discussed to address the geoprocessing demands in the emerging information infrastructure, and the role of GeoPW in establishing the Geospatial Processing Web is identified. The services in GeoPW are implemented by developing middleware that wraps legacy GIS analysis components to provide a large number of geoprocessing utilities over the Web. These services are open and accessible to the public, and they support integrated geoprocessing on the Web.     

离子迁移谱现场观测渤海和北黄海二甲基硫的研究

渤海、黄海是高产二甲基硫（Dimethyl Sulfide, DMS）的大陆架海区。该海区DMS的现场调查研究有助于准确评估海洋DMS释放量及其对全球气候变化的负反馈作用。目前,无论是基于模型还是直接测量法的通量估算均以表层海水或低层大气DMS浓度为基础,因此,先进的检测技术对其通量估算的准确度具有决定性作用。气相色谱法、质谱法、化学发光法以及卫星遥感技术是现在常用的观测技术,而本文则基于苯辅助光电离离子迁移谱技术进一步提出了一种可在海域现场观测海水中DMS的方法。通过结合动态气提-Nafion管在线除水进样系统,消除环境水汽的干扰;在最优条件下,基于DMS两个产物离子峰,可以实现0.10～120 nmol/L之间DMS的定量分析,检测限低至0.065 nmol/L;然后将所建方法应用于2019年秋季渤海、北黄海海水中DMS的现场观测。结果表明,表层海水中DMS的浓度为0.080～0.96 nmol/L（平均值为（0.44±0.34）nmol/L）,其海气通量为0.12～17.75 μmol/（m²·d）（平均值为（ 3.23±4.02）μmol/（m²·d））;通过结合实验室检测结果、环境因子和浮游植物群落结构讨论了海水样品低温储存条件下DMS的变化和影响因素,结果显示,营养盐成分及浮游植物群落结构是影响储存样品中DMS浓度显著增加的主要因素,进一步表明了现场观测方法的建立对海洋DMS释放量的准确评估具有重要意义。     

滩涂水边线Landsat-5影像提取方法以长江口崇明东滩为例

快速获取遥感影像水边线时空变化信息具有重要意义,滩涂水边线提取一直是遥感技术应用的难点问题。水边线在遥感影像上具有独特的空间关系与光谱特征。本文综合使用颜色模型变换法、信息熵计算法、最大类间方差法及边缘检测方法。以长江口崇明东滩为研究区,研究了Landsat-5卫星影像海陆对比度增强及不同尺度下的边缘提取,重点给出了基于热红外波段的水边线空间特征与光谱特征的计算方法,在面向对象技术框架下提出了一种顾及空间关系和光谱特征的遥感影像水边线快速提取方法。实验结果表明：(1)基于最大类间方差法的局部阈值分割法能够自动提取band 6的水边线,水边线连续、完整,空间信息丰富;(2)综合使用最佳指数法、离散度方法及颜色模型变换方法,能够有效增强海陆对比度,基于最大类间方差法的局部自适应Canny算子能够自动检测出增强后遥感影像高精度边缘;(3)利用水边线的空间关系和光谱特征,能够由计算机自动实现水边线的识别与连接工作;(4)本文提出的水边线提取方法速度快、自动化程度高,分别继承了阈值分割法的连续性强的优点和Canny算子定位精度高、细节呈现能力强的优势。研究结果对于海岸带动态变化、陆海相互作用机制、海岸带资源保护与开发及近海工程管理等研究具有重要的参考价值。     

微扰Kepler系统轨道微分方程的近似Lie对称性与近似不变量

把极角θ视为独立变量,得到Kepler系统的轨道微分方程. 首先讨论Kepler系统轨道微分方程的Lie对称性和不变量,微扰Kepler系统轨道微分方程的精确Lie对称性和精确不变量,其次讨论微扰Kepler系统轨道微分方程的近似Lie对称性和近似不变量,并得到了微扰Kepler系统的9个一阶近似Lie对称性和6个一阶近似不变量,其中1个实为精确不变量,而其余5个分别等于微扰系数ε乘以Kepler系统相应的5个不变量.