均方误差意义下的正则化参数二次优化方法

Tikhonov正则化法是大地测量中应用最为广泛的病态问题解算方法之一。影响正则化法解算效果的重要因素是正则化参数,然而,最优正则化参数的确定一直是正则化解算的难题,如L曲线法确定的正则化参数具有稳定性好、可靠性高的优点,但存在过度平滑问题,导致正则化法对模型参数估值精度改善较小。本文从均方误差角度分析了正则化参数对模型参数估计质量的影响。基于奇异值分解技术,提出了由模型参数投影值分块计算均方误差的方法,避免了均方误差迭代计算,并基于均方误差最小准则给出了正则化参数优化方法,实现了对L曲线正则化参数的优化。数值模拟试验与PolInSAR植被高反演试验结果表明,正则化参数优化方法有效改善了正则化法解算效果,提高了模型参数估计精度。     

基于截断牛顿法的VTI介质声波多参数全波形反演

不同类别参数间的相互耦合使多参数地震全波形反演的非线性程度显著增加, 地震波速度与各向异性参数取值数量级的巨大差异也会使反演问题的性态变差.合理使用Hessian逆算子可以减弱这两类问题对反演的影响, 提高多参数反演的精度, 而截断牛顿法是一种可以比较准确地估计 Hessian 逆算子的优化方法.本文采用截断牛顿法在时间域进行了VTI介质的声波双参数同时反演的研究.不同模型的反演试验表明, 在VTI介质声波双参数同时反演中, 截断牛顿法比有限内存 BFGS (Limited-memory Broyden-Fletcher-Goldfarb-Shanno, L-BFGS)法能更准确地估计 Hessian 逆算子, 进而较好地平衡两类不同参数的同时更新, 得到了比较精确的反演结果.     

地震波频散效应与反Q滤波相位补偿

本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度.     

基于选权拟合法的电离层电子密度层析重构

附加约束的电离层层析算法是解决电离层电子密度反演中不适定问题的主要方法,为避免此类方法中约束权阵的选取不当对电子密度分布重构产生的不良影响,本文将选权拟合法应用到电离层层析成像技术中.该方法特别设计了依据电子密度空间分布特性构造参数权矩阵的方案.新方法有明确的物理意义,挖掘了隐含的信息量,为解决电离层电子密度反演中由于观测数据的不足等因素引起的不适定问题提供了一种新途径,可以得到符合客观实际的结果.数值模拟实验和实测数据的反演结果证实了该算法有效性、可靠性和优越性.     

主体功能区建设面临的问题及调整思路

主体功能区建设是推进我国区域协调发展的重大战略举措,但要从构想走向实际操作却面临一系列问题:一是自身需要完善和发展,在划分标准、划分层次、划分的基本单元等方面需要进一步探讨;二是面临行政区利益的冲突;三是面临与各种空间规划或区划的衔接和协调;四是在实施保障方面,面临差别化区域政策的制定、实施主体的明确和生态补偿机制的完善等问题.因此,主体功能区需要从以下5个方面进行建设:①选择合理的行政区单元;②逐步推进主体功能区规划的编制;③平衡主体功能区与行政区的关系;④建立主体功能区规划的核心地位;⑤完善主体功能区划的生态补偿机制.     

大地问题的截面椭圆弧长新解

分析了以前大地问题中截面椭圆弧长计算方法的缺点,在运用Mathematic的基础上,提出了在这一问题上的计算方便、过程简单的新方法,并用该方法进行了一个实例计算,效果良好。     

应用改进并行遗传算法反求水文地质参数研究

利用单纯形法局部搜索速度快和遗传算法全局寻优的特点,同时为了克服两种方法各自的弊病,提出采用混合单纯形技术的并行遗传算法(Hybrid Simplex Parallel Genetic Algorithm,SPGA)进行水文地质逆问题的求解。详细阐述了SPGA算法的具体操作和实现,并将该方法应用于水源地的地下水模拟反演中。计算结果表明,SPGA算法在水文地质参数寻优计算中具有比较好的可靠性和计算效率。     

不同围压下岩石中泊松比的各向异性

由弹性波纵横波速比计算得到的泊松比,在利用地震波反演地下结构和物质组成时,可以提供比单纯地利用纵波或者横波波速更强的约束. 为了研究不同岩石中泊松比及其各向异性随围压的变化,选取19块典型的岩石样品,在不同围压下沿相互正交的3个方向同时测量纵、横波速度,并通过计算得到了泊松比及其各向异性. 结果表明,在大部分样品中泊松比各向异性与纵波速度各向异性有一定的相关性,但在有的样品中泊松比各向异性与纵波速度各向异性没有明显的关系,甚至表现出完全相反的变化规律. 由于泊松比也表现出很强的各向异性,在排除微裂隙的影响后有的仍可以达到30髎以上,在利用实际地震波传播数据通过泊松比反演地下物质结构时必须考虑其各向异性的影响.     

柔度法纤维模型在方钢管混凝土柱滞回仿真中的应用

本文编制了基于二阶柔度法纤维模型梁柱单元的矩形钢管混凝土框架结构非线性分析程序,采用虚拟阶步法和改进的滞回路径追踪方法消除了其中截面切线刚度矩阵奇异时无法求逆的隐患,提高了算法的稳定性.运用此程序对矩形钢管混凝土柱进行滞回仿真,并与试验结果对比,结果表明:柔度法纤维模型梁柱单元可以有效解决框架结构的几何非线性问题和压弯耦合的材料非线性问题,使用较少的单元数量就可以得到较准确的计算结果;本文所采用的改进措施切实有效,保证了程序的稳定性.     

用于速度图象重建的层析成象法

本文提出了一种用于速度图象重建的层析成象法。特别注意到地震学上的ST与医学上的CT的不尽相同之处,并处理了由此引起的困难。与速度图象重建有关的正问题的分析表明,不同频率范围的资料其分辨能力是不同的,因此,区域地震与远震资料的解释应予注意。本文有关速度图象重建的讨论不仅涉及到与成象方法有关的反演理论和方法,还给出了进行数值计算的内存量和运算量的分析。关于重建图象的质量评价,还给出了反演解的可靠性分析,为分辨特性的描述提供了一种既简便又形象直观的方案。 本文的方法与现行主要方法的比较表明:1.本文的方法允许存在速度间断面,有利于揭示不同构造区域的地壳厚度的差异;2.在常用的方法中,作为正问题假定的常速度块同反演结果解释时非常速度块之间不协调。由于我们把给定网格点的速度值的插值函数作为速度的空间函数,这一困难已被解决;3.相对于Aki等的一般方法,本文的内存量节省约一个量级、运算量减少约一半;4.相对于ART类方法本身不能给出成象结果的可靠性分析,本文的方法则具有显著的优点。 作者在分析了天然地震资料的状况之后相信,对充分利用我国地震台网的现有观测资料,本文的方法是有效的。