利用共轭梯度算法的电阻率三维有限元正演

引入对称超松弛预条件共轭梯度(SSOR PCG)迭代算法求解电阻率三维有限元计算形成的大型线性方程组，并有机结合系数矩阵的稀疏存储模式，使得三维有限元正演计算的速度大大提高而内存需求则大大减少. 该算法可方便地应用于求三维异常电位，在保持快速计算的基础上，正演计算的精度显著提高.     

多重网格在二维泊松方程有限元分析中的应用

本文简要介绍多重网格(MG)算法的基本原理及基本步骤,然后将多重网格算法引入有限单元中,对二维泊松方程进行求解.单元数尺度从8×8逐次增加至1024×1024,并与单重网格中高斯-赛德尔迭代法(GS)、共轭梯度法(CG) 在程序运行时间以及迭代次数方面进行比较.结果表明MG在计算速度和迭代次数都明显优于GS、CG方法.在1024×1024网格中,MG不仅比GS快500多倍,比CG快60多倍,而且与理论解的误差更小.     

反假频非均匀地震数据重建方法研究

研究基于Fourier变换的数据重建方法,既能进行非均匀采样数据重建,又可以去除空间假频. 将不规则采样数据重建问题归结为信息重建的地球物理反演问题,采用最小二乘方法从观测的稀疏或不规则数据反演模型空间完全信息. 在求解信息重建反演问题时,引入DFT 加权范数规则化策略,采用预条件共轭梯度法（PCG）求解,保证解的稳定性和收敛速度. 处理线性同相轴假频问题时,根据采样定理,引入线性预测方法,采用Yule Walker方程由带限信号的无假频低频功率谱预测高频功率谱,达到反假频目的. 本文研究了均匀采样数据内插,非均匀采样数据重建,非均匀分布高频信息重建等方面问题,数值试验取得较好效果.     

重力梯度全张量数据三维共轭梯度聚焦反演

随着全张量重力梯度测量技术的日趋成熟和应用领域的不断扩大,重力梯度全张量数据的三维反演技术越来越受重视.本文利用剖分单元之间几何架构等效性,实现了重力梯度全张量场三维正演快速计算和导数矩阵优化存储.并将积分灵敏度、粗糙度和最小支撑泛函约束以及参考模型和模型参数界限约束引入到目标函数中,采用共轭梯度法进行反演迭代,实现了重力梯度全张量数据三维快速正反演计算.多种模型的反演试算表明:本文提出的反演算法的可靠性和稳定性较好,并且算法速度快、占用内存低且易于并行化.     

穿越南沙礼乐滩的海底地震仪广角地震试验

本文对穿越礼乐滩东北部向西北方向延伸进入中央海盆长369 km的广角地震剖面OBS973-2进行了反演研究,以期了解南海南部陆缘的地壳结构,同时探讨南、北陆缘的共轭问题.结果表明OBS973-2剖面的速度模型中三个沉积层的速度分别为1.8~2.0 km/s、2.0~2.7 km/s和3.5~4.0 km/s;沿剖面沉积层总体上较薄或缺失.礼乐滩上地壳厚约9~10 km , 速度为5.5~6.4 km/s,顶部存在小型火山;下地壳厚约11 km,速度为6.6~7.1 km/s. 过渡区和海盆的上地壳厚约4~5 km,速度为5.9~6.1 km/s;下地壳厚约2~4 km,速度为6.6~6.9 km/s.从总体上看,海盆和过渡区的地壳厚度偏小,显示了拉伸减薄作用,速度分层显示陆壳比较典型,而洋壳和过渡壳的上地壳速度比标准洋壳偏高.莫霍面总体从陆坡向海盆方向较快速地抬升,在礼乐滩埋深约23 km, 在海盆中的埋深8~12 km,海盆中的莫霍面顶部速度8.0 km/s,要明显小于礼乐滩下方的8.2 km/s.将OBS973-2剖面与北部陆缘的OBS2006-1剖面进行了对比,结果表明两者速度结构极为类似,讨论认为礼乐滩与中沙块体互为共轭.进而根据西北次海盆的宽度、礼乐滩与中沙块体的距离及前人提出的南海海盆的扩张时代,计算了扩张速率.     

山东姜家窑金矿岩体内成矿构造波动原理分析

姜家窑金矿位于东西向胶东复背斜西端,北东向招平衡裂带中段、矿区地 以断裂为界,下盘为玲珑化岗岩,上盘为太古宇胶东群变质岩。人的招平断裂带下盘入宇型次级断裂较发育,且延伸到花岗岩岩体内部,在二者所构成的锐角扇形区的岩体中,形成一系列具有等间距分布的岩石破碎块段,靠民岩石破碎的直接原因是岩体现人的共轭节理。同时根据地质力学原理,其锐角扇形区为应力拉张区,有利于成矿,所以姜家窑金矿的矿体除了沿招平断裂带     

塔河油田奥陶系走滑断裂体系特征及成因机制

塔河油田奥陶系“断溶体”油藏是塔里木盆地油气勘探开发的重要对象。本文对塔河油田10区西、于奇西、于奇8井三块三维地震资料开展了精细的走滑断裂构造解析,并在此基础上对塔河油田走滑断裂体系的形成演化过程与成因机制有新认识。研究结果表明：(1)除了在塔河油田南部10区西、托普台等区块发育典型的“X”型共轭走滑断裂体系外,首次发现塔河油田北部于奇地区发育非对称式共轭走滑断裂体系,据此划分出塔河油田南、北两大走滑断裂体系,但两者在形成时间、最大主应力方向、动力学来源等方面均存在明显差异;(2)两大走滑断裂体系的形成及演化过程大致分为3个阶段：加里东期中期I幕形成南部NNE、NNW两组走向的“X”型共轭走滑断裂体系,最大主压应力方向为N3°E,主要受到西昆仑洋俯冲消减作用产生的由S向N的挤压应力;加里东期晚期—海西期早期形成北部以近NE走向为主、少量SN走向组成的非对称式共轭走滑断裂体系,最大主压应力方向为N23°E,主要受到南天山洋俯冲消减作用产生的由N向S的挤压应力,海西期晚期继承性活动,主应力方向保持不变;印支期—燕山期两大走滑断裂体系继承性活动,最大主压应力方向为N37°E,主要受到南部古...     

台湾地区近期7级以上大地震前的中强地震条带现象

本文分析了台湾地区近期 5 5级以上地震的时空图象 ,发现台湾地区从 93年开始5 5级以上地震逐渐形成一条中强地震共轭条带 ,其中主条带长仅 30 0公里左右。强震条带形成后 ,1 999年 9月 2 1日在主条带的西南段 1 /3处条带的西侧边缘南投发生了7 6级地震。南投地震后 ,5 5级以上地震条带仍然存在 ,2 0 0 2年 3月 31日在共轭条带交汇处附近苏沃海外发生了 7 5级地震。在共轭条带内相隔两年半时间里发生两次 7级以上地震实属罕见     

复杂地形边界层流场分析的共轭函数方法

将变分资料同化方法用于复杂地形边界层流场分析，利用兰州的一次观测资料进行了数值试验。结果表明，这样得到的分析场比简单客观分析场或数值模拟场都要合理。