南海海底地震仪异常数据的分析和处理

海底地震仪(Ocean Bottom Seismometer,OBS)数据处理至关重要,是获取深部地壳结构的基础与前提.2006年实施OBS2006-2测线时,有2台OBS(OBS03,OBS06)数据出现异常,无法使用.由于海上航次花费巨大,采集到的数据弥足珍贵.本文采用数据格式检查、邻近台站对比分析、重采样等方法,成功地对这2台OBS数据进行了解编处理,得到了这两个台站的综合地震记录剖面;利用上述方法对2011年实施的OBS973-3测线中的异常台站OBS03进行了分析处理,同样得到了OBS03台站的综合地震剖面;通过查看两次海上实验班报发现,OBS2006-2测线之OBS06与OBS973-3测线之OBS03内部Sedis编号相同,为同一台记录仪器,再一次验证上述处理方法正确可行;然后对OBS2006-2测线2个台站进行震相识别与走时拾取后,利用前人纵波速度模型开展了射线追踪与走时模拟.此次对异常OBS数据的重新处理工作,不仅为OBS探测提供了宝贵的数据处理经验,而且将提高OBS2006-2测线地壳结构的可靠性和约束性,具有重要的研究意义.     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

陇东黄土高原下垫面不均匀性指标的建立及其对大孔径闪烁仪(LAS)观测感热通量的影响

目前定量研究下垫面不均匀性对大孔径闪烁仪(LAS)、涡动相关仪(ECS)观测感热通量差异的影响还比较少.本文利用黄土高原庆阳观测站2012年6、7月典型晴天两主风向范围E-SE和SW-W的陆面过程数据,在建立了下垫面不均匀性定量化指标的基础上,分析了下垫面不均匀性与地表温度变率的关系及其对LAS、ECS观测感热通量差异的影响.结果表明:庆阳站下垫面不均匀性大小η和地表温度变率г具有很好的相关性,相关系数达到0.566以上,证明了新建立指标的合理性.下垫面不均匀性大小η和LAS、ECS观测的感热通量差异HLAS-HECS具有很好的一致性,相关系数达到0.634.下垫面不均匀性对LAS和ECS的感热通量差异影响显著,下垫面越不均匀LAS和ECS测量的感热通量差异越大.对两主风向分别进行分析,在E-SE风向范围η和HLAS-HECS的相关系数为0.430,HLAS和HEC拟合的线性趋势系数为1.279,在SW-W风向范围η和HLAS-HECS的相关系数为0.680,HLAS和HEC拟合的线性趋势系数为1.297.下垫面不均匀性的影响程度越大,LAS和ECS观测的感热通量差异越大.     

西沙群岛西科1井晚中新世-上新世生物礁沉积的磁性地层学初步结果

本文对西沙群岛西科1井钻孔岩心晚中新世-上新世生物礁沉积进行了详细的岩石磁学、磁性扫描与磁性地层学研究.结果显示,西科1井生物礁相沉积中的载磁矿物主要是磁铁矿.我们推测,这些磁铁矿的微小颗粒主要来自海水中含有的陆源物质,在生物生长过程中通过珊瑚体或其它寄生生物对海水的过滤与吸附作用保存在生物礁沉积中.进一步的磁性地层学研究及其与地磁极性年表的对比发现,在上新统莺歌海组、上中新统黄流组内部可获得多个年龄控制点,并对莺歌海组和黄流组的底界位置给出了初步制约.其中,莺歌海组记录了C2An.3n和C3n.4n,黄流组记录了C3An.2n至C5n.2n.虽然这一对比方案存在一定不确定性,但是在现阶段生物地层年代及其它年代学资料相对匮乏的情况下,我们认为本项研究的磁性地层学结果能为西沙群岛晚中新世以来的生物礁沉积提供更多可靠的年龄控制点,并为今后的区域地层对比提供磁性地层年代学依据.     

2.5D direct‐current resistivity forward modelling and inversion by finite‐element–infinite‐element coupled method

