基于反射地震记录变子波模型提高地震记录分辨率

本文给出了地震记录变子波模型的一种近似数学表达式．基于该表达式研究了反射系数序列不满足白噪假设和子波在地下传播时发生变化这两种情况下地震道谱的组成及结构,讨论了谱白化及反褶积方法在这两种情况下效果不佳的原因．然后基于变子波模型,提出了一种新的提高地震记录分辨率的方法：第一步,用自适应于地震记录的Gabor分子窗把地震记录恰当地划分成若干片断,每段内信号近似平稳,然后将地震记录变换到时间-频率域;第二步,在变换域对每个分子窗内信号的振幅谱进行处理以拓宽频带;最后把处理后的时间-频率域函数反变换回时间域得到提高分辨率后的结果．本文提出的方法具有能较好地适用于反射系数不满足白噪假设的情况及提高分辨率后的地震记录能较好地保持原地震记录的相对能量关系等优点,模型和实际资料算例结果均表明,本文方法在拓宽地震资料频带及保持地震记录局部能量相对关系方面均明显优于谱白化方法．     

自适应噪声消除器在地电场观测数据处理中的应用

在地电场和电阻率同场地观测的台站, 地电场观测必然会受到地电阻率观测人工供电的干扰, 尤其是当地电场观测的采样率提高到1次/s以上时, 这种干扰更是无可避免, 而且不易消除。 该文阐述了应用自适应噪声消除器对受种干扰的地电场观测数据进行处理, 消除观测数据中的干扰噪声, 从而获得地电场真实信号波形。 对实际观测数据的处理结果表明, 采用此方法可以很好的消除秒采样地电场观测数据中地电阻率供电产生的干扰。 这种方法还可用于消除其他观测方法中存在的有规律的干扰。     

可控源电磁场三维自适应矢量有限元正演模拟

本文实现了可控源电磁(CSEM)场三维自适应矢量有限元正演算法,该算法采用非结构四面体单元进行三维网格剖分,能够真实模拟地形起伏和复杂电性异常体.采用一次场和二次场分离的方式计算电磁场响应,能够有效解决有限元模拟中的源点奇异性,提高场源附近电磁场数值精度,其中一次场利用CSEM一维正演算法解析求得,二次场采用矢量有限元方法求得.并利用基于后验误差估计的自适应网格细化算法指导网格细化,以减少人为设计网格导致的误差.通过一维和三维模型的数值模拟,验证了本文算法的有效性:一维模型有限元解与解析解吻合得很好,电磁场振幅相对误差在1%左右,相位差整体小于1°;三维模型有限元解与有限体积解吻合得也很好.模拟了一个含三维倾斜板状异常体的可控源电磁响应,表明了本文算法模拟复杂地电结构电磁场的能力和有效性.     

基于多次波和自适应预测误差滤波器的地震近炮检距数据重建方法

地震数据采集过程中,受野外施工条件的制约,往往很难获得完整的地震波场,尤其近炮检距数据的缺失尤为严重.当前,很多地震数据处理方法的应用都依赖于近炮检距数据,如何对近炮检距的缺失数据进行重建,是一个重要的研究课题.本文通过多次波和一次波的互相关构建虚拟一次波,利用数据本身的波场信息,对缺失的近炮检距数据进行插值重建.由于通过多次波构建的一次波与真实的一次波存在振幅和相位方面的差异,提出通过基于非平稳自回归过程的自适应预测误差滤波器来表征虚拟一次波的能量谱,利用最小二乘反演方法重建近炮检距缺失数据,自适应预测误差滤波器通过求解正则化约束下的数学欠定问题来实现局部自适应特征.通过对Sigsbee2B模型和实际数据的测试结果表明新方法可以合理地重建复杂的近炮检距缺失数据.     

基于自适应有限元的复电阻率法2.5维复杂构造电磁场模拟

复电阻率法在矿产、油气勘探调查中发挥着重要作用.为了认识复杂构造的复电阻率法电磁场的变化规律,本文基于自适应有限元方法,采用非结构化网格,引入Cole-Cole模型,实现了电偶源2.5D复电阻率法电磁场正演,可以模拟复杂地形和地电结构,正演结果更符合野外实际地质情况.通过将本文的计算结果与半空间模型解析解、层状介质和起伏模型结果进行对比,验证了本文算法的正确性.最后,基于复杂地电模型,通过正演模拟,系统分析了地形、激电参数、复杂构造对复电阻率法电磁场的影响特征.     

垂直井眼中三维阵列感应测井资料自适应井眼校正

三维感应测井受井眼环境影响很大,如何消除这些不利影响是有效获取各向异性地层电导率真值并提高资料处理和解释质量的一项重要工作.本文以一维柱状各向异性地层模型为基础,通过最小平方反演技术从实际测量资料中快速反演井径、偏心距、泥浆电导率、地层水平电导率和各向异性系数等模型参数,设法减小井眼环境对测量结果的影响.首先,利用仪器偏心情况下不同旋转角电导率张量的关系方程,建立旋转角提取与校正方法,获取无旋转角三维感应测井数据.在此基础上,根据事先计算好的井眼校正库并结合多维非线性有限元逼近技术,建立快速计算三维感应测井响应以及Jacobi导数矩阵的插值公式,并利用自适应正则化迭代反演技术不断修改模型参数,实现理论合成资料与输入资料的最佳拟合.最后,通过反演出的模型参数计算三维感应测井资料所有分量的校正量,实现三维感应资料的井眼校正.理论模型和实际资料的处理结果均证实,在低阻井眼泥浆情况下,一维柱状模型中的反演方法仍然可以提取出较可靠的地层水平和垂直电导率,且井眼环境校正对于受井眼环境影响相对较小的测井曲线具有非常好的校正效果.     

