利用重力归一化总梯度及相位法研究断裂构造

用波数域的位场转换和快速傅立叶变换理论研究重力归一化总梯度和归一化相位;在波数域向下延拓和导数计算的滤波算子中,分别引入圆滑滤波因子,抑制了向下延拓计算中对噪声干扰的放大作用,增加了计算的稳定性;利用改进的方法,计算并阐述了4种断裂构造GH场等值线和相位曲线的变化特征;别列兹金重力归一化总梯度法中极大值理论在研究断裂构造中具有局限性,对于无限延伸断裂构造,应根据GH场等值线的走向、形状的特点,结合归一化相位曲线的转折点确定断裂空间位置.     

利用DCT法计算重力归一化总梯度

利用重力归一化总梯度方法可以寻找场源位置,通过“两高夹一低”的典型场特征可以检测油气构造,传统的方法主要采用的是别列兹金法（即傅立叶级数法）和傅立叶变换法。利用DCT法计算重力归一化总梯度是最新研究方法。从理论上给出了基于DCT的重力异常转换公式,进而研究了无限长水平圆柱体、三度体非含油气球冠以及储油球冠模型G_H场特征。模型实验证实：随着地质体埋深的增大,谐波数会逐渐减小,G_H场奇点极大值向下移动,且呈现线性衰减的特征;储油球冠G_H场等值线表现出明显的“两高夹一低”的储油气藏典型标志;用DCT法计算的G_H场具有优良的稳定性,而且能够增加向下延拓的计算深度。     

重力归一化总梯度及其相位的叠置方法研究

通过理论推导和模型实验证明：将重力归一化总梯度及其相位二者相叠置的方法，用于确定壳内深大断裂的空间位置是准确可行的。在实际资料处理地，将该方法应用于新疆横穿塔里木盆地及其周缘造山带的大红柳滩－老虎台条带域，划分出盆地的边界和多条深大断裂。对各种参数的影响和选择做了深入的探讨。     

重力归一化总梯度计算法在苏州西部地区的应用

通过在苏州西部地区地质剖面重力资料计算,笔者认为归一化总梯度法可以定性地提供场源的深度及分布状态,为了解深部构造、寻找隐伏岩体提供重力场的判断依据。     

基于导数迭代法的重力归一化总梯度

针对传统重力归一化总梯度方法对油气藏分辨率不够和深度计算精度较低的问题,本文利用导数迭代法的向下延拓结果,计算重力归一化总梯度。将传统方法和本文方法同时用于模型试验,分析谐波数N、埋深和充满系数对G~H场特征的影响。结果表明,相比于传统方法,本文方法受谐波数N的影响较弱,G~H场对油气藏的分辨率更高,反演的油气藏深度更加准确。通过实际应用,证实了本文方法的有效性。     

利用正则化方法计算重力归一化总梯度

针对重力归一化总梯度计算过程中向下延拓的不稳定性,将正则化因子引入到向下延拓中,提出了利用正则化方法计算重力归一化总梯度的方法。在模型试验中,对非均匀球冠模型中不同厚度油气藏进行试验分析,结果表明:选取合适的谐波数和正则化因子,可以较好识别出"两高夹一低"的储油特征,且识别能力强于常规的Fourier变换法,这说明本文方法能够提高重力归一化总梯度计算结果的稳定性和分辨率。在实例应用中,利用此方法不但清晰地显示了已有油气藏的位置,还推测了5个油气构造区。     

重力归一化总梯度法勘探浅层油气藏浅析

随着仪器采集精度的提高,分析、反演解释技术的长足进步以及油气藏重力勘探的需要,出现了重力归一化总梯度法。结合测井资料,通过重力模型正反演计算,得出高精度重力勘探浅层油气藏是可行的结论。重力归一化总梯度法自提出以来获得了不断的完善发展,并在某气田的实践中被证明是勘探浅层油气藏的有效方法。     

丹青河林场-道台桥剖面的断裂及构造分区

在波数域应用重力归一化总梯度及相位法研究了4种典型的台阶模型,讨论了GH场等值线和相位曲线与台阶断面的对应关系;应用该方法对大庆探区外围丹青河林场—道台桥剖面（简称DB4剖面）野外实测的高精度重力数据进行处理,划分出8条断裂;通过与项目组大地电磁测深的反演结果进行对比,研究了依兰—依通断裂带、牡丹江断裂和勃利盆地的西缘边界断裂;综合重力与大地电磁测深数据,划分了剖面的构造分区并作出地质解释。由于断裂构造与生油凹陷的发育和油气圈闭的形成密切相关,笔者旨在研究DB4剖面的主要断裂和构造特点,从而为在大庆探区外围进一步优选油气区提供地球物理佐证。     

波数域重力归一化总梯度法中的圆滑滤波因子

讨论了波数域重力归一化总梯度法,提出在位场转换的向下延拓和求导过程中分别引入圆滑滤波因子,该方法增强了计算的可靠性和稳定性,扩大了可延拓深度;通过模型实验和数据分析,确定了向下延拓和导数因子中圆滑滤波因子指数参数的取值规律,指出了欲求地质体越浅(深),采用的谐波数越大(小)的地质体质心埋深与谐波数N之间的变化规律;发现了均匀密度三度背斜模型和非均匀密度储油背斜模型与二度模型相似的G^H场等值线特征,为地质体质心埋深计算及储油构造研究提供了思路。     

基于泰勒级数迭代的重力归一化总梯度

