纪念中国大陆科学钻探实施15周年、国际大陆科学钻探委员会成立20周年

大陆科学钻探是“入地”的重要手段,是“深入地球内部的望远镜”。中国大陆科学钻探事业开展15周年以来,取得重要进展,获得全球地学界的高度关注,特别是2001年实施的中国第一口大陆科学深钻 （5158m）,成果辉煌,影响巨大。继后,又开展了青海湖环境科学钻探、松辽盆地白垩纪科学钻探、柴达木盐湖环境资源科学钻探,汶川地震断裂带科学钻探以及中国大陆科钻资源集成计划,总共钻进约 35km,显示了中国科学钻探方兴未艾的景象。为纪念国际大陆科学钻探20周年（1996~2016）和中国大陆科学钻探实施15周年（2001~2016）,本文回顾中国大陆科学钻探实施15年来的艰辛和奋斗的历程,展望中国大陆科学钻探的未来。     

北淮阳构造带的属性及其对板块缝合线位置的启示

北淮阳构造带为一条沿大别造山带北麓展布的狭长构造带,是华北板块与扬子板块碰撞造山的产物,其物源组成、构造属性等研究内容是揭示南北板块在中国东部的大地构造单元格局以及时空演化的关键信息。结合构造变形、岩石学、岩石地球化学、锆石U-Pb年代学等方法对北淮阳前中生代岩石地层单元苏家河群/卢镇关群和信阳群/佛子岭群进行对比分析,就其构造属性形成这样的基本认识:北淮阳构造带是杨子板块向华北板块俯冲-碰撞的结合带,其物质组成具有双源性-主体部分来源于杨子板块,由代表了杨子板块北缘基底的晚元古代苏家河群(784±19Ma)/卢镇关群(756±12Ma)和杨子板块被动陆缘复理石沉积建造-中泥盆世至晚泥盆世早期的信阳群南湾组组成,后者相当于商丹缝合带以南的刘岭群;信阳群龟山组亲缘于华北板块并构成华北板块南缘的弧前楔形沉积体,是沿商丹缝合带分布的桃花铺弧前楔形沉积体的东延。佛子岭群总体可与信阳群相连,全岩成分显示,佛子岭群沉积构造环境既有活动大陆边缘,又有被动大陆边缘。由于中生代以来多期构造变形叠加(至少三期),佛子岭群岩石地层单元的叠置关系被改造或重置,使其缺少"构造解体"的证据,而无法在构造单元上细分。对于两大板块之间的缝合带位置,笔者认为:商丹缝合带进入北淮阳构造带后,应沿着北淮阳的北界,与龟山-梅山-六安-磨墩水库一线延伸的断裂带相连。     

辽东湾坳陷东部地区走滑双重构造的发育特征 及其石油地质意义

在将构造发育特征与走滑双重构造理论模型进行类比分析的基础上,对辽东湾坳陷东部地区的走滑双重构造的发育特征进行了研究,并结合油气勘探实践成果,就其石油地质意义进行了探讨,研究结果表明:在受控于走滑与拉伸并存的多动力源区域地质背景下,辽东湾坳陷新生代构造演化体现了断陷与拗陷、拉张与走滑的叠加效应,NNE走向的辽中1号断裂、辽中2号断裂及辽东断裂的发育不仅控制了辽东湾坳陷东部地区的凸凹格局,其走滑侧接关系也导致了研究区走滑双重构造的发育;辽东凸起是由于辽中2号断裂、辽东断裂走滑弯曲、右行左阶侧接而形成的挤压走滑双重构造,而辽中凹陷在早期拉张断陷基础上,经历了后期走滑作用的改造,其断裂体系的发育体现了辽中1号断裂与辽中2号断裂的右旋右行侧列所产生的拉张走滑双重构造效应;在走滑双重构造的增压带,断裂侧向封堵性强,有利于形成有效的断层圈闭,是油气藏发育的有利区带。     

奇异值分解法在病态问题中的应用

用截断奇异值分解法和修正奇异值分解法对大地测量病态问题进行了处理,并与最小二乘的结果进行了比较,最后将两者方法进行结合同样对病态方程进行了处理,得到了结合奇异值分解的解,并与截断奇异值、修正奇异值分解的解以及真值进行比较,发现结合奇异值分解的解即修正奇异值截断法比截断奇异值和修正奇异值的解更加靠近于真值,修正奇异值截断法相比于截断奇异值和修正奇异值法对于病态方程抗干扰能力更强,更具有实际意义。     

利用多光谱卫星遥感和深度学习方法进行青藏高原积雪判识

