基于大炮初至的近地表层析反演方法研究及应用

详细研究了射线追踪、层析反演方法原理以及层析反演数值求解算法,提出了利用大炮初至直接建立初始地震模型,由小折射和微测井进行约束的层析反演近地表速度模型的处理流程,提出了基于多次回溯的高精度快速射线追踪算法和层析反演数值求解算法。采用双重网格技术提高正演模拟计算精度和层析反演的成像质量,另外,充分考虑多种质量控制方案,确保最终的反演速度模型精确。设计了10种不同地表地质条件的理论模型,进行模型试算,并在实际生产中应用,理论模型试算和实际资料反演结果证明了本文方法的正确性。     

建立涌水量曲线方程的最小残差法

本文在综述了目前稳定流涌水量经验方程建立方法的基础上,根据最小二乘法原理,提出了一种新的方法—最小残差法,并给出了计算实例。计算表明,最小残差法是建立稳定法涌水量经验方程的一种较好方法,计算过程在微机上易于实现。     

1996年上半年世界地震灾害综述

１９９６年上半年全世界共发生６级以上地震５８次，其中７级以上地震１０次，成灾的地震１０次，死亡人数总５１０人，伤约１．５万人，经济损失为８３亿－１００亿元人民币。其中我国境内灾害地震３次，死亡人数共计３４９人，伤约１．４５万人，经济损失约７８亿元人民币。我国因地震死亡的总人数约占全世界地震     

识别天然地震和人工爆破的判据选择

从快速识别事件性质的要求出发,分析了天然地震和人工爆破的发震时间,P波初动方向,P波、S波振幅比值,P波、S波最大振幅与尾波持续时间比值等判据,得到P波初动方向和P最大振幅与S最大振幅比值是识别爆破和地震的两个有效判据,为研制“识别天然地震和人工爆破的分类决策支持系统”提供了应用依据.     

一类信号恢复问题的投影迭代硬阈值算法

本文给出了基于L₀模求解该问题的非凸模型,借助于稀疏正则化方法来克服问题的不适定性。该模型利用紧小波框架对信号进行稀疏逼近,并利用L₀模度量稀疏性。提出了求解该模型的投影迭代硬阈值算法,并证明了算法的全局收敛性。该算法每一步都有闭式解,计算过程简洁高效。数值实验表明,方法在重建信号的视觉质量和量化指标方面均优于所对比的pFISTA方法。     

基于高阶滑模控制器的水下潜器运动控制

水下环境复杂多变,由于水流的不可预知性和多变性,潜器的水动力系数往往无法准确获取,使得依赖这种参数的潜器运动控制算法的应用受到了很大的局限。为了解决控制器对模型参数的依赖,设计了一种基于高阶滑模控制算法的模型无关控制器,并通过设置合理的过渡过程,解决了这种控制算法依赖初值的弊端。仿真结果表明,位置和姿态的控制能够快速的收敛,误差很小并且不依赖于初始条件,控制器需调节参数很少,并且算法简单,适用于工程的实际需要。     

深埋泥岩地层的最小主应力测定方法探讨

结合理论分析与实际资料统计,根据泥岩地层在水力压裂垂直裂缝剖面中所占的比例大小和瞬时停泵压力的变化趋势,获得了鄂尔多斯盆地西峰地区长8层(2 100 m)和靖边地区盒8地层(3 000 m)的泥岩最小主应力,这可为油气钻采、压裂设计以及测井、地震资料解释提供参考依据。间接证据表明,利用该方法求取的泥岩地层地应力更符合实际情况。     

基于三维层状介质模型地震层析成像正演网格研究

地震层析成像技术是研究地层速度结构的一个重要手段,该技术是利用震源到接收点的地震波走时或波形的观测数据,结合建立的数学物理模型来重建地层速度结构。而层析成像正演模拟的精度将直接影响反演出的速度与地层真实速度的拟合度,因此找出一种精度较高的层析成像正演模拟算法是非常重要的。本文针对三维层状介质,通过网格化模型,采用最小走时射线追踪方法展开层析成像正演研究,分别采用三角形网格、矩形网格和六边形网格进行模型参数化,并对比分析各自正演结果的精度,总结出针对三维层状介质模型的最佳网格划分方式。     

基于交替方向乘子算法的二维磁异常稀疏反演

磁异常反演是获取地下场源磁化率分布的重要手段之一,在地球勘探中扮演着重要角色.在磁异常反演中,对比光滑反演,稀疏反演的结果具有边界分明,物性参数分布集中的特点,更符合实际情况.针对稀疏反演,本文首先构建了具有代表性的基于L₁范数目标函数,利用交替方向乘子算法可分离凸函数的特点,将极小化L₁范数的优化问题分解为一系列的子问题,通过对子问题求解获得原问题的解;为了增强交替方向乘子算法的适应性,本文结合广义软阈值函数将交替方向乘子法推广于L_p（0 < p < 1）范数的反演中.为了验证本文提出的算法的有效性,采用了三种常规模型进行模拟实验.与基于L₂范数的反演算法进行实验对比,结果表明,本文算法得到了边界清晰,磁化率分布更集中的反演结果.最后,将基于交替方向乘子算法的L₁和L_p（0 < p < 1）范数的反演应用到青海省尕林格铁矿保护区获得的实际磁异常数据中,获得了较为符合实际地质情况的稀疏反演结果.

     

Numerical analysis of the symmetric methods

Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.