有色反演技术在少井区岩性体预测中的应用

阐述了有色反演的基本原理,并给出了反演流程。通过模型试算,明确有色反演关键参数的物理意义及取值范围,并在孤南洼陷实际工区应用,较好地刻画出沙三下早期低位扇体的横向展布范围。与常规地震属性相比,反演结果更加符合沉积规律,与实钻井岩性分布情况吻合程度高。表明有色反演在保留地震数据原始现象的基础上,分辨率明显提高,适合无井或少井区的地层或岩性解释。     

甘青地区马家窑文化遗址的地貌环境分析及其土地利用研究

马家窑文化是甘青地区新石器时期的主流文化,其时空演化和生业模式转变过程的研究取得了较大进展,但其不同发展阶段的地貌环境和水文条件的变化却仍不清晰。本文借助GIS方法和遗址域的概念,采用定量手段分析了马家窑文化不同时期遗址分布的地貌背景和汇水累积量的变化,并以此推论这些不同的地貌环境所反映的土地利用方式及其与生业模式之间的关联。分析结果表明,马家窑文化在不同时期的遗址域内有不同的平地等级和汇水累积量分布,这些指标所反映的土地开发、占用模式与当时的生业模式有很大关系。这项工作为理解不同气候背景下史前人类对生存空间和生业模式的选择机制提供了新的视角。     

半径-斜率域加权反假频地震数据重建

为减小地震数据缺失给地震后续处理工作带来的影响,需要对地震数据进行插值重建.针对反假频插值重建这个难点问题,进行了相关研究,并由此提出了一种改进的R-P（半径-斜率）域加权反假频地震数据插值重建方法.该方法将F-K（频率-波数）谱变换到R-P域,在R-P域设计一个权函数并将其作用于每次的迭代插值过程.通过模型数据和实际数据的测试,证明了该方法具有较好的反假频插值重建能力.     

特征高斯波包叠前深度偏移方法

高斯波包（Gaussian packet）传播算子可在局部时空域高效地计算局部波包的传播.高斯波包叠前深度偏移的基础是在Gabor变换域描述观测数据,再利用高斯波包传播算子计算炮点波场和检波点波场,两者相关即可得到偏移结果.利用炮道集的局部τ-p特征在Gabor变换域表达观测数据,可以仅关注部分高斯波包框架函数上的数据投影,这样既实现了波场的压缩存储,同时可利用高斯波包传播算子反传框架函数以实现整个炮道集的快速反传.这些综合了观测数据局部τ-p特征的高斯波包函数称为特征高斯波包（characteristic Gaussian packet,CGP）,相应的波场反传称为特征高斯波包反传.理论及数值分析证明了上述特征高斯波包反传方法是有效且快速的.炮点正传波场也利用高斯波包传播算子模拟.利用互相关成像条件可实现特征高斯波包叠前深度偏移（characteristic Gaussian packet pre-stack depth migration,CGPM）.由于高斯波包传播算子描述了局部方向及局部空间的波场,所以CGPM可以自然地提取角度域成像道集（ADCIG）,并易于实现面向目标叠前深度偏移,从而作为偏移引擎为偏移速度分析（MVA）服务.数值实验证明了CGPM和面向目标CGPM的有效性和实用性.     

用扩散电场法估算CHTEM-I系统的探测深度

时间域航空电磁系统的探测深度是一项关键技术参数,在系统研制设计和资料解释中都有重要作用.当目标层与围岩电性差异不大或目标层引起的异常响应不明显时,常规的估算探测深度方法不再适用,为克服这种缺陷,本文提出了一种系统探测深度的估算方法,该方法通过模拟电场在均匀半空间模型和层状模型中的扩散过程,记录电场最大幅值在地下的瞬时位置,将系统的探测深度定义为观测到的响应值等于给定的噪声水平时,对应时刻的地下介质中感应电场幅度的最大值对应的深度.论文以我国自主研制的直升机时间域航空电磁系统CHTEM-I为例,根据对不同条件下的计算结果的分析,给出了探测深度与模型电导率、飞行高度和噪声水平之间的关系曲线,并总结出一些能够提高系统探测深度的建议.结果表明半空间模型电导率在0.000295~0.0422 S/m范围内,CHTEM-I系统可以在噪声水平与t－0.5成正比条件下达到300 m的探测深度.本文方法的估算结果不受模型层厚薄,或目标层与围岩电性差异小等因素影响,因而具有较高的实用性.文中的结论对时间域航空电磁系统设计有理论指导意义,也可用于其它时间域电磁勘探系统探测深度的估算.     

基于核主成分分析的时间域航空电磁去噪方法

时间域航空电磁数据往往在测量过程中受到天然和人文噪声的干扰.如果不能很好滤除这些电磁噪声,那么将会降低资料质量、影响反演的精度,甚至获得错误的解释结果.本文提出了一种基于核主成分分析的去噪方法,通过核主成分分析提取叠加后衰减曲线的主成分,然后使用能量占比方法分离反映地下介质的有效信号和噪声,最后使用反映地下介质的特定成分进行重构.本文所推荐的去噪方法不仅能剔除天然噪声,例如天电产生的尖脉冲或者振荡,而且能有效地抑制人文噪声.分别使用基于核主成分分析的去噪方法,以及AeroTEM软件的处理方法对同样的吊舱式时间域直升机航空电磁勘查系统实测数据进行处理,并比较其结果.处理结果表明：所推荐的去噪方法要优于AeroTEM软件.     