青藏高原积雪对全球气候变化十分重要,针对已有积雪遥感判识方法中普遍采用的可见光与红外光谱数据易受复杂地形与高海拔影响,导致青藏高原地区积雪判识精度较低的问题,提出了一种基于多光谱遥感与地理信息数据特征级融合的积雪遥感判识方法:以风云三号卫星可见光与红外多光谱遥感资料与多要素地理信息作为数据源,由地面实测雪深数据与现有积雪产品交叉筛选出样本标签,构建并训练基于层叠去噪自编码器(SDAE)的特征融合与分类网络,从而有效辨识青藏高原遥感图像中的云、积雪以及无雪地表。经地面实测雪深数据验证,该方法分类精度显著高于使用相同数据源的FY-3A/MULSS积雪产品,略高于国际主流积雪产品MOD10A1与MYD10A1,并且年均云覆盖率最低。试验结果表明该方法可有效地减少云层对积雪判识的干扰,提升分类精度。     

利用引力梯度数据计算径向梯度的优化方法

根据重力梯度观测各分量的方差及协方差信息,提出了利用GOCE梯度数据计算径向重力梯度的优化方法。首先给出了径向重力梯度的计算方法,并深入分析了误差传播规律,通过建立相应的条件极值问题,给出了计算径向重力梯度最优组合因子的方法;通过模拟数据验证了本文所提出的优化因子的优越性。实际数据计算表明:相对于传统方法,采用优化组合因子可使反演所得引力位模型的累积大地水准面精度在250阶时提高约2 cm。由于径向重力梯度不仅可以用于地球引力场模型的求解,也可直接应用于地球物理问题的讨论,因此本文所提出的优化方法也可对部分地球动力学问题的讨论提供方便。     

利用压密方程计算分析西安市地面沉降特征

利用西安市地下水位监测资料,基于水文地质三维结构模型,在不考虑粘土层滞后压缩变形的理想情况下,采用压密方程计算获得西安市抽取地下水可引起的理论地面沉降量及沉降分布特征,结果表明:计算所得沉降显著区位于西安市西南部的高新区及东南部的曲江新区,沉降量呈由北向南递增的特征,与InSAR监测结果整体趋势具有较好的一致性;地面沉降分布特征受到活动地裂缝影响,沉降曲线呈近NE向偏转展布。研究结果可为沉降灾害预防研究提供宏观的参考信息。     

一种快速二维 Delaunay 三角网点定位算法

在构建二维Delaunay三角网的逐点插入法中,定位待插点所在三角形的快慢是影响整个算法构网速度的关键因素。针对目前已有算法存在的搜索路径长、搜索路径求解计算量大等问题,结合三角形重心的几何性质,对点定位算法进行改进,避免求三角形重心和相交边的过程。实验结果表明,文中算法较目前其他点定位算法能够有效地缩短搜索路径,减少点定位的计算时间,提高Delaunay三角网构网过程中点定位的效率。     

改进的距离约束最小二乘模糊度搜索算法

针对短基线最小二乘模糊度搜索算法搜索效率低的缺点,该文提出了一种改善模糊度搜索空间以提高模糊度搜索效率的方法。最小二乘搜索算法通过基线长度范围确定模糊度搜索空间,完成模糊度的搜索和确认。距离约束的短基线模糊度搜索空间是一个空心椭球,各模糊度之间具有相关性,影响模糊度搜索效率。该文通过最小二乘去相关调整方法降低模糊度之间的相关性,减少搜索次数。实验证明,该算法的搜索效率提高了25%以上,且随着基线长度的增加,搜索效率有更多的提升。     

结合投影变换和相位相关的海冰运动监测

为提高对渤海海冰旋转和平移运动的监测精度,提出一种基于投影变换的相位相关跟踪方法。选取连续8景静止水色成像仪(GOCI)图像序列,根据特征图像窗口的投影变换构造辅助函数,通过寻求函数最优解得到旋转角度集合,选择修正相关系数确定最佳旋转角度,同时根据相位信息实现海冰样本间亚像素级别的平移跟踪,消除传统相位相关法中因忽略图像相频特性所造成的匹配误差。实验结果表明,以手动测量的旋转角度为基准,该方法和传统相位相关法的旋转监测均方根误差的最大值/平均值分别为0.59/0.50与1.41/0.94,跟踪速度提高了50.6%;海冰平移运动的速度矢量与辽东湾的现场实测数据及历史资料记录数据基本一致。该方法对渤海海冰旋转和平移运动具有较好的跟踪效果。