GPS地壳形变监测网应变解算方法与程序设计

分析了在平面直角坐标系和大地坐标系内解算应变的方法及其特点,并以Matlab为开发平台,研制出一套在平面直角坐标系和大地坐标系中解算应变的程序GPSstrain。利用GPS监测网测站位移(速率),该程序能够解算剪切应变(γ1、γ2)、面应变、块体旋转、主应变及其方向、最大剪切应变及其方向,同时也可解算出这些应变分量的误差。另外该程序还具有高斯投影解算功能,根据各测站坐标能够自动形成所有的Delaunay三角形,并以图形方式显示。最后以实际算例给出了该解算程序的部分结果,并进行了初步的分析和讨论。     

高光谱遥感影像波段选择算法评价方法研究

波段选择在高光谱影像数据降维处理中尤为重要,其算法较多,但针对波段选择算法评价方法的研究却很少。该文提出一种波段选择算法的评价方法。引入信息量、类别可分性和相关性3个指标,对所选波段子集进行度量,通过比较不同波段选择算法所选波段子集的度量结果,实现对波段选择算法的评价。在试验中,将该评价方法应用于4种常用的波段选择算法,分析了评价结果并验证了该评价方法的可行性。     

流动单元划分新方案及其在临南油田的应用

把沉积学与储层物性相结合, 从流动单元体系的角度出发, 提出了流动单元划分的新方案.在流动单元体系内部划分出流动单元、亚流动单元和渗流区3个不同层次.在储层精细小层对比的基础上, 首先根据区域内连续分布的隔层把储层分成几个独立的流体压力系统, 即流动单元; 然后再根据不连续分布的隔层, 把一个流动单元进一步分成若干个亚流动单元; 最后根据储层物性的差别把流动单元/亚流动单元划分成不同的渗流区.按照这个思路, 选取临南油田的典型高产区块———以三角洲前缘亚相沉积为主的夏52块砂三中段三砂组, 进行了流动单元、亚流动单元和渗流区的划分, 共划分出7个流动单元、7个亚流动单元和63个渗流区, 这样划分出来的流动单元体系同时包括了油藏整体与局部细节的特征, 为油藏开发提供了详细的地质依据, 也在实际应用中取得了良好的效果.      

基于模糊预测对象的证据权方法及其在土地沙漠化评价中的应用

为了使在地理信息系统中被广泛用于点事件预测的证据权方法能对面事件进行评价和预测, 提出了一种新的基于模糊训练层的证据权方法.它是一种更广泛的证据权方法, 与普通证据权方法所不同的是, 它的训练层是模糊集合, 其取值是它的隶属度.通过适当的变换也可以把点训练层转换为模糊集合.因此, 该方法可以对面事件、点事件和线事件进行评价和预测.该方法可以处理训练层和证据层均为模糊集合的情况, 被称为双重模糊证据权方法.作为该方法的一个应用实例, 本文介绍毛乌素沙漠边缘的晋陕蒙地区土地沙漠化评价的应用实例.      

南海西部夏季冷涡的季节变化特征

基于美国海军的空间分辨率为0.5°×0.5°月平均的GDEM(Generalized Digital Environmental Model )三维温盐资料,采用P Vector方法,计算了南海西部夏季冷涡及附近海域的三维环流,分析了此冷涡的水平结构和垂向结构及其季节变化规律。结果表明南海西部夏季冷涡所占据的海域为一上升流区,对应的是低温区和气旋式环流,其温度场具有明显的垂向层化结构,冷涡随深度增加而减弱。冷涡中心在上层靠近越南沿岸,在温跃层以下有离岸的趋势。此冷涡具有明显的季节性变化：在5、6月份冷涡开始形成,其范围达到110.75°E,最大切向流速为8 cm/s;在7、8月份冷涡达到强盛,其范围向东扩展至112.50°E,冷涡中心也向东移至(110.50°E,13.25°N),最大切向流速增加为30 cm/s;在9、10月份冷涡开始衰减。     

不同尺度地形的SH波频率域响应特征研究

地形是影响实际地震记录的一个重要因素,理论计算已证明地形影响地震波的走时、能量、震相、波形等特征,同时地形也是产生尾波信号的重要原因.复杂且具有多尺度成分的地形对地震波传播的影响不仅与地震信号的频率有关,与其自身尺度成分也表现出强烈的依赖关系. 本文利用局域离散波数法模拟计算了不同尺度地形SH波的频率响应,得到地形尺度与频率之间似共振的响应关系. 这一关系可以更好地解释复杂地形对不同频率地震波传播的影响方式,同时可以指导数值模拟模型的构建,在关注的频率范围内合理取舍构造尺度成分,达到以最小的计算量得到最大的计算效益.     

改进的等效半空间法及瑞雷波频散曲线反演

层状模型的理论频散曲线表明,层状模型表面瑞雷波同一频率对应的相速度具有多值特点,这就是所谓的瑞雷波多阶性. 面波的多阶性导致了实测频散曲线的复杂性,当地层中存在软弱夹层时这一问题尤为突出,从而给实测频散曲线的解释带来了困难. 本文从工程实用的角度出发,根据等效半空间理论,提出了一种计算理论频散曲线的新算法——改进的等效半空间法,避开了面波多阶性这一复杂问题. 由此方法计算得到的面波相速度并不同于某个具体阶数的面波,而是对多阶面波的综合反映. 按照改进的等效半空间法编制了拟合反演程序,并在工程中进行了运用,取得较好的效果.     

非均匀介质热蠕变流动的数值求解

针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.     

无单元法用于地震波波动方程模拟与成像

波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.     

直接消除电磁耦合的斩波去耦方法

应用激发极化法时经常遇到电磁耦合的干扰.基于耦合效应和激电效应在双频波测量波形上具有不同的表现形式,提出一种在野外采集激电数据时,直接剔除电磁耦合而又基本上不损失激电信息的“斩波”去耦方法.本文从“斩波去耦”的方法原理、斩波对正常双频波产生的畸变、以及存在EM和IP效应时的双频波斩波等方面进行了详细的推导和分析.数值计算结果表明,该方法去耦效果明显,可以在现场条件下直接去耦,值得在激发极化法中推广应用