欧拉梁动力反应初边值问题的微分求积解

动荷载作用下欧拉梁动响应的计算是一个初边值问题,通常很难得到解析解,传统数值方法一般是把空间和时间分别离散进行求解,计算相对复杂,效率也不高.针对分布动荷载作用下欧拉梁的振动偏微分方程,采用传统微分求积法,在空间和时间上同时进行离散;对于所有非0阶的初/边值条件,采用嵌入法在权系数计算中予以考虑.算例的数值结果与精确解的对比证明采用传统微分求积法处理此问题是可行的,而且是高效的.对于实际工程中的其他类似问题,该方法同样适用.     

航磁梯度数据解释新方法在迁安铁矿勘探中的应用

为更加清晰和准确地利用磁梯度数据获得地质体的分布特征,笔者采用解析信号总水平导数法和梯度欧拉反褶积技术组合实现磁梯度数据的解释,解析信号总水平导数的极值可有效地圈定出矿体的分布范围,梯度欧拉反褶积法可有效地避免背景异常的干扰,从而准确地给出磁性体的位置和深度信息。理论模型证明解析信号水平导数法和梯度欧拉反褶积法的组合可准确地获得磁性体的水平范围和深度。应用于河北迁安地区磁梯度数据的解释,圈定出5个呈条带状分布的较大型铁矿藏,埋深在200~450 m。     

中国地区水平地壳形变模型的建立

我国大陆板块是欧亚大陆板块的重要组成部分,其板块运动剧烈,地质构造复杂,是全球研究板块运动的重要地区。本文通过对各个块体分别建模和解算,得出结论:中国地区除东北块体相对于欧亚板块不动外,其他块体均有不同程度的形变运动,东部形变较小,西部形变较大;华北块体总体上以大约6mm/a的速率向东运动;华南块体以约5~8mm/a的速率向东南方向运动;新疆块体以约10~15mm/a的速率向北东方向运动;西藏块体也向北东方向运动,但运动速率从南到北逐步减少。     

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

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

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

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

对欧拉矢量附加不等式约束的RELSM模型及其应用

借鉴选权拟合思想,对区域地壳运动的整体旋转和线性应变模型中欧拉矢量部分附加不等式约束,将区域地壳运动的欧拉参数约束到由多年观测资料确定的平均欧拉参数值附近,建立附有不等式约束的整体旋转和线性应变模型。对中国大陆构造环境监测网络在环渤海区域的近几期GPS观测速度场数据进行拟合分析,结果表明,本文方法通过不等式约束可适当增强块体整体运动趋势的平滑性,获得对速度场更高的拟合精度,可进一步增强模型的物理解释能力。     

地震记录的时间域反褶积、仿真及在地震计方位角相对测量中的应用

介绍了一种用冲击不变法从速度响应平坦地震计二阶传递函数设计得到的递归滤波方法,比用双线性变换得到的递归滤波系数少,因此使用更方便.作者用此方法作了反褶积与仿真的数值实验和实际地震记录的短周期仿真试验,仿真结果与实际仪器记录符合较好;并将此法应用于地震计方位角的相关分析法相对测量,结果与不仿真的相关分析测量相比更接近寻北仪的绝对测量结果.     

基于非结构有限元的时间域海洋电磁三维反演

海洋电磁法是一种有效的钻前储层评价手段,可识别出海底构造是否储油气,从而减少干井率降低勘探成本.近年来,频率域三维海洋电磁反演得到了快速发展,但受到空气波的影响,其在浅海环境中无法取得很好的效果.为解决这一问题,本文研究时间域海洋电磁数据三维正反演.正演模拟和计算中,我们选择基于非结构网格的矢量有限元方法.该算法中空间离散采用非结构四面体网格,可很好地拟合复杂海底地形条件和地下复杂结构;而对于时间离散,我们采用无条件稳定的后推欧拉方法,以确保任意时间步长数值计算的稳定性.反演计算中,灵敏度信息采用伴随正演隐式进行计算;同时,依据时间域反演方法的特点采用L-BFGS方法计算模型修正量.通过利用合成数据反演结果证明了本文提出的三维时间域反演方法可用于复杂海底环境,特别是在浅海环境下的有效性.

     

Compression of the North Hemisphere derived from space geodesy

Introduction With the development of global tectonics and overall detections for global tectonics with multi-geophysical methods, ones can roundly study on the geological tectonics of sampling and magnetic stripe image, so as to summarize and interpret the geometrical and kinematical charac-teristics for the distribution of the ocean and the land, and spreading state of the global tectonics in a global scale. From a comprehensive view, the South and North hemispheres are clearly unsym-metrical…     

金融领域的随机建模与基于软件R的Monte Carlo模拟(4):随机微分方程模型

主要研究线性随机微分方程模型,为此定义Itô随机微积分,建立Itô公式.鉴于研究的重点是利用R软件进行数值模拟,所以详细讨论了过去10多年来随机微分方程数值解的研究.