兰州地磁台GPS差分方位角的测量

为了适应高采样率、高分辨率、高精度的新型地磁仪器对磁偏角和磁倾角观测的需要,应用GPS测定了改造后的兰州地磁台的7个观测墩（一个临时观测墩）至2个方位标志的方位角,共取得了13个方位角参数,使兰州地磁台在地磁偏角测量或仪器验收比测中可供6套仪器同时观测,并可做到全天候观测。     

倾斜各向异性地层中多分量电磁波测井响应三维时域有限差分（FDTD）算法

利用三维时域有限差分技术（FDTD）研究建立倾斜各向异性地层中多分量电磁测井响应的数值模拟算法.首先,利用Yee氏交错差分格式和线性内插公式对各向异性介质中时域Maxwell方程进行离散化,得到各向异性介质中各个交错节点上电磁场的时域递推公式,同时将面积加权平均与旋转矩阵技术相结合确定各个网格节点的等效电导率张量,此外,单轴各向异性吸收边界条件（UPML）被用于降低求解区域外边界反射.在此基础上,利用单频正弦磁偶源计算时域电磁场空间分布以及特殊的双方程算法（2E2U）从数值结果中快速提取电磁场幅度和相位,确定多分量电磁测井仪器的三维响应.最后,通过数值结果证明该算法的有效性,并系统考察不同层厚、地层倾角以及侵入带半径等变化对多分量电磁测井响应的影响.     

交错网格与旋转交错网格对VTI介质波场分离的影响分析

利用交错网格有限差分和旋转交错网格有限差分进行各向异性介质弹性波场数值模拟时, 质点振动速度分量与应力张量的网格节点定义方式均不相同, 从而对各向异性波场分离效果产生不同的影响. 针对这一问题, 本文以具有垂直对称轴的横向各向同性(VTI)介质的波场分离为例, 首先分析了两种网格的参数定义方式以及VTI介质波场的分离过程; 其次, 详细研究和分析了这两种网格的参数定义方式对各向异性介质波场分离的影响, 并依据波前面连续性以及波场分离效果等方面, 通过数值模拟实验对该影响进行分析验证. 结果表明, 旋转交错网格的参数定义方式更有利于进行各向异性介质波场数值模拟和波场分离.      

三维任意各向异性介质中海洋可控源电磁法正演研究

由于海底介质受沉积环境的影响,层理发育呈现明显各向异性特征.对于海洋可控源电磁法各向异性的研究以往主要局限于一维和二维模型,为更深入了解复杂情况下海底各向异性对海洋可控源电磁响应的影响规律,本文开展三维任意各向异性介质中海洋可控源电磁法正演研究.采用交错网格有限差分技术,通过对任意各向异性介质电导率张量实行体积和空间电流密度平均,完成海洋可控源电磁二次散射电场的离散化,成功实现任意各向异性介质中海洋可控源电磁正演模拟.通过对几种典型各向异性电性模型条件下海洋电磁电场多分量响应及分布特征和各向同性情况的对比分析,总结电各向异性对海洋电磁响应的影响规律和识别方法.本文算法研究及算例可为海洋可控源电磁数据精细化处理解释提供技术支撑.     

Lebedev网格改进差分系数TTI介质正演模拟方法研究

本文采用一种新的交错网格-Lebedev网格(LG)进行TTI介质的正演模拟研究,避免了Virieux标准交错网格(SSG)算法在处理TTI、单斜等各向异性介质时波场插值引入的数值误差,提高了模拟精度.在方法实现过程中,本文针对有限差分正演模拟面临的网格频散与边界反射两个关键性问题分别做了优化,并通过模型试算验证了它们的有效性与可行性:(1)结合最小二乘思想推导出新的频散改进差分系数(DIC),该系数比Taylor系数更能有效地压制粗网格引起的数值频散,可以节约内存,提高计算效率;(2)将分裂的多轴完全匹配层(M-PML)吸收边界条件引入到LG算法中,解决了传统PML边界条件在某些各向异性介质中的不稳定现象并且具有较好的边界吸收效果.     

局域双曲线Radon变换压制散射波研究与应用

