InSAR复影像配准方法与实践

利用合成孔径雷达干涉测量技术获取地表三维信息的研究，受到世界各国的广泛关注。高精度复影像配准是InSAR数据处理的关键步骤，其结果直接影响生成干涉相位的质量，进而影响获取数字高程模型的精度。本文采用相干系数法进行粗配准，采用基于相干系数拟合的方法进行精配准，可以实现复影像数据的高精度配准。经机载和星载数据的试验验证，该方法有效可行且运算效率高，具有实用性。     

基于速度分析的探地雷达阻抗介电常数反演

阻抗反演是利用波阻抗与介电常数关系开展地下介质参数估计的重要技术,在探地雷达以及叠后地震资料解释中具有广泛的应用。常规阻抗反演需要钻孔或测井曲线作为约束项,约束项信息直接影响阻抗反演的估计精度。在缺少钻孔数据的实际应用中,如何开展探地雷达阻抗反演是该方法研究的重要内容之一。基于上述问题,本文提出了基于速度分析的探地雷达阻抗反演方法。其基本思想是基于多偏移距雷达数据开展速度谱分析和Dix反演,以获得不同深度的速度信息作为阻抗反演的约束项;同时,采用K-means方法自动拾取速度谱信息,大大降低了常规人工拾取误差,提高了计算效率。通过典型随机土壤介质模型,验证了本文方法在无钻孔条件下仍然可以获得较好的介电常数估计结果,并测试噪声适应能力强。最后通过美国密歇根州Wurtsmith AFB,in Oscoda区域的探地雷达数据测试了本文提出方法在探地雷达实测数据参数估计中具有较好的应用效果。     

内蒙古宁城地热资源的构造制约与勘查标志

内蒙宁城地热区位于华北板块北缘,是我国北方重要的矿带和矿田。本次研究对宁城地区成热控热构造系统进行厘定和研究,用成生联系的思想综合分析了不同尺度的控热断裂带之间的联系,包括嫩江-八里罕断裂、黑里河断裂和大庙-锦山断裂,阐述了不同系统的主干构造体系及其复合对地热分布的控制。嫩江-八里罕断裂南段"山"字形构造的NW向断裂,向SE向发生近40°错移,偏转后呈"S"型弧形弯曲,在转折端西侧构造应力集中的深成隐伏侵入岩蚀变带中出现地热异常。同时,SN向的热水镇-大营子断褶带后期继承性改造使断层向SN两翼外扩,为深成岩浆侵入及热能上升提供了通道,可能是地热形成的主要构造因素。由于岩浆侵入活动加剧,造成地热区差异性上升运动,使壳源深部热源上移;地下水在深部循环加热后,在高海拔补给带的静水压力下沿断裂通道上升,形成地热流体。     

一种基于扩展卡尔曼滤波的多基线相位噪声抑制及展开方法

针对传统多基线干涉合成孔径雷达存在系统高度模糊数与缠绕相位解缠可靠性之间不可兼顾的矛盾,该文提出一种基于扩展卡尔曼滤波(EKF)的多基线SAR干涉相位展开方法。该方法利用全方位局部相位梯度联合估计技术直接从多幅复干涉图抽样谱中提取相位梯度信息,同时利用EKF强大的数据融合能力以信息融合的方式同时完成二维噪声抑制及相位展开,既有利于解决相位展开中的"坡度欠估计"问题,又避免了传统方法在相位展开之前须进行噪声滤波的不足。多基线模拟数据和单基线实测数据处理结果验证了方法的有效性。     

附不等式约束最小二乘平差的一种新算法

在大地测量数据处理中,很多情况下可根据先验知识建立合理的不等式约束,能够改善平差结果,提高精度。首先简要总结了附不等式约束平差的各种方法及存在的问题。根据有效约束和库恩塔克条件,提出了解决不等式约束平差的新算法,把不等式约束平差转化为等式约束平差问题,从而得到解的显示表达。最后用一数值算例证明了该算法的可行性。     

沉积记录对白垩纪中国大陆古气候模拟的约束

摘要：古气候模拟已经成为白垩纪古气候研究中不可替代的技术手段。但是目前中国白垩纪陆相古气候模拟的一些关键边界条件尚未完全确定。通过白垩纪沉积记录对中国大陆白垩纪古海陆边界和古地形做出约束。早白垩世,新疆西南部喀喇昆仑地区可能存在喀喇昆仑海湾,在东北松辽盆地东侧有乌苏里湾;晚白垩世,受全球海平面上升的影响, 特提斯海向北延伸,新疆西部地区存在塔西南海湾和喀喇昆仑海湾;而受燕山运动的影响,乌苏里湾向东退出松辽盆地,影响变小。在古地形方面,认为东北大兴安岭,西北天山、祁连山、昆仑山等山脉,西南云贵高原、川西高原、横断山,秦岭-大别山在白垩纪时均已形成并隆升到一定高度;指出“东部高原”和“沿岸山脉”的存在与否可能对白垩纪古气候状况产生深远的影响。此外, 还讨论了沉积记录与古气候模拟结果的对比、验证对中国大陆白垩纪古气候重建的重要作用。     

2D cooperative inversion of direct current resistivity and gravity data: A case study of uranium bearing target rock

