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

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

基于改进的预处理共轭梯度法和三维插值法精确和快速解算GRACE地球重力场

本文首先对比分析了基于直接最小二乘法(DLSP),预处理共轭梯度法(PCCG)和三维插值法(TDIM)解算卫星观测方程的优缺点;其次,基于能量守恒法分别利用改进的PCCG和TDIM反演了120阶GRACE地球重力场,在120阶处累计大地水准面精度为17.316 cm和15.421 cm;最后,通过和德国波兹坦地学研究中...     

井地大地电磁非线性共轭梯度二维反演研究

本文提出并验证了地下大地电磁测深数据的正、反演算法和程序,为井地大地电磁测深和井地可控源音频大地电磁测深这一前沿研究课题提供了地球物理解释方法.使其更加完善和有效.在Rodi(2001)发表的大地电磁场二维NLCG反演理论和算法的基础上,推导了井地大地电磁场二维NLCG反演算法,改编反演的代码,开发出二维NLCG反演程...     

新疆巴楚地区共轭膝折带的物理模拟研究

本文以塔里木地区巴楚隆起地层及构造为原型,设计了3层结构模型:底部为黄油层,作为滑脱层;中间为黏土层,模拟该地区膝折发育的地层;上部为细砂层,模拟上部松散沉积.研究了挤压速率、黄油层厚度、砂层厚度等因素对共轭膝折带形成的影响,对共轭膝折带的形成和演化阶段进行分析.结果表明:形成膝折带的有利条件是较厚的基底滑脱单元、上覆...     

基于正则化条件的地震数据局部信噪比估计方法

信噪比是衡量地震数据质量的重要指标之一,在地震数据处理和解释中有着重要的作用.目前已有的地震数据信噪比估计方法往往得到的是整个数据的全局信噪比,这种方法只能说明地震数据总体质量的好坏,无法直观细致地刻画地震信号的局部质量.本文提出一种基于正则化条件的局部信噪比估计方法.该方法的基本原理是使用正则化共轭梯度法求解局部信噪比最优解,正则化算子的参数将控制地震信号各点数据局部信噪比的平滑性.其中应用一种基于“过滤波”的级联信号估计方法来计算有效信号,该方法利用有效信号和噪声的相关性特征计算局部信噪比中的有效信号.局部信噪比估计方法利用了信号中每个数据点及其邻域各点的局部信息,避免了使用单个数据点而可能出现的信噪比不合理值,而且局部处理能够减少全局噪声对信噪比估算的影响,该方法可以更准确地表征地震资料信噪分布特征.另外,局部信噪比对去噪方法的评估也具有重要意义.理论模型测试和实际资料处理结果表明,局部信噪比估计方法能够准确反映任一给定地震信号剖面的局部信噪比特征,为非平稳地震数据质量评估提供了直观的评判标准.

     

基于不等式约束的最小二乘法三维电阻率反演及其算法优化

基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景.     

基于耦合标势与矢势的有限体积法模拟非均匀各向异性地层中多分量感应测井三维响应

为快速有效地研究、考察各向异性地层条件下多分量感应测井的响应特征,本文利用电场标势与矢势的有限体积法研制出三维频率域电磁场响应的数值模拟算法,克服由低频发射或高阻地层产生的低感应数问题,有效提高了三维电磁数值模拟算法的应用范围和计算效率.首先利用电场的标势与矢势将Maxwell方程转化为满足库仑规范条件的耦合势Helmholtz方程,以Yee氏交错非均匀网格中不同位置上的节点为中心建立四种控制体积单元,通过对控制体积单元中电磁场与电导率的积分平均实现耦合势方程和磁偶极子旋度的离散,并得到一个对角占优的大型稀疏复线性代数方程组,然后,通过不完全LU分解预处理和稳定双共轭梯度法快速求解离散方程.数值结果证明了该算法的有效性,并进一步考查了仪器偏心、倾斜井、垂直裂缝等复杂条件下多分量感应的响应特征.     

地面可控源频率测深三维非线性共轭梯度反演

讨论了地面可控源电磁勘探三维非线性共轭梯度反演的可行性以及反演过程中考虑场源的必要性.反演采用非线性共轭梯度反演方法.反演过程中,模型响应利用交错网格有限差分技术计算.反演数据采用与发射源平行的电场x分量Ex.利用层状导电模型作为背景,设计了两个理论模型进行数值试验:第一个模型中包含两个电阻率异常,以检验反演的有效性;第二个模型中,在测区外设置了一个低阻异常,以考察源的信息在反演中的作用.两个模型的反演分别从层状背景模型开始,迭代120次后终止.数值试验结果表明,(1)非线性共轭梯度反演所获得的电阻率分布和理论模型吻合较好;(2)非线性共轭梯度算法收敛速度较慢,需要较多的迭代次数完成反演;(3)对于可控源频率电磁勘探,必须考虑源位置信息.因此,本文采用考虑场源信息的地面可控源非线性共轭梯度反演方法能完成真正意义上的可控源频率电磁测深数据的反演.     

南海西南次海盆与南沙地块的OBS探测和地壳结构

跨越南海西南次海盆南部陆缘和南沙地块中部的OBS973-1测线是南海南部首次采集的海底地震仪(OBS)广角反射与折射深地震测线,本文通过震相分析和走时正演拟合,获得了沿测线的二维纵波速度结构模型.模拟结果显示表层沉积物速度2.5~4.5 km/s,厚度1000~3000 m,局部基底面起伏较大.结晶基底的速度从顶部的4.5~5.5 km/s增加到地壳底部的6.8~6.9 km/s,中地壳有一个小的速度不连续面(0.1~0.2 km/s),而地壳底部的莫霍面有较大的速度反差(1.2 km/s),上地幔顶部的速度为8.0~8.1 km/s.莫霍面埋深和地壳厚度在测线的北段和南段有很大的不同,在测线北段的海盆区,莫霍面埋深约11 km,结晶地壳的厚度仅为5~6 km,表现为典型洋壳的特征,而在测线南段的陆块区,莫霍面埋深最大达24 km,地壳厚度可达20 km,表现为减薄陆壳的特征,从海盆区到陆块区莫霍面埋深和地壳厚度迅速增加.陆块区上下地壳的厚度和变化趋势相似,下地壳没有看到高速层(HVL),可能说明地壳内部是以纯剪拉张的均匀减薄为主,地壳下部的岩浆底侵不发育.对比OBS973-2和OBS973-3测线的结构模型,可以推测南沙地块的中部和东部具有相似的构造性质,西南次海盆两侧是一对非火山型的不对称共轭陆缘.