学科分类
/ 2
29 个结果
  • 简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利EFGM计算复杂模型精度高的优势,提出了EFG—PIMEFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIMRPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。

  • 标签: 无单元Galerkin法 点插值法 径向基点插值法 泊松方程 线源二维正演 无网格耦合法
  • 简介:边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原

  • 标签: 地球物理异常 归一化边界识别 归一化总水平导数 正则化倾斜角 theta图
  • 简介:本文分析了海陆两栖地带地震勘探中同时使用的速度检波器和加速度检波器存在的差异及其对地震信号的影响;在此基础上,设计研制出了陆用压电检波器,并对其性能特点进行了分析;通过试验资料分析,消除了海陆两种不同机理的检波器资料的相位差的问题,陆用压电检波器和水中压电检波器记录信号的频带和能量达到一致,提高了地震资料的分辨率;实现了海陆地区可以同时采用相同机理的检波器进行地震信号的接收,解决了滩海地区速度检波器和加速度检波器长期混用的问题。更多还原

  • 标签: 海陆两栖 振动系统 陆用压电检波器 频率特性 频宽
  • 简介:阵列声波信号是典型的非线性、非平稳信号,Hilbert-Huang变换(HHT)是处理非平稳信号的一种比较新的时频分析方法。通过对信号进行经验模态分解(EMD)和对瞬时频率的求解,可以获得声波信号的时-频谱。其关键技术就是进行经验模态分解,任何非平稳的信号都可以分解为有限数目并且具有一定物理意义的固有模态函数。EMD方法可以理解为以声波信号极值特征尺度为度量的时频滤波过程。滤波器充分保留了声波信号本身的非线性和非平稳特征,在声波信号的滤波和去噪中具有很大的优势。文中介绍了HHT时频滤波的实现过程,并列举了一些声波测井波列实例,说明了该方法的有效性。

  • 标签: HHT 经验模态分解 固有模态函数 时频滤波
  • 简介:分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。

  • 标签: 分数阶S变换 FASTICA 分数阶时频分析 谱分解
  • 简介:不同类型的页岩,微观物性特征差异明显,本文针对四川盆地龙马溪组页岩气储层进行岩石物理建模VTI各向异性参数反演。首先,基于前人对粘土矿物的定向排列是产生页岩固有各向异性主要原因这一地质认识,在岩石物理建模过程中引入粘土矿物压实指数CL参数描述粘土矿物的弹性各向异性。之后,基于岩石物理模型开发反演算法,计算页岩储层CL参数Thomsen各向异性参数,解决了由于无法测得与井壁垂直方向上的声波速度,各向异性直接测量存在困难的问题。计算结果表明,通过在岩石物理建模中引入粘土压实参数,反演方法能够合理估计龙马溪页岩储层的弹性各向异性,反映了龙马溪页岩的微观物性特征。进一步分析发现,龙马溪页岩中粘土含量与参数CL相关性较弱,表明粘土矿物的多少对其压实或各向异性程度影响较小。同时,参数CL在目标层龙马溪组底部和五峰组具有高异常值,反映了储层微观结构与含油气特征具有关联性。最后,基于模型构建了岩石物理模板,可用于储层测井数据与多物性参数关系的定量解释。测井数据在岩石物理模板上的合理分布也验证了岩石物理建模方法的有效性。

  • 标签: 龙马溪 页岩 各向异性 岩石物理 压实指数
  • 简介:通过最新的全球地磁模型——CALS10K.1b,结合CALS3K.4与IGRF11模型,计算并分析了10000BC~1990AD期间中国大陆邻近地区非偶极子(ND)磁场Z分量的时空变化。为了深入了解ND场的变化,从场源的角度,对2n(n=2-10)极子ND场及其对应的能量进行了分析。结果显示在研究期间ND场的变化可分为3个阶段。在10000BC-2500BC期间,ND场以正值为主并持续了近7500a,在2500BC-1500AD期间强度转弱为以负值为主并持续至1500AD,自此快速增强为以正值为主。东亚地区ND场异常基本在截断阶数(n)为3时即形成,且该异常区已在大陆地区内形成了封闭的圆形区域,这意味着前3阶的ND场占据了总ND场强度的大部分。ND场在核幔边界(CMB)处衰减较快,在地表处则趋于稳定。

  • 标签: 地磁场 非偶极子 CALS10K.1b IGRF11
  • 简介:分数阶S变换(FRST)具有较强的时频聚集性。利用FRST处理地震数据,通过合适的分数阶参数将频率轴旋转到适当位置,即可实现目标地质特征信息的最佳识别。由于不同的地震信号的最优分数阶参数可能不同,因而对整体的分数阶参数的最优估计不利于对多道地震数据的处理。本文首先利用FRST分离出共频率数据体,并利用共频率数据体进行了低频伴影分析,然后提出FRST和盲分离结合的方法,不需要对地震数据的最优分数阶参数进行估计,即可提取识别有效地质特征信息的独立频谱,提高对地震数据的解释效率。仿真实验表明在分数阶时频域内此方法能有效分离出独立的频率信息。将该方法用于实际的地震数据,并与已知井信息进行比对,验证了其有效性。

  • 标签: 分数阶S变换 FASTICA 分数阶时频分析 谱分解
  • 简介:四川盆地及其周缘地区的筇竹寺组和五峰-龙马溪组页岩是目前国内页岩气勘探的主要层位之一,但其地震弹性性质响应规律的区域性特征需要开展相关的实验和理论研究工作予以明确。本研究对干燥状态下的筇竹寺组和五峰-龙马溪组页岩的露头样品进行了超声波速度测试,系统地分析了地震弹性性质随页岩岩石特征的变化规律。研究结果表明,孔隙度与粘土矿物含量呈正相关、与脆性矿物含量呈负相关;粘土、石英、长石和碳酸盐构成页岩岩石基质,与孔隙共同构成页岩岩石骨架,而干酪根和黄铁矿主要赋存于孔隙中,与页岩骨架的耦合较弱。通过将全部连通性孔隙近似等效于仅存在于骨架粘土矿物之内和采用Gassmann流体替换类似的思路处理干酪根和黄铁矿,可以较为简单地将自相容近似(SCA)理论、微分等效介质模型(DEM)和Gassmann方程组合起来构建研究区页岩的地震岩石物理等效模型。该模型通过关键参数等效孔隙纵横比(采用研究区样品平均值或者由碳酸盐含量进行估算)可以较为准确地预测筇竹寺和五峰-龙马溪组页岩的纵波速度,验证了等效模型的有效性和较广的适用性,可为筇竹寺组和五峰-龙马溪组页岩气储层的测井解释和地震“甜点”预测提供依据。

  • 标签: 筇竹寺组页岩 五峰-龙马溪组页岩 地震岩石物理 等效模型 等效孔隙纵横比