论地球重力场奇积分计算

定义斜率Ｙ／Ｘ作为新的积分变量,将地球重力场中与距离倒数１／ｒ有关的奇异积分转化为非奇异。斜率变量使用前后的对比表明,在新变 量的形式下,数值积分收敛很快。     

基于三轴加速度的波浪测量技术研究

以浮标为载体,以TMS320C6416型DSP处理器、三轴加速度传感器、电子罗盘为硬件电路核心,利用加速度频域二次积分的波高计算方法以及三轴加速度倾角补偿电子罗盘的波向测量方式,进行海洋波浪要素测量技术研究,获取了某一海域不同时段的波浪特征值。     

HYSPLIT模式轨迹计算误差分析

基于NCEP 1°×1°再分析资料GDAS和NCEP2.5°×2.5°再分析资料GBL,利用美国NOAA 空气资源实验室ARL研发的一种用于计算和分析大气污染物输送、扩散轨迹的混合单粒子拉格朗日综合轨迹模式HYSPLIT,对2010年3月12日12时和田站距地100 m高度空气质点的3 d前向轨迹进行了模拟,分析了轨迹计算中的积分误差和分辨率误差.此次轨迹模拟试验中,积分误差对轨迹计算误差的贡献很小.随着积分时间的延长,积分误差略有增大.积分误差并未随着驱动数据分辨率的提高而减小.分辨率误差在积分的各个时刻并不相同,它与地形高度与天气系统有关.在轨迹模式中,轨迹计算基于风速在时间和空间上的线性插值,轨迹模拟的不确定性也与插值有关.使用不同分辨率的驱动数据对轨迹计算结果影响显著,分辨率误差对轨迹计算误差的贡献远大于积分过程中截断近似等带来的积分误差.     

一种改进的TSI3563积分浊度误差校正方法

TSI3563型积分式浊度计是一种性能出色的气溶胶散射系数观测仪器,然而由于仪器设计所固有的限制,TSI3563型浊度计观测结果包含有角度截断和非朗伯体光源两项系统性误差,会使观测结果较真值偏小10%左右。因此,需要对TSI3563型浊度计的观测结果进行校正才能得到较为精确的散射系数观测值。该研究利用2009年华北平原HaChi气溶胶外场观测数据测试了现有校正方法,结果显示,传统的校正方法在我国华北平原这样的高气溶胶污染地区并不适用。为此,提出一种改进的校正方法,利用同时观测的PM₁和PM₁₀数据,在校正方案中加入超微米粒子体积比这一参量,对于不同体积比采用不同的校正函数。利用实际观测数据检验后发现,改进方法的校正效果相对于传统方法有很大改善。     

双三次数值模式与理想Lamb波、龙卷扰动对比模拟分析

通过"双三次曲面拟合—时间步长积分—双三次曲面拟合……",实现成为一种双三次曲面拟合新算法数值模式。它对模式大气作非线性(空间二阶可导)描述,是用数值分析(三次样条/双三次曲面)与计算方法近似求解大气运动非线性偏微分方程组。用Z坐标系高分辨率双三次数值模式,采用原始大气运动方程组和采用欧拉时间积分方案,设置理想扰动,对大气声波Lamb波做数值模拟;采用大气不可压假定,除去声波,对大气重力波龙卷扰动做数值模拟。对比模拟发现:理想Lamb波以声速在水平方向传播,扰动能量频散快,仅能在原地持续数秒钟;理想龙卷扰动则在原地可持续约1 min,且波动内圈(外圈)自下至上始终处于风场辐合(辐散)状态,具冷心结构和气压驻波特征,可因水汽进入扰动与借助凝结潜热释放而得到发展。     

汶川地震映秀-北川地表破裂带虹口乡段精细填图、位移特征和地震构造分析

同震地表破裂带是地震破裂最直观的表现。地表破裂带的精细填图对理解地震破裂过程及发震机理有着重要意义。我们对5.12汶川地震映秀-北川同震破裂带南部的虹口乡区段开展了较为精细的填图工作,该区段地表破裂带的平面几何和同震位移特征沿走向均显示出公里尺度的横向变化和复杂性。根据同震破裂不连续性特征可把地表破裂分为4段,自南向北依次为深溪沟段、庙坝段、高原新村段和八角庙段。其中,深溪沟段和八角庙段同震垂直位错较大,达.5～6m;庙坝段同震垂直位移量普遍较低,1～2m或更低;  在高原新村附近,地表破裂带分为近平行的南北两支,北支以右旋走滑为主,南支以南东盘逆冲抬升为主,这与与映秀-北川同震破裂总体北西盘抬升的性质相反。本文对高原新村附近的反向陡坎进行了深入的探讨,并综合前人的研究结果,认为是破裂在从深部往上扩展的过程中,在浅表部位顺着发育在虹口东南部的大鱼洞-龙溪飞来峰后缘的具南东倾向的断面继承而成。与北川沙坝附近的反向陡坎一样,说明老断裂带构造几何特征对单次地震破裂的扩展有牵引作用。最后,讨论了八角庙有切割关系的两组断层擦痕的地震学和地质学意义,指出可能指示了映秀-北川断裂面上局部区段震前的初始应力水平较低。     