To reduce the numerical errors arising from the improper enforcement of the artificial boundary conditions on the distant surface that encloses the underground part of the subsurface, we present a finite‐element–infinite‐element coupled method to significantly reduce the computation time and memory cost in the 2.5D direct‐current resistivity inversion. We first present the boundary value problem of the secondary potential. Then, a new type of infinite element is analysed and applied to replace the conventionally used mixed boundary condition on the distant boundary. In the internal domain, a standard finite‐element method is used to derive the final system of linear equations. With a novel shape function for infinite elements at the subsurface boundary, the final system matrix is sparse, symmetric, and independent of source electrodes. Through lower upper decomposition, the multi‐pole potentials can be swiftly obtained by simple back‐substitutions. We embed the newly developed forward solution to the inversion procedure. To compute the sensitivity matrix, we adopt the efficient adjoint equation approach to further reduce the computation cost. Finally, several synthetic examples are tested to show the efficiency of inversion.     

基于球面QTM的大范围有源淹没算法研究

面对海平面上升后沿海区域大范围淹没模拟的应用需求,提出了一种基于球面QTM的大范围有源淹没算法,首先构建了研究区域的QTM三角格网集,并利用海陆分离的思想确定了初始淹没单元;然后基于QTM的12邻近搜索设计了面向多目标三角格网的淹没范围扩张方法,将初始淹没单元递归扩张得到最终的淹没范围;最后应用Visual C#语言及Microsoft DirectX三维图形接口设计开发了相应的试验原型系统,并利用天地图高分辨率影像和SRTM地形数据对试验区域进行了海平面上升后的大范围淹没模拟试验。结果表明,该算法具有较好的实用性,可以有效地对海平面上升后的大区域范围进行动态地淹没模拟,对辅助防灾减灾具有一定意义。     

火星快车HRSC线阵影像投影轨迹法近似核线重采样

火星快车高分辫率立体相机(high resolution stereo camera, HRSC)线阵推扫式影像具有分辫率高、全球覆盖性好等特点,对火星快车HRSL影像进行摄影测量处理可获取火星全球高分辫率地形数据,而影像匹配是自动获取地形数据的关键利用核线几何约束可以将影像匹配搜索范围由二维降至一维本文在构建火星快车HRSC线阵影像严格几何模型的基础上,提出了一种近似核线重采样方法首先,利用投影轨迹法分析 HRSC线阵影像核线的直线性与平行性;然后,将原始HRSC影像进行微分纠正生成水平影像,分析投影轨迹法得出的核线在水平影像上的几何关系并确定近似核线方向,再沿核线方向进行核线重排生成近似核线影像实验结果表明,本文算法生成的近似核线影像同名点上下视差优于1个像素,且在水平影像上进行近似核线重采样减少了由外方位元素变化引起的几何畸变,更有利于影像匹配。     

主动源与被动源地震数据插值及联合数据成像

本文提出了两种情况下主动源数据和被动源数据的插值方法,并研究了两种数据在偏移成像中的互补效果.基于互相关法被动源数据重构原理,本文提出了结合多域迭代去噪技术的重构方法.提出了两种时间域主动源和被动源数据的插值方法,分别是共炮点域能量匹配插值和共检波点域最小平方匹配插值.然后对获得的主动源和被动源联合地震数据进行叠前深度偏移成像.在被动源活跃度不是很高的地区进行被动源地震勘探时,少量的主动源地震数据可以有效控制和补充被动源数据的成像效果.在稀疏炮点的主动源勘探中,有效利用被动源的信息能够在成像中增加更多的细节信息,提高成像质量.     

CRInSAR与PSInSAR融合探测地表线性变形速率

提出利用CRInSAR与PSInSAR融合解算地表变形参数的算法。将CR点上获取的线性变形速率与高程改正结果作为PS基线网络的起算数据,依据最小二乘原理求解PS点上待求线性变形速率与高程改正值的最优解。这一算法可融合PSInSAR与CRInSAR两种算法的优势,起到很好的互补作用。实验显示,这一算法获取的线性变形速率精度可达±0.37 mm/a,高程改正值精度可达±0.5 m,证实该算法在实际地表变形监测中具备可行性。     

陆态网络腾冲台连续重力观测质量分析

对陆态网络腾冲台2012-07-01~2013-06-30连续重力观测资料进行处理,获得该台站观测精度及噪声水平。结果表明,腾冲台1~100s噪声水平在9×10^-8ms^-2左右,100s~100min噪声水平在0.3×10^-8ms^-2左右,白噪声水平为0.509×10^-8ms^-2,对O₁波、M₂波等主要潮波振幅因子测量精度达到SG-053超导重力仪的观测量级,其中M₂波振幅因子测量精度最高（为0.000 18）。腾冲台可为川滇地区强震和腾冲火山活动监测、震前重力变化及相关研究提供高精度、可靠的连续重力观测资料。