基于电流密度连续性条件的直流电阻率各向异性问题自适应有限元模拟

直流电阻率法被广泛应用在工程和环境及水文地球物理、野外采矿、地热探测等领域.地下岩石常具有层理面和裂缝等具有方向依赖性的结构,岩石电导率常常呈现各向异性特征,因此研究复杂直流电阻率各向异性问题的高精度正演算法具有迫切的理论和学术需求.本文利用面向目标的自适应有限元算法和非结构化网格相结合的方式,解决了带地形任意复杂直流电阻率各向异性问题的高精度正演这一难题.有别于前人的研究成果,本文提出了一种特别的二次虚拟场算法实现带源的任意起伏地形问题模拟;另外,本文第一次基于电流密度连续性条件构建适合直流电阻率各向异性问题的后验误差估计算法,有效地驱动面向目标有限元网格自适应加密过程.最后,通过三组电阻率各向异性模型验证本文提出算法的正确性和适应性,测试结果表明：对于任意复杂直流电阻率各向异性问题,本文提出的算法具有精度高、适应性强等特点;另外,我们还发现电流密度连续性条件可用于设计直流电阻率问题的有效后验误差估计算法.     

中国西南地区GPS大气水汽转换系数模型精化研究

利用中国西南地区19个探空站2011～2014年数据,通过积分法计算大气水汽转换系数K。采用2011～2013年的K值对Emardson模型进行精化,分别构建西南地区不顾及和顾及高程因子的K值模型Emardson-I和Emardson-H。利用2014年积分计算的K值检验这两种模型的预报精度,结果表明:1)相对于Emardson-I模型,Emardson-H模型表现出更高的K值预报精度和更好的适应性;2)在高海拔地区,Emardson-H模型预报精度明显优于Emardson-I模型,表明高程因子是影响高海拔地区K值计算的重要因素。将两种新模型用于拉萨站GPS大气水汽反演,Emardson-H模型表现出更优的反演精度,两种模型的反演精度均优于2 mm。     

Structure and Mössbauer spectroscopy of barbosalite Fe2+Fe3+ 2 (PO4)2(OH)2 between 80 K and 300 K

Natural barbosalite Fe²⁺Fe³⁺ ₂ (PO₄)₂(OH)₂ from Bull Moose Mine, South Dakota, U.S.A., having ideal composition, was investigated with single crystal X-ray diffraction techniques, Mössbauer spectroscopy and SQUID magnetometry to redetermine crystal structure, valence state of iron and evolution of ⁵⁷Fe Mössbauer parameter and to propose the magnetic structure at low temperatures. At 298?K the title compound is monoclinic, space group P2₁/n, a _o ?= 7.3294(16)?Å, b _o ?=?7.4921(17)?Å, c _o ?=?7.4148 (18)?Å, β?=?118.43(3)°, Z?=?2. No crystallographic phase transition was observed between 298?K and 110?K. Slight discontinuities in the temperature dependence of lattice parameters and bond angles in the range between 150?K and 180?K are ascribed to the magnetic phase transition of the title compound. At 298?K the Mössbauer spectrum of the barbosalite shows two paramagnetic components, typical for Fe²⁺ and Fe³⁺ in octahedral coordination; the area ratio Fe³⁺/Fe²⁺ is exactly two, corresponding to the ideal value. Both the Fe²⁺ and the Fe³⁺ sublattice order magnetically below 173?K and exhibit a fully developed magnetic pattern at 160?K. The electric field gradient at the Fe²⁺ site is distorted from axial symmetry with the direction of the magnetic field nearly perpendicular to V_zz, the main component of the electric field gradient. The temperature dependent magnetic susceptibility exhibits strong antiferromagnetic ordering within the corner-sharing Fe³⁺-chains parallel to [101], whereas ferromagnetic coupling is assumed within the face-sharing [1?1?0] and [?1?1?0] Fe³⁺-Fe²⁺-Fe³⁺ trimer, connecting the Fe³⁺-chains to each other.     

改进的Sage—Husa自适应滤波在GPS动态导航定位中的应用

本文利用Kalman滤波方法对动态测量进行数据处理,由于高动态的GPS测量,不易确定系统动态噪声和观测噪声.同时标准的Kalman滤波在应用过程中由于状态模型确定的误差存在,滤波效果不佳.因此本文结合动态导航的实时性和高动态性,建立了动态导航系统中滤波状态方程和观测方程,采用改进的Sage-Husa自适应滤波对来进行实时定位数据处理,利用已有测量数据进行了实例分析.改进的Sage-Husa自适应滤波在计算过程中计算量小,结果稳定,有较强的自适应性.