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的物理意义是指新相的产生和旧相的即将消失或者前一段的连续变化过渡到后一段的连续变化,几何意义是指一点,两蒸发阶段的交点或者过程的终点。简单地讨论了过程的单向性和连续件。     

探地雷达技术在铁路隧道监测中的应用

周边工程施工威胁着铁路隧道运营安全,为了实时了解铁路隧道结构变化情况,需要进行隧道监测.探地雷达具有较好的探测和显示功能,常用于建筑结构物体的无损探测.为了弥补测绘技术不能探测隧道内部结构体的不足,将探地雷达技术应用于隧道结构监测.着重介绍了隧道探测的技术特点和探测参数设置,根据隧道介质反射图像特征,确定隧道内部结构信息.结果证明了隧道探测技术的有效性,研究内容可为相关工程技术应用提供参考.     

祁连山南部地区三叠纪弧后前陆盆地——基于沉积填充与大地构造的证据

本文通过对祁连山南部地区三叠纪沉积充填特征的研究,建立了盆地形成过程中沉积充填与区域构造演化的关联。依据沉积充填特征、构造接触关系及区域地质背景,认为祁连山南部地区三叠纪盆地具有弧后前陆盆地的构造属性;综合分析年代学、沉积物性和盆地内部填充单元的差异划分出与柴达木岛弧带和祁连山南部弧后前陆盆地相关的区域大地构造单元。认为早三叠世-中三叠世中期,南部遭受俯冲作用持续向北挤压推覆导致盆地挠曲沉降;中三叠世中期-晚三叠世中期发生弧(东昆仑-柴北缘)-陆(巴颜喀拉地块)碰撞,祁连山南部盆地挠曲沉降增大,其间形成区域性断裂和不整合等盆地记录的地质事件;晚三叠世晚期祁连山南部盆地发生过补偿陆相填充,弧后前陆盆地消亡。     

内蒙古中西部宝音图群石榴角闪岩的地球化学、变质作用及年代学特征

宝音图群分布于内蒙古中西部狼山、图古日格和达茂旗一带，主要由云母片岩、石英片岩、石英岩、角闪岩和大理岩组成。在狼山格尔敖包沟和图古日格西南地区发育十分典型的石榴角闪岩，均以似层状或透镜状产自云母片岩中，但两地具有不同的矿物组合，前者主要以石榴石+角闪石+斜长石+含钛矿物（金红石和钛铁矿）为特征，后者主要以石榴石+角闪石+绿帘石+含钛矿物（金红石和钛铁矿）为特征。地球化学研究显示两个地区的石榴角闪岩原岩具有相似的化学组成，均为亚碱性玄武岩的拉斑系列。宝音图群的石榴角闪岩具有右倾的稀土配分模式，REE总量为83.31×10^-6~125.9×10^-6，（La/Yb）_N比值为2.17~6.48，δEu=0.87~0.98。Ta、Nb、Ti没有明显的负异常，其配分型式类似E-MORB特点。微量元素构造判别图解表明这些石榴角闪岩的原岩产于板内拉张环境中。本文通过模拟格尔敖包沟样品LS01和图古日格西南地区样品LS35的P-T视剖面图研究其变质作用及矿物演化过程，结果显示这两个样品均经历了早期进变质，峰期以及峰后近等温减压的顺时针型P-T轨迹。样品LS01的峰期温压条件为~11kbar/~735℃，模拟得到的峰期矿物组合可能为石榴石+角闪石+斜长石+黑云母+透辉石+金红石+熔体，P-M（H₂O）视剖面图显示，由于峰后演化过程中存在饱和流体渗透，导致峰期透辉石无法保留。样品LS35的峰期温度压力条件为~8kbar/675℃，其峰期组合为石榴石+角闪石+绿帘石+金红石。两个地区的石榴角闪岩的温度压力条件明显不同，推测宝音图群中呈透镜体或似层状产出的变质基性岩，存在着递增型的变质作用。锆石LA-ICP-MS U-Pb年代学的研究结果显示石榴角闪岩样品LS01的变质锆石年龄为394±8Ma。由此推测这些石榴角闪岩的原岩可能形成于华北北缘在新元古代发育的裂陷盆地，在泥盆纪中期，宝音图群作为华北克拉通的一部分，被卷入到与古亚洲洋闭合有关的造山过程，发育中压型递增变质作用。     

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

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

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

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

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.     

Characterizing heterogeneous permeable media with spatial statistics and tracer data using sequential simulated annealing

Characterizing heterogeneous permeable media using flow and transport data typically requires solution of an inverse problem. Such inverse problems are intensive computationally and may involve iterative procedures requiring many forward simulations of the flow and transport problem. Previous attempts have been limited mostly to flow data such as pressure transient (interference) tests using multiple observation wells. This paper discusses an approach to generating stochastic permeability fields conditioned to geologic data in the form of a vertical variogram derived from cores and logs as well as fluid flow and transport data, such as tracer concentration history, by sequential application of simulated annealing (SA). Thus, the method incorporates elements of geostatistics within the framework of inverse modeling. For tracer-transport calculations, we have used a semianalytic transit-time algorithm which is fast, accurate, and free of numerical dispersion. For steady velocity fields, we introduce a transit-time function which demonstrates the relative importance of data from different sources. The approach is illustrated by application to a set of spatial permeability measurements and tracer data from an experiment in the Antolini Sandstone, an eolian outcrop from northern Arizona. The results clearly reveal the importance of tracer data in reproducing the correlated features (channels) of the permeability field and the scale effects of heterogeneity.