Interpretation of a single geophysical data set is not sufficient to get complete subsurface information. Cooperative or joint inversion of geophysical data sets is the preferred method for most case studies. In the present study, we present the results of the cooperative inversion approach of direct current resistivity and gravity data. The algorithm uses fuzzy c-means clustering to determine the petrophysical relationship between density and resistivity to obtain the similarity. Synthetic data set has demonstrated that the cooperative inversion approach can produce more reliable and better resistivity and density models of the subsurface as compared to those obtained through individual inversions. To utilize the presented cooperative inversion algorithm, the number of geologic units (number of clusters) in the study region must be known a priori. As a field study, the cooperative inversion approach was used to identify the extension of uranium-bearing target rock around the Beldih open cast mine. We noted the inconsistencies in both resistivity and density models obtained from the individual inversions. However, the presented cooperative inversion approach was able to produce similar resistivity and density models while maintaining the same error level of the respective individual inversions. We have considered four geologic units in the presented cooperative inversion as a field case study. We have also compared our cooperative results with drilled borehole and found to be a reliable tool to differentiate between the target rocks (kaolinite and quartz–magnetite–apatite rocks) and the ultramafic rock (host rock quartzite/alkaline granite). However, this study is subject to certain limitations such as the inability to differentiate between closely spaced kaolinite and quartz–magnetite–apatite rocks.     

基于一种反余切紧凑约束的电磁波走时层析研究

电磁波走时层析可用于空洞和孤石等离散块体异常的探测中.在传统电磁波走时层析中,通常采用一个低阶差分泛函来稳定反演过程,而这类泛函都具有光滑约束的性质,在成像结果中难以对此类离散块体异常进行准确的解释.本文基于Tikhonov正则化思想,提出了一种反余切泛函,该泛函具有紧凑约束的性质且无需选择一个聚焦因子.结合弯曲胖射线层析理论和重加权正则化共轭梯度反演算法,对两个理论模型进行了成像效果研究.与低阶差分泛函和最小支撑泛函的成像结果相比,该方法能够更好地突出异常的边界,成像结果中的虚假异常也要更少.同时,还分析讨论了激发电磁波的中心频率以及数据噪音对成像结果的影响.此外,针对一个地下连续墙墙体内部缺陷的现场模型,利用该方法取得了理想的成像结果,良好地揭示了缺陷的形态和位置.     

Three‐dimensional potential field data inversion with L0 quasinorm sparse constraints

The quantitative explanation of the potential field data of three‐dimensional geological structures remains one of the most challenging issues in modern geophysical inversion. Obtaining a stable solution that can simultaneously resolve complicated geological structures is a critical inverse problem in the geophysics field. I have developed a new method for determining a three‐dimensional petrophysical property distribution, which produces a corresponding potential field anomaly. In contrast with the tradition inverse algorithm, my inversion method proposes a new model norm, which incorporates two important weighting functions. One is the L0 quasi norm (enforcing sparse constraints), and the other is depth‐weighting that counteracts the influence of source depth on the resulting potential field data of the solution. Sparseness constraints are imposed by using the L0 quasinorm on model parameters. To solve the representation problem, an L0 quasinorm minimisation model with different smooth approximations is proposed. Hence, the data space (N) method, which is much smaller than model space (M), combined with the gradient‐projected method, and the model space, combined with the modified Newton method for L0 quasinorm sparse constraints, leads to a computationally efficient method by using an N × N system versus an M × M one because N ? M. Tests on synthetic data and real datasets demonstrate the stability and validity of the L0 quasinorm spare norms inversion method. With the aim of obtaining the blocky results, the inversion method with the L0 quasinorm sparse constraints method performs better than the traditional L2 norm (standard Tikhonov regularisation). It can obtain the focus and sparse results easily. Then, the Bouguer anomaly survey data of the salt dome, offshore Louisiana, is considered as a real case study. The real inversion result shows that the inclusion the L0 quasinorm sparse constraints leads to a simpler and better resolved solution, and the density distribution is obtained in this area to reveal its geological structure. These results confirm the validity of the L0 quasinorm sparse constraints method and indicate its application for other potential field data inversions and the exploration of geological structures.     

顾及多特征约束的建筑物组合化简方法

建筑物是大比例尺地图中的重要要素,其自动化简是制图综合领域的重要研究内容。针对当前化简方法仅侧重于某一方面的问题,本文提出一种顾及多特征约束的建筑物组合化简方法。首先,对建筑物进行冗余点删除、尖角去除等预处理操作;然后根据化简阈值筛选出需要化简的短边,并对其局部结构进行自动分类识别;然后,充分考虑建筑物的位置、面积、方向和形状等特征约束的影响,针对不同的结构类型采用不同的化简操作,在化简的过程中采用阈值逐渐增大的渐进式确定方法,直到满足化简要求为止;最后,从多个特征角度对化简结果进行评估,对不合理的化简实体进行调整修改。以OpenStreetMap数据为例进行实验,并与其它化简方法进行对比,实验结果表明,本文所提的方法可以有效保持建筑物的面积、形状、直角等基本特征,可以实现建筑物的自动连续化简,具备一定的通用性。