OBC与拖缆数据的子波域匹配滤波方法研究

海底电缆（OBC）采集方式的出现,导致OBC采集与拖缆采集数据进行拼接成为三维连片地震资料处理中的一个重要环节.由于采集方式不同,数据在振幅、频率、相位等方面均存在较大差异.为了消除和减小这些差异,为后续的地质解释和油气预测提供便利,本文提出了一种在子波域的1范数匹配滤波方法来计算匹配滤波因子,以用来提高OBC数据和拖缆数据的匹配精度,改善连片处理的数据品质,为后续的地质解释和油气预测提供便利.通过实际资料检验,子波域匹配滤波方法在OBC与拖缆数据拼接应用中取得了较好的效果.     

GPS -InSAR合成方法用于地面沉降监测的可行性展望

GPS技术在地面沉降领域的应用已日臻成熟并完善;InSAR技术在地面沉降领域的应用已经开始但还未成熟,目前是研究热点。GPS—InSAR合成方法是将GPS高时间分辨率和InSAR高空间分辨率进行有效的统一,使两种技术达到互补,从而发挥各自的优势。本文指出了GPS—InSAR合成方法用于地面沉降监测的可行性,并对两种实施方案进行了探讨。最后对西安地区开展GPS—InSAR研究进行了展望。     

空间-波数域三维大地电磁场积分方程法数值模拟

地球物理勘探中, 地下复杂地质体电磁场三维数值模拟计算量大、存储要求高.针对这一问题, 本文提出一种空间-波数域三维电磁场数值模拟方法, 该方法利用电磁场积分方程为卷积的特点, 将积分沿水平方向进行二维傅里叶变换, 将电磁场三维空间域卷积问题转换为不同波数相对独立的多个一维积分问题, 由此计算量和存储需求大大减少, 易于实现并行计算.采用有限单元法中的形函数进行一维数值积分, 一维积分离散为多个单元积分之和, 每个单元采用二次形函数表征电流变化, 可得出单元积分的解析表达式; 保留垂向为空间域, 优势之一在于可根据实际情况合理调整单元疏密程度, 准确模拟任意复杂地形和导电率异常, 兼顾计算精度与计算效率, 优势之二是用形函数拟合求得积分的解析解, 计算精度和效率高; 最后引入压缩算子, 采用迭代求得电磁场的数值解.本文方法充分利用不同波数之间一维积分高度并行性、一维形函数积分的高精度及快速傅里叶变换的高效性, 实现电磁场三维高效高精度数值模拟.设计模型将本文方法数值模拟结果和软件INTEM3D的数值模拟结果对比, 验证了方法的正确性; 设计高阻和低阻异常体, 研究了异常电导率与背景电导率差异对迭代收敛速度的影响; 改变计算规模, 随着计算网格增多, 算法耗时与存储呈近似线性增长; 设计复杂模型与目前主流数值模拟方法对比, 本文算法速度快1个数量级以上, 且计算规模越大, 算法效率的优势越明显.研究结果表明, 本文提出的空间-波数域三维电磁场数值模拟方法理论和方法正确, 计算效率高, 对计算机存储要求低, 算法高度并行, 适合任意复杂条件大规模三维电磁场高效、高精度数值模拟.

     

各向异性介质空间-波数混合域直流电阻率法三维数值模拟研究

将各向异性介质分成各向同性背景介质和各向异性异常介质,并提出一种空间-波数混合域方法实现了各向异性介质下直流电阻率法的三维数值模拟.不同于传统直流电阻率数值模拟方法,本文算法直接对空间域异常电位满足的偏微分方程沿水平方向进行二维傅里叶变换,使水平方向转换成波数域,保留垂直方向为空间域,可根据地下介质电流密度变化的快慢灵活剖分.这样可把空间域异常电位满足的三维偏微分方程转化成不同波数满足的一维常微分方程,把一个大规模三维数值模拟问题分解为多个一维数值模拟问题,利用一维有限单元法求解方程组,并通过采用压缩算子迭代计算,最终获得较为精确的数值解.与自适应有限单元法对比验证了本文算法的正确性;测试了算法的收敛性,结果表明在满足精度要求的情况下,算法的收敛性只与异常体和围岩之间的电导率差异相关,而与异常体大小和埋深无关;分析了算法计算效率,结果表明算法的计算效率与剖分网格节点成线性关系,算法可在微型计算机中较快计算出剖分节点总数超过千万的各向异性模型的结果;设计简单观测系统并验证其具有反映地下各向异性结构特性的能力;最后模拟异常体沿着不同方向旋转不同角度时的响应特征,对比分析可知异常体为各向异性情况下,沿不同方向旋转相同角度或沿同一个方向旋转不同角度,在同一个位置上测量得到的视电阻率有差异.