基于ALOS DEM的乌蒙山国家级自然保护区地貌特征研究

以乌蒙山国家级自然保护区为对象,融合野外地质地貌调查和ALOS DEM数字地形分析,定量与定性相结合阐述区内地形分异、地貌格局与发育演化等特征。结果表明：1）保护区基本地貌类型以大起伏中山（57.45%）为主,其次为大起伏亚高山（30.66%）和中起伏中山（11.17%）,海拔、起伏度、坡度、面积-高程积分值等定量地形指标的平均值从三江口片区至朝天马、海子坪片区逐渐降低,夷平剥蚀面级数逐渐减少;2）保护区及附近地区共发育7级夷平剥蚀面,其中Ⅰ～Ⅲ级为山顶面和残余夷平面,Ⅶ级为金沙江侵蚀阶地面,Ⅳ～Ⅵ级为剥蚀面,与区域3次构造旋回相对应;3）保护区面积-高程积分值在0.197～0.815之间,平均为0.503,超过80%的区域处于地貌演化的壮年阶段,表明区域地表侵蚀动力较强,在人类活动增强背景下存在较大的生态风险,未来的保护规划应考虑各个基本地貌类型区的地表动力差异。综合而言,独特、丰富的地貌结构使得保护区生物多样性及其所依附的生态环境要素得以保存,加强地貌多样性研究有助于从整体性和系统性的角度实现保护区资源环境管理与可持续利用。     

向下延拓在深部矿产勘探中的应用——以青海某矿区为例

勘测深部盲矿体时,矿体距离观测点较远,磁异常的相互叠加常常在地面形成具有一个磁异常中心的宽缓磁异常,分辨不清矿致异常和背景异常。通常用向上延拓、方向导数等处理的效果差强人意。笔者采用向下延拓方法处理了青海某矿区的宽缓磁异常,向下延拓深度为30倍点距,获得了精细的各矿体异常特征。通过与大量的见矿、未见矿钻孔比较,向下延拓场反映的磁异常位置与已见矿钻孔有较好的对应关系,同时指示出全部的未见矿钻孔都处于向下延拓场反映的磁异常外围。基于当前向下延拓方法的稳定计算问题已经取得突破,笔者认为针对一些大埋藏深度的宽缓磁异常采用向下延拓处理,区分叠加异常,圈定异常位置,可以获得事半功倍的效果,为钻孔布置提供佐证。同时对勘探精度也作了一些有意义的讨论。     

滩坝砂地震精细解释软件开发及实现

滩坝砂储层埋藏深,单个储层较薄(坝砂大多几米厚,而滩砂多在2m以下),横向变化大,地震反射弱,很难用常规软件进行精细描述与预测。这里在研究滩坝砂精细解释方法基础上,使用C#语言,在基于.NET的Visual Studio2010平台上,开发出一套用于滩坝砂精细解释的软件(Seis_BBS:Beach-Bar Sandstone)。Seis_BBS具有以下特点:1软件集滩坝砂特征分析、正演模拟、目标处理、精细描述于一体,具有很强的实用性;2软件中采用了三参数小波薄互层检测、古地貌与地震属性融合、多子波消除强屏蔽效应、基于HHT的目标处理等多项滩坝砂解释新方法、新技术,弥补了商业软件在这方面的不足;3模块由C#语言开发,具有很好的移植性,可扩展性高,可以随时将后续新的解释方法添加到软件中。     

SWAN雷达拼图VIL产品中鄂西南地物回波特征及其剔除方法

利用2011年6~8月的SWAN雷达拼图垂直累计液态含水量(VIL)产品,分析VIL产品中鄂西南地物回波特征,提出直接定位剔除的方法:首先建立鄂西南地物位置索引表,然后根据索引表在VIL拼图产品中剔除地物,最后对剔除地物后的空白区域通过局部递归插值来进行图像平滑处理。经过实际资料检验,在一定样本数的统计结果中取得了较好效果,能够更真实地反映VIL水平分布特征。