基于SVD-TLS算法的非预测欧拉反褶积

针对地下有多个异常源时,单一预测构造指数难于表征多个异常源。采用非预测欧拉反褶积以避免可能错误确定构造指数使得欧拉解过度发散的问题;同时针对欧拉反褶积超定方程组的条件数很大,致使欧拉反褶积解集中良解占优率低等解的非唯一性和解的不稳定性等局限性,采用奇异值分解总体最小二乘法（SVD-TLS算法）,以降低由于奇异值分析不当造成计算欧拉解非唯一性和解的不稳定性的问题,并利用SVD-TLS的截断误差构造阈值函数对解集进行过滤。数值结果表明了算法的有效性和可靠性。      

Calculation of seismic transmission speed of HTI medium

During the seismic wave propagation process, as for the anisotropic lower medium, the speed is a function of the propagating direction. This article focuses on solving the problem how to get the transmittance angle and speed, knowing the upper seismic wave propagation velocity and the angle of incidence conditions. The main theories used Snell law, Christoffel equation and speed characteristics. Taking the HTI medium as an example, the authors give the detailed solving process and draw the space velocity characteristic curve.     

塔里木盆地航磁场分析与磁源体结构

本文研究塔里木盆地区域磁异常图的反演及磁源体结构.由于众多异常的叠加和反演固有的多解性,区域磁异常图的准确解释是非常困难的.三维欧拉反褶积是一种确定地质体位置和埋藏深度的自动定量反演方法,比较适用于计算区域磁异常源的埋藏深度.由于大型克拉通沉积盆地地层具有上新下老的规律性,将磁异常源分解为三个深度层次,圈定它们各自的分布区域,便可将它们与形成的地质作用及时代联系起来,为准确解释区域磁异常图提供可靠的依据.本文应用三维欧拉反褶积反演方法,计算出的塔里木盆地深度为2~5 km、5~10 km、10~20 km三个等级的磁异常源,它们与形成的地质作用及时代分别为: 中生代构造运动,海西期玄武岩侵位和太古代结晶基底的变质作用;圈定了它们各自的分布区域.     

旅游网站信息流距离衰减的集中度研究

通过百度、Google、Alexa、Cnzz等多种网络搜索引擎和网站访问量统计工具获取了24个旅游网站分省访问量资料,运用Origin软件评估了旅游网站信息流距离衰减的曲线模型,在此基础上深入研究了旅游网站信息流距离衰减的集中度问题。研究经过3个阶段,逐级递进分别应用Zipf定律、地理集中度指数和指数模型3种方法,描述了空间分布的集中性,揭示了本地集中性与经济集中性特征,讨论了本地集中性与指数模型符合性的关系。研究发现:①旅游网站信息流的位序—规模分布符合Zipf定律,以单分形特征为主,其信息流的规模结构主体呈Pareto分布模式,具有明显的空间分布集中性,该集中性表现为随分维值大小而相应变化。②各旅游网站信息流距离衰减的区位商值大于1的省份主要为网站所在省或经济发达省份,距离衰减具有明显的本地集中性和经济集中性特征;各旅游网站信息流距离衰减的空间洛伦茨曲线呈内凹型,且多数网站基尼系数大于0.5,距离衰减的集中度较高、空间分布不均衡。③旅游网站信息流距离衰减的本地集中性受旅游网站性质影响,且其本地集中性越强,指数模型拟合优度指数越高,指数模型与网站数据点的匹配效果越好,各省份访问量越接近标准曲线,与拟合曲线的吻合度越高。④旅游网站信息流距离衰减的集中性特征,为旅游目的地确定目标市场和旅游网站的建设与营销提供了理论支持。     

川滇地区次级地块及其主要断裂带现今活动研究

以位于川滇区域的陆态网和四川省CORS站速度场为约束,结合研究区已知的断裂分布情况,通过使用欧拉矢量计算方法和反欧拉矢量计算方法求出每个次级块体的欧拉矢量,进而对川滇区域次级块体作进一步的调整。同时对该区域进行1°×1°的网格插值,并利用各次级地块的欧拉矢量反演出每个次级块体的长期运动方向及其速度;对相邻块体的边界处,结合各自块体的欧拉矢量求出同一点在2个块体的速度并作差,进而研究分析各主要断裂带的现今活动速率。     

用GPS位移场进行中小区域变形分析方法探讨