一般而言,由于地下非均质体的存在所产生的二次波源,由它再生成新的波场,叫散射波场。目前,将散射波作为有效波来成像,已开始在溶洞和裂缝等特殊地质体的识别中得到应用。但对于野外资料采集来说,地表复杂地区,如戈壁、砾石区和山前带,大量存在的散射波却是干扰波,它们的存在会严重影响资料的品质,而其研究与实际应用国内外还很少。因此,通过正演模拟,分析散射波的基本特征,在此基础上研究散射波的去噪方法显得十分必要。本文从地震波运动学时距关系出发,研究了反射波和散射波的几何特征;然后用有限差分正演,模拟了散射波场,用理论模型研究并测试了局域双曲线Radon变换散射波去噪新方法。对于实际炮集资料,分析了F-K滤波方法压制散射噪声的局限,采用局域双曲Radon变换有效地去除了炮集中存在的散射噪声,取得了较好的应用效果。     

海底天然气水合物OBS多分量地震正演模拟（英文）

根据Ecker的水合物沉积物的三种微观模式,计算含水合物沉积层和含游离气沉积物的弹性模量,分析对比了水合物的不同微观模式、不同水合物饱和度以及不同游离气饱和度对沉积物弹性模量的影响;从纵横波分离的弹性波动方程出发,采用交错网格空间有限差分方法模拟地震波在海底天然气水合物沉积地层的传播,得到纵、横波的海底地震（OBS）共接收点道集。数值算例表明,当水合物作为流体的一部分或胶结颗粒骨架时,仅纵波记录上存在BSR;当水合物胶结颗粒接触,纵、横波记录上均存在BSR。并且,OBS会接收到上行纵波和上行横波在海底界面形成的转换波,干扰横波记录上BSR的识别。     

Atmospheric correction of HJ-1 A/B CCD over land: Land surface reflectance calculation for geographical information product

This paper proposed a method to retrieve the land surface reflectance from the HJ-1A/B CCD data. The aerosol optical depth（AOD）, the most important factor affecting the atmospheric correction of CCD images at all bands, is proposed to retrieve from the CCD imagery by the approach of dense dark vegetation（DDV） method. A look-up table in terms of the transmittances, the path radiances and the atmospheric spherical albedo as functions of the AOD was established for a variety of sun-sensor geometry and aerosol loadings. The atmospheric correction is then achieved with the look-up table and the MODIS surface reflectance output（MOD09） as the priori datasets. Based on the retrieved AOD and the look-up table of atmospheric correction coefficients, the land surface reflectance was retrieved for the HJ-1A/B data according to the atmospheric radiative transfer equation. Some in-situ measurement Data for Yanzhou of Shandong province in East China and MODIS land surface reflectance products MOD09 are used to preliminarily validate the proposed method. The results show that the proposed method can remove effectively the atmospheric contributions, and the overall accuracy of the retrieval land surface reflectance can be improved substantially.     

DGPS水下立体定位系统的解析求解算法

差分GPS水下立体定位系统通过测量目标到两个浮标的距离差来进行空间定位,由于这种双曲线定位模型在进行迭代计算时容易发散,所以本文根据空间浮标和目标的空间几何关系,提出了一种空间双曲线交会的直接算法,该方法能够直接求解目标点的空间坐标,且不存在定位结果的发散问题。     

双线偏振雷达差分传播相移的质量控制

利用中国气象科学研究院灾害天气国家重点实验室的车载C波段双线偏振多普勒雷达 (C-band Polarimetric Doppler Radar on Wheel,CPDRW) 的外场试验,在统计分析降水、地物回波差分传播相移Φ_DP数据的差别与信噪比关系等基础上,提出了一套数据分析和处理的方法。该方法通过Φ_DP的异常波动并结合回波的强度Z_H和速度V_r信息将地物回波信号分离出来,在降水估测或衰减订正等定量应用时将其剔除。对于气象回波则根据信噪比及零滞后互相关系数ρ_HV(0) 将Φ_DP资料分为较好、较差和差3类。对于较好数据直接进行后续的预处理,对于较差数据先订正后处理,而对于差数据将其剔除以保证Φ_DP资料的整体质量。经过大量资料的验证,该方法在最大程度上保留气象信息的同时也保证了Φ_DP资料的质量,并能估算出质量较高的差分传播相移率K_DP资料。