Interpolation by regularized spline with tension: I. Theory and implementation

Bivariate and trivariate functions for interpolation from scattered data are derived. They are constructed by explicit minimization of a general smoothness functional, and they include a tension parameter that controls the character of the interpolation function (e.g., for bivariate case the surface can be tuned from a membrane to a thin steel plate), Tension can be applied also in a chosen direction, for modeling of phenomena with a simple type of anisotropy. The functions have regular derivatives of all orders everywhere. This makes them suitable for analysis of surface geometry and for direct application in models where derivatives are necessary. For processing of large datasets (thousands of data points), which are now common in geosciences, a segmentation algorithm with a flexible size of overlapping neighborhood is presented. Simple examples demonstrating flexibility and accuracy of the functions are presented.On leave from the Department of Physical Geography and Cartography, Comenius University, Mlynská dolina, Bratislava, Czechoslovakia.On leave from the Institute of Physics, Dúbravská cesta 9, Bratislava, Czechoslovakia.     

相律的应用:蒸发结晶过程自由度

在连续原理和浓度限制条件的基础上,分析了卤水蒸发过程的自由度,提出了“蒸发结晶过程自由度”这一概念。发现蒸发结晶过程自由度在整个过程中不是1便是0,在1和0之间交替变化,与选择的组分数和析出固相数目无关。蒸发结晶过程自由度为1,物理意义是在不引起新相产生旧相不消失的前提下过程连续地变化(水分的连续蒸失),几何意义是指在适当的坐标系中的一条线,一条直线或曲线;0的物理意义是指新相的产生和旧相的即将消失或者前一段的连续变化过渡到后一段的连续变化,几何意义是指一点,两蒸发阶段的交点或者过程的终点。简单地讨论了过程的单向性和连续件。     

珠江口盆地西江主洼油气差异分布机制

西江主洼是珠江口盆地一个低勘探程度洼陷，油气差异聚集特征明显，但其机理不清。本文在断陷盆地油气成藏理论指导下，利用研究区地质、地球物理和地球化学等资料，开展了烃源岩、断盖组合、储集体系和运聚模式等方面的研究。结果表明，裂陷期沉降、沉积中心有序迁移，造成主力烃源岩自东向西由文四段迁移至文三段、文一+二段，控制东、西部油气差异分布。恩平组区域性泥岩发育且晚期断裂缺乏，导致油气纵向上更易聚集于下构造层；区域泥岩减薄尖灭或晚期断裂切开盖层的区域，上构造层有一定油气分布。"源-汇"类型从宏观上控制储层优劣，影响下构造层油气富集程度；"仓储"运移是控制上构造层油气规模聚集的主要模式。下构造层古近系应作为重点部署方向，东部围绕文四段、西部围绕文三段、文一+二段烃源岩构成的含油气系统展开；上构造层勘探需关注珠海组，在隆起周边寻找具备"仓储"运移模式的有利区带。     

钻孔灌注桩施工对邻近运营隧道变形影响研究

桩基施工将引起周边地层位移和邻近隧道结构变形,以南京龙津桥改建项目桩基工程为背景,通过现场实测,分析研究钻孔灌注桩施工全过程各工况条件对地层位移场和邻近既有隧道结构变形的影响规律。结果表明:由于钢套管的护壁作用,钻孔灌注桩施工过程中产生的最大地层位移和隧道结构变形较小,说明采用“钢套管边旋压边取土”、“群桩间跳施工”等工艺,对周围土体扰动影响程度较小,建议近隧桩基工程采用全套管灌注桩的施工方法,监测分析结果可为类似工程提供技术参考。     

30 m global impervious surface area dynamics and urban expansion pattern observed by Landsat satellites: From 1972 to 2019

Using more than three million Landsat satellite images, this research developed the first global impervious surface area (GISA) dataset from 1972 to 2019. Based on 120,777 independent and random reference sites from 270 cities all over the world, the omission error, commission error, and F-score of GISA are 5.16%, 0.82%, and 0.954, respectively. Compared to the existing global datasets, the merits of GISA include: (1) It provided the global ISA maps before the year of 1985, and showed the longest time span (1972–2019) and the highest accuracy (in terms of a large number of randomly selected and third-party validation sample sets); (2) it presented a new global ISA mapping method including a semi-automatic global sample collection, a locally adaptive classification strategy, and a spatio-temporal post-processing procedure; and (3) it extracted ISA from the whole global land area (not from an urban mask) and hence reduced the underestimation. Moreover, on the basis of GISA, the long time series global urban expansion pattern (GUEP) has been calculated for the first time, and the pattern of continents and representative countries were analyzed. The two new datasets (GISA and GUEP) produced in this study can contribute to further understanding on the human’s utilization and reformation to nature during the past half century, and can be freely download from http://irsip.whu.edu.cn/resources/dataweb.php.

     

利用主成分分析改善土地利用变化的遥感监测精度——以珠江三角洲城市用地扩张为例

近年来,珠江三角洲由于经济的快速发展,城市用地急剧增加,利用多时相的遥感图,可以定量地监测这种城市化的现象。但理,由一般的遥感动态监测方法所得的结果往往夸大变化的程度,以及获得一些不合理的结论．该文提出主成分分析的方法来改善遥感动态监测的精度。将该方法应用应用于珠江三角洲发展最快的东莞市,获得了较满意的结果。     

北京云蒙山地区花岗岩穹隆及伸展构造的探讨

北京云蒙山花岗岩为一中生代侵入的花岗岩穹隆,花岗岩穹隆的叶理普遍发育,叶理轨迹基本平行于穹隆的外部边界,并显示出从核部到边部逐渐增强,东南侧明显强于西北侧的特点。变形构造研究显示,花岗岩穹隆的边部及围岩中普遍存在不同层次及不同运动方向的伸展构造。东南侧以具河防口－水峪伸展型韧性剪切带为特征,剪切运动标志显示为从ＮＷ－ＳＥ的正剪切运动,有限应变分析估算其剪切位移量在１０ｋｍ以上,剪切带上部被河防口正     

Research on temporal and spatial distribution, evolutionary character and mechanism of crustal deformation field before and after the Tangshan earthquake

Ｒｅｓｅａｒｃｈｏｎｔｅｍｐｏｒａｌａｎｄｓｐａｔｉａｌｄｉｓｔｒｉｂｕ┐ｔｉｏｎ,ｅｖｏｌｕｔｉｏｎａｒｙｃｈａｒａｃｔｅｒａｎｄｍｅｃｈａｎｉｓｍｏｆｃｒｕｓｔａｌｄｅｆｏｒｍａｔｉｏｎｆｉｅｌｄｂｅｆｏｒｅａｎｄａｆｔｅｒｔｈｅＴａｎｇｓｈａｎ...     

参数的期望估计及其在形变分析中的应用

利用“参数的期望估计”能准确定位、定量粗差及参数的期望估计不受粗差影响折特殊性质,寻找地壳变形区域和不变形区域,以确定拟稳点,进而进行拟稳变换,分析地菜变,将是一种有效方法。     

Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.