因局部区域相对全球框架不可避免地存在整体刚性平动分量,常影响对局部范围内地壳相对变形量级与特征的判读.对分离局部整体刚性平动分量的方法进行了探讨,通过实例分析认为,对京津地区28个GPS连续观测站给出的位移速率场,整体平移法和欧拉矢量法都能得到较好的效果,但相比之下,对于类似的中小区域,整体平移法更加简单明了.     

我国石油资源流流量位序-规模分布特征变化

以我国石油(特指原油)资源流动为研究对象,以我国省级行政单元作为石油资源空间流动的节点,从1985～2007年中选取5个特征时点.借助齐夫定律,研究了我国石油资源输出量和输入量的位序-规模分布特征.结果表明:(1)各特征时点石油输出量和输入量的位序-规模分布均满足齐夫法则.(2)从无标度区范围来看.石油输出量的无标度区范围变化不大,而输入量的无标度区范围呈现出先缩小后又逐年扩大的趋势.(3)从空间分形结构来看.石油输出量的无标度区没有分段,始终呈现为单分形,而输入量则经历了由双分形向单分形的演化.(4)从空间分布形态来看,各特征时点石油输出量均呈帕累托分布模式,而输入量的空间分布形态变动较大,经历了对数正态模式-帕累托模式-对数正态模式-帕累托模式的演变.我国石油输出量和输入量位序-规模分布的不同步变化,与大油田分布、石化产业布局以及石油流通体制等有着密切联系.     

考虑流固耦合的桥墩地震反应方法

跨海桥梁地震反应中,作用于桥墩上的动水压力具有明显的流固耦合特征。依据张量理论,推导时变区域散度变换关系及微分形式的几何守恒律;基于任意拉格朗日一欧拉描述,从采用欧拉描述的流体运动Navier—Stokes方程出发,推导时变区域的流体运动控制方程;给出流固耦合问题中的结构计算模型、耦合面接触条件、耦合场计算方法以及流场网格运动控制方法。以某跨海大桥为例说明桥墩地震反应方法,重点突出地震动输入、流场初始条件模拟等问题。计算结果表明：流固耦合理论能够模拟桥梁墩台地震反应中的流场和结构特性;流场初始条件的正确模拟可保证计算稳定性,并减少运算量;横向地震动激励下,桥墩基底剪力较大,纵向地震动激励下,结构运动剧烈;流固耦合系统中的线弹性结构在地震反应中具有明显的非线性特征。     

基于耦合欧拉-拉格朗日法充灌袋挤淤过程的有限元研究

在含水率大、承载能力低的软土地基上,利用充灌砂袋进行挤淤置换的方法工艺简单、施工方便、处理效果明显,已经广泛地应用在建造围埝和修筑道路中。充灌袋下沉到一定深度后,其重量与地基承载力相等,从而达到极限平衡状态。充灌袋挤淤所达到的深度与其高度存在一定的对应关系,在工程中确定此关系十分重要。通过3组不同宽度的挤淤模型试验,得到充灌袋下沉量与充灌袋高度的曲线关系、淤泥的位移特性及隆起量。利用ABAQUS中耦合欧拉-拉格朗日法（CEL法）,将充灌袋和淤泥分别用拉格朗日体和欧拉体描述,模拟充灌袋挤淤下沉的过程,将模拟得到的充灌袋下沉深度和高度的关系、淤泥的位移特点与试验结果以及理论公式进行对比。结果表明,CEL法模拟得到的结果与试验结果、理论解具有很好的吻合性,说明利用CEL法可以较好地模拟充灌袋的挤淤问题。     

求解模型对欧拉矢量的影响

利用GNSS水平运动场求解欧拉矢量时,待估参数的个数会对欧拉矢量结果产生影响,针对该问题本文分别从理论推导和具体算例角度进行了分析.首先从理论上推导出不同的求解模型得到的欧拉矢量的差异;然后以2004 ～2007年中国大陆GNSS水平运动场为基础,选用两种常用的求解模型(块体整体旋转模型和块体的整体旋转与均匀应变模型),讨论了求解模型对欧拉矢量及后续研究的影响.结果表明,两种模型下块体整体旋转的差异最大可达2.60mm·a-1,是不能忽略的.因此认为选用不同的块体运动模型会得到不同的地壳水平运动图像,在地壳水平运动分析中对此需加以重视.