针对重力归一化总梯度常规计算方法中向下延拓计算的不稳定性,提出了基于泰勒级数迭代法进行快速稳定向下延拓的重力归一化总梯度场( GH 场) 计算模式。水平圆柱体模型试验表明: 该方法的GH 场表现为单一高值特征,而且极大值位置与模型体质心位置相吻合。非均匀( 含油气) 球冠模型表明: 在油气含量较低时,Fourier 级数法和Fourier 变换法无法识别出“两高夹一低”的储油特征,而本文方法可以清晰地识别。这充分说明基于泰勒级数迭代的重力归一化总梯度计算方法不仅提高了计算结果的稳定性,而且还提高了异常的分辨率。     

三度体球冠模型重力异常计算方法及其重力归一化总梯度场特征

提出利用铅垂柱体薄片作为面元的面元积分法,计算三度体球冠重力异常的计算公式,用三度体球冠模型模拟三度背斜体,计算了均匀密度球冠模型、非均匀密度储油球冠模型的重力异常。最后通过误差分析,验证了该方法的可行性。采用波数域重力归一化总梯度计算方法,计算了均匀密度球冠模型、非均匀密度储油气球冠模型G^H场,发现似三度背斜体与二度背斜体有相似的G^H场等值线特征,似三度贮油气藏背斜构造的G^H场亦表现出明显的“两高夹一低”的储油典型标志。所以,可以用非均匀密度三度体储油球冠模型模拟三度贮油气藏背斜构造。     

重力及其梯度异常正演的Moving-footprint大尺度模型分解方法

重力及其梯度异常正演计算效率决定了反演的可行性,也是高效构建足量、多样深度学习样本数据的基础。为了进一步提高重力及其梯度异常正演的计算速度,受航空电磁领域"Moving-footprint"快速正演技术的启发,本文在基于网格点几何格架函数空间域快速正演的基础上,提出了一种"Moving-footprint"大尺度模型分解的重力及其梯度异常正演计算方法。此方法在地下半空间内选择观测点正下方一定有效范围的子空间,该观测点异常近似为其对应子空间内长方体单元的异常和,而忽略子空间外长方体单元产生的异常;当观测点移动,其对应的子空间随之移动,从而可以将地下半空间长方体模型进行大尺度模型分解,为每一个计算点所对应的子空间进行正演计算。模型实验表明,在256x256x15个长方体模型的地下半空间内选取32x32x15的子空间进行计算,重力异常及部分梯度异常的相对平均误差小于10%,计算速度提高19倍;文中方法针对1024x1024x15个长方体模型计算时间约为32min,相比已有算法中存在的超常规计算量的瓶颈问题具有显著的计算优势。     

移动平台全张量重力梯度数据的噪声抑制

全张量重力梯度仪器测量数据中包含了大量的白噪声和有色噪声。传统的数字滤波器只能滤除某一频段外的噪声,对于混叠在重力梯度有用信号频段范围内的有色噪声不能很好的对其进行分离。为了同时滤除白噪声和有色噪声,笔者利用卡尔曼滤波器采用增广矩阵法将全张量重力梯度数据中的有色噪声进行估计,在抑制白噪声的同时将有用信号和有色噪声分离,并利用数字滤波器与卡尔曼滤波器的优点,将其结合生了更好的滤波效果,得到了更高质量的梯度信号。通过模型试验验证了本方法对噪声的滤波能力,并满足高精度重力梯度数据处理要求。     

基于重力及重力梯度数据的位场偏移成像

为了验证位场偏移快速成像方法的有效性以及对比基于重力和重力梯度数据的偏移成像的效果,笔者利用MATLAB编写了重力及重力梯度数据位场偏移成像的算法,并进行了模型正演计算,发现基于梯度数据的偏移成像的横向分辨率、抗噪性和抗背景场能力好于重力数据的偏移成像,并且能够较好的确定模型的中心埋深。将此方法应用于葫芦岛市黑鱼沟的某铅锌矿洞的深度计算中,验证了该方法的有效性。     

基于重力全张量数据的Parker-Oldenburg算法研究

相比重力数据,重力张量数据通常包含更多的异常信息。本文根据重力数据与重力张量数据的关系,利用位场转化技术,将重力张量数据应用于传统的Parker-Oldenburg密度界面反演算法中。通过模型试验,证明了在网格间距较大或者数据存在一定噪音时,使用本文算法进行反演能得到更好的效果。实验结果说明利用重力张量数据可以有效地提高密度界面反演的分辨率。     

全张量重力数据水平分量噪声滤除方法

为获得更好的噪声滤除方法,本文建立一个球状模型,对全张量重力梯度三个水平分量加入高斯白噪声,并利用FFT和PSD分别对加噪后三个分量(Vxx1,Vxy1,Vyy1)进行频谱分析,根据频谱分析结果分别设计了巴特沃斯低通滤波器和高斯低通滤波器对加噪数据进行滤波,取得了较好效果,并利用复杂模型验证滤波效果。     

全张量重力梯度数据垂直分量噪声滤除方法研究

全张量重力梯度数据中含有高频成分,为了在去除高频噪声的同时,保留浅部小异常体信息,本文利用最小二乘原理,求解重力梯度垂直分量与重力数据联合约束下的超定方程组。结果表明,随着计算窗口的增大,去噪效果越来越好,对于本文所设计的模型,当计算窗口增大到32×32时,重力数据中99%以上的噪声被去除,张量数据Vxz、Vyz中,分别有大约40%、50%的噪声被去除;同时,浅部小异常体信息几乎全部保存。