基于基本多项式的GPS精密星历插值方法研究

采用基本多项式插值方法对GPS精密星历进行了插值,该方法利用Householder变换来避免对病态的法方程系数矩阵直接求逆而导致的插值计算错误,并对高阶基本多项式插值过程中出现的龙格现象进行了讨论,给出了多项式拟合区间长度、有效区间长度和基本插值多项式阶数的选取准则,最后对插值的结果进行比较分析。结果表明,对于15分钟采样率的精密星历,基本多项式插值精度达到毫米级,能够满足精密定位的需要。     

资源三号卫星不同空间分辨率图像斑状植被检测比较研究

黄河三角洲滨海区域大量分布着由柽柳、芦苇、翅碱蓬、二色补血草、白茅等等组成的斑状植被,多呈棋盘状规则分布,胜利油田地震勘探被认为是产生这一现象的原因之一。深入研究其分布格局并探讨其形成机理,对于区域湿地生态系统植被恢复至关重要。为此,首先就需要对斑状植被进行有效地检测和制图。本文评价了我国首颗民用高分辨率光学传输型立体测图卫星（资源三号卫星）不同空间分辨率图像检测黄河三角洲斑状植被的能力,结果表明：在利用霍夫曼圆检测方法下,前视、后视、正视全色图像和正视多光谱图像都能很好地检测斑状植被,空间分辨率高的图像有更高的检测精度。     

北斗非GEO卫星信号码并行捕获

2012年12月,北斗二代的接口控制文件正式公布。北斗B1I频点的非GEO卫星上调制20bit的NH码,因此每1ms信号可能存在比特跳变。借鉴传统GPS并行码相位捕获算法,本文研究一种基于补零FFT的并行码相位算法。该算法能够实现1-20ms相干积分,利用非相干累加提高信噪比,并采用一种对信号强度和采样率不敏感的判决变量。仿真表明,通过相干积分和非相干累加,该算法使用400ms数据可以捕获24dBHz的弱信号。     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

多次波压制技术在南海北部陆缘深水区的应用

南海北部洋陆过渡带水深变化剧烈,海底地形崎岖,构造复杂多变,地震勘探中多次波十分发育,常规处理难以取得很好的成像效果.由于方法的局限性,任何一种多次波压制方法都很难彻底地压制全部多次波.本文通过对国家基金委共享航次在南海北部陆缘采集的二维反射地震资料的分析,总结了不同多次波压制技术的优缺点和适用范围,结合滤波类方法和波动理论类方法,综合应用了表面相关多次波衰减法(SRME,Surface Related Multiple Elimination)、τ-p域预测反褶积、高分辨率抛物线Radon变换以及叠后压制残余多次波等方法去除多次波,有效突出了一次反射波,证明了方法的有效性和实用性,同时也形成了针对南海深水资料多次波压制的一种有效的处理流程和方法组合.     

用折线逼近法计算余弦变换的频率范围讨论

要想提高余弦变换的精度,可以减小折线区间,增加折线的数量,最大限度地达到折线代替函数曲线的目的,然而在实际应用过程中我们发现,合理正确的选择变换的频率范围对余弦变换的计算精度会产生决定性的影响,本文目的就是为了减小因频率区间选取的不同而造成余弦变换结果的误差,使瞬变响应的计算精度达到最高.本文以大回线和磁偶极子为例,在电阻率、线圈半径、收发距不同因素的影响下,求其频率域响应,利用余弦变换计算其瞬变响应,通过与解析解的对比,来讨论使余弦变换结果达到最高精度所需要的最佳fH、fL.结果发现,随着电阻率的增加,余弦变换精度对频带的依赖性增强;线圈半径、收发距增加,余弦变换精度对频带依赖性减弱.无论哪种情况,最佳的频率范围基本上都是在10-6~107 Hz,而且当频宽不变时(logfH-logfL≥12),固定频带fH,大回线源的计算精度可以得到保证,固定fL,磁偶极子的计算精度可以得到保证.从而我们总结出一个在不同因素影响下的余弦变换的频率选取规律,大大降低了余弦变换结果因频段的选取而造成的误差.     

λ-f域加权抛物Radon变换地震数据重建方法研究

大多数地震处理技术都要求地震数据具备完整性和规则性,然而,由于诸多因素的影响,地震勘探所采集到的资料普遍存在数据缺失,需要对地震数据进行重建.本文在传统Radon变换重建的基础上提出一种λ-f域加权抛物Radon变换的地震数据重建方法.通过引入新变量λ,消除了Radon变换算子对频率的依赖,使得Radon变换算子及算子的逆只需计算一次,显著提高了计算效率.同时,在λ-f域抛物Radon变换迭代计算过程中引入变化的权系数,更好地实现了λ-f域的能量聚焦.理论模型及实际数据试算表明,文中方法对地震数据重建精度较高,单道对比吻合较好.     

运用分形方法模拟地震勘探中的三维随机介质

实际的地震勘探遇到的是复杂随机介质,常规层状介质地震波场理论已经不能适应于高精度地震工作,需要进一步发展非均匀介质地震波场理论.本文首次给出了运用分形方法建立三维随机介质的具体算法;其次引入了三维锥形函数来压制建模时产生的离散误差;最后,通过建立几种三维随机介质模型来讨论尺度系数和豪斯特指数对模拟结果的影响,并且简单构建了区域多尺度随机介质模型.模拟的结果表明,分形方法可以更方便、灵活地模拟三维随机介质.针对南海某地区生物礁储层非均质性强的特点,在建立生物礁速度模型时,对目标区充填了多尺度的分形随机介质,波动方程正演模拟结果表明,分形随机介质可以用来描述生物礁等非均质性储层.     

基于Curvelet变换的探地雷达资料噪声衰减方法(英文)

噪声衰减是探地雷达信号处理中的关键问题之一。当探测目标埋藏深度比较浅时,其反射信号与直耦信号和地面回波信号相互重叠,直接影响目标反射波到达时刻的检测及目标的正确定位。针对这个问题,本文提出了一种基于Curvelet变换的噪声衰减方法。通过对理论数值模拟数据和实测数据的处理,以及与平均消去法和二维连续小波该方法处理结果的对比,验证了该方法的可行性和有效性。处理结果显示,该方法不仅可以去除背景噪声、同时可以衰减倾斜相关的相干干扰和数据中的随机噪声。与二维连续小波变换方法相比有更高的计算效率。     

“摩尔空间”及其在先存构造活动性预测中的应用

应用"先存构造活动性准则",结合数值计算分析,建立了"摩尔空间"中先存构造面的极点(σn,τn)与其方位(θ′和φ′)的一一对应关系,从而可以十分直观和便捷地确定三维空间上任意先存构造面上的正应力(σn)和剪应力(τn).在此基础上,阐述了先存构造面的活动性在"摩尔空间"中的表达方式,以及预测分析不同先存构造面相对活动性的技术方法和步骤.最后,应用"摩尔空间"预测了砂箱实验模型中先存构造面的相对活动性,以及活动断层的相对活动性.结果表明,利用摩尔空间可以十分直观和便捷地表达三轴应力状态下任意面上的正应力和剪应力、确定先存构造的活动性系数、预测不同先存构造的相对活动性;应用"摩尔空间"的预测结果得到了实验结果和实际测量资料的证实."摩尔空间"为任意三轴应力状态(主应力轴倾斜)的应力分析、先存构造活动性分析提供了便捷的途径,由于它的直观性和方便性,预期具有广阔的应用前景.