简介：有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层（PerfectMatchedLayer,PML）吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的＂反射波＂,以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。

简介：海底电缆双检采集是利用水检和陆检对海水层鸣震的响应极性相反,通过水、陆检数据的合并处理来消除海水层的鸣震,进而拓展了数据的频带范围,提高了资料的保真性和分辨率,而得到了广泛应用。但海底电缆双检资料中的陆检数据,由于受洋流、海底地形和施工工艺等的影响,很难保证其与海底的良好耦合,因而受耦合系统传输函数的影响比较大,造成水、陆检数据在振幅、频率、相位等方面的差异,降低了陆检数据的信噪比,影响了双检地震数据的合并处理效果。而海底电缆双检中的水检资料检测的是海水压力的变化,与海底不存在耦合问题,因而水检数据信噪比往往比较高。本文首先给出了陆检与海底耦合系统传输函数的数学表达式,然后利用水检数据作为约束条件估算出了陆检与海底耦合系统的传输函数,并利用估算出的传输函数对陆检数据进行了耦合校正处理,解决了陆检与海底的耦合效应对陆检数据的振幅、相位等的影响,提高了陆检数据的信噪比,取得了较好的应用效果。

简介：提出一种对整幅高光谱图像的稀疏表示结果进行直接显示的方法,图中不仅包含了稀疏表示中保留的光谱信息,还可显示整幅图像的空间信息。稀疏表示后,将字典中的各有效原子根据光谱特性选择颜色标签,之后根据稀疏系数进行混合颜色显示,此时的图像能够同时满足可分性及距离保持特性。针对局部地物时,提出的单像素混合阵列表示法及改进的裂片纹理技术能够直观且完整的显示出每个像元的具体组成情况,还能够根据所生成图像中的信息对原始HSI进行重建,进而提高数据的利用率。该模型不仅能够良好地显示地物的空间特性,同时能够显示稀疏系数的组成,同时单像素混合阵列表示法及裂片纹理技术弥补了混合像素彩色显示中颜色表达混乱的弊端。对真实地物数据进行实验,结果证明该模型产生的彩色图像具有良好的视觉效果及可分性,满足距离保持特性。

简介：地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层（PML）边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件（HABC）具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。

简介：地层品质因子Q反映了地层对地震波吸收衰减的强弱,是油气检测的重要指示因子,准确提取Q值十分重要。本文根据地震子波能量衰减特征,提出基于子波衰减前后能量比的Q值提取方法。该方法是综合利用多点信息估算子波能量,规避了常规Q值方法（如质心频移法）对地震子波谱的假设,适用于任意波谱形状,但对资料的保真度要求较高。VSP正演模型测试发现,在非固有衰减（如几何扩散、反射透射损失等）得到准确补偿的前提下,能比法能较好地反演地层Q值。实际零偏VSP资料应用表明,能比法提取的Q值趋势与谱比法基本一致,具有一定的可信度。

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

简介：埋藏在6000m以下的地震反射异常体是礁体还是火成岩体?或者是海底滑塌体?在高质量重磁电资采集的基础上,通过对研究区域内的地层岩石物性细致研究,利用CEMP资料界面约束二维反演、重磁力异常剥离技术和LCT重磁力地震联合反演解释技术细致研究了异常体的物性参数,从密度、磁化率、电阻率、速度等多种特征推测判断了深层异常体地质属性,为钻井部署提供了间接依据.

简介：震源定位是微地震监测关键技术之一。本文提出用于微地震定位的弹性波和加权弹性波（WEW）干涉成像方法。该方法在保证定位精度的同时，还可避免震源假象。通过各向同性水平层介质状模型的数值试验，初步表明该方法可适应低信噪比微震信号、速度随机扰动、较稀疏的检波器分布等情况，并在速度模型存在一定的系统误差时也仍保持较高的定位精度。由于干涉成像方法不需要进行初至拾取，定位效率相对传统走时方法也得到了提高。采用二维断层模型试算Nnumericalresultsofusingatwo-dimensionalfaultmodel，表明方法还能实现多震源定位，且比逆时成像有更高的定位精度。

简介：相干体是一种能够有效反映断层和裂缝等地质上非连续性的地震资料解释工具。然而由于受到附近地层,尤其是具有较强反射系数的地层的影响,微小断层和裂缝在相干体的时间切片和纵剖面上仍然难以辨别。本文提出了一种基于局部直方图规定化的相干体增强方法,实现相干体中的微小断层和裂缝的增强。与传统直方图规定化方法不同,本文方法处理三维相干体时不需要对数据进行离散化以统计其分布情况。在将相干体划分为若干子块后,以相干体整体的分布函数作为每个子块内的目标分布函数进行直方图规定化。另外,相邻子块的部分区域需要重叠,以克服分块产生的边界效应。对实际数据的处理结果表明,该方法可以提供更多的微小断层和裂缝的细节信息。

简介：本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.

简介：常规的重力梯度法对地质体位置划分很模糊,当深度越大,地质体越小时,对Vzz和Vzx求导的误差也越大,导致结果与实际偏差很大。本文提出状态判别因子对角点位置进行优化确定,比常规高阶重力导数法更精确,使重力梯度法在对小地质体和断裂上有更好的分辨率,理论模型试算和实际资料处理效果很好,为找矿或划分地下地质体产状提供了更有利的依据。

简介：随着煤矿的深部开采，多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁（简称TEM）法因场结构与地层的耦合关系，对低阻层探测有利。但低阻层的屏蔽作用不仅使得探测同样的深度需要更长的观测时间，而且还会减弱下伏地层的异常响应。本文通过直接时域数值模拟和水平分层大地的模型正演，估算了探测目的层所需要的时间长度，根据噪声对观测数据造成的影响，给出了多层积水采空区可分辨的标准。山西大同达子沟煤矿水文勘探中获得的TEM实测曲线，表现了对多层积水采空区的探测能力。在实测曲线的定量反演解释中，利用电测井资料作为初始参数进行定量反演，约束了煤系薄层的等值性。所推断的三层和两层积水采空区，为钻孔所证实。研究结果表明，当观测时间有足够的长度、下伏地层的异常显示大于观测误差，中心回线TEM法探测多层采空区积水是可行性的。

简介：Theattenuationfactororqualityfactor(Q-factororQ)hasbeenusedtomeasuretheenergyattenuationofseismicwavespropagatinginundergroundmedia.ManymethodsareusedtoestimatetheQ-factor.WeproposeamethodtocalculatetheQ-factorbasedontheprestackQ-factorinversionandthegeneralizedS-transform.TheproposedmethodspecifiesastandardprimarywaveletandcalculatesthecumulativeQ-factors;then,itfindstheinterlaminarQ-factorsusingtherelationbetweenQandoffset(QVO)andtheDixformula.TheproposedmethodisalternativetomethodsthatcalculateinterlaminarQ-factorsafterhorizonpicking.Becausethefrequencyspectrumofeachhorizoncanbeextractedcontinuouslyona2Dtime–frequencyspectrum,themethodiscalledthecontinuousspectralratioslope(CSRS)method.ComparedwiththeotherQ-inversionmethods,themethodoffersnearlyeffortlesscomputationsandstability,andhasmathematicalandphysicalsignificance.WeusenumericalmodelingtoverifythefeasibilityofthemethodandapplyittorealdatafromanoilfieldinAhdeb,Iraq.TheresultssuggestthattheresolutionandspatialstabilityoftheQ-profileareoptimalandcontainabundantinterlaminarinformationthatisextremelyhelpfulinmakinglithologyandfluidpredictions.

简介：为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区（近区、过渡区和远区）资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。

简介：基于岩石物理模型和广义各向异性Zoeppritz方法在频率域计算裂缝型反射体模型反射波方位地震AVO响应。反射体模型为粘弹各向异性、有限厚度的地质体，其地震反射波形序列包含如下动力学信息，即分界面处介质的波阻抗和非弹性差异、反射体内部波的各向异性传播、在传播路径上的频散与衰减，以及来自顶底界面的反射波的调谐与干涉等。计算表明，速度频散和衰减增顶界面反射波大入射角反射时的振幅，而减弱底界面反射振幅。对于固定人射角的方位地震方位地震响应，PP波反射特征表现为随方位角的增加反射波形序列延续时间变长，而PSV和PSH转换类型反射波的方位各向异性变化特征稳定且受储层厚度影响较小，表现为PSV波反射振幅随方位角增加而增加，PSH波在0°和90°方位无反射能量，在45°方位反射振幅最强。

简介：在直流电法有限元数值模拟中，针对常规结构化网格源点附近网格节点数值精度低的问题，并考虑计算效率和反演成像的需求，本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度，降低源点奇异性的影响范围，并提升对地形的模拟能力；通过在计算区域垂向方向上将多个网格收缩到一个网格中，降低网格节点的规模，进而提高数值计算的效率。理论模型检验表明，网格加密-收缩方法能有效地解决源的奇异性问题，与均匀加密网格相比，网格节点总数降低了约80％。

简介：本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。

简介：受沉积环境影响,海底地层层理发育,各向异性介质更接近于真实海洋地质条件。对海洋各向异性研究可更好认识海底构造特征,有效地进行海底资源勘查。本文从各向异性麦克斯韦电磁方程出发,根据电磁场的无源特征引入标量位函数;利用电场和磁场的连续性分别向海底深部和海水中延拓,并将延拓后的电磁位函数在海底耦合到发射源上,从而实现海底电场和磁场的递推求解。我们首先探索如何利用海洋电阻率法识别和求解海底电各向异性特征。研究发现海底各向异性可从视电阻率测深曲线和海底视电阻率极性图进行求解。进而,我们通过海底各向异性地层中高阻体(油气藏)模型进行正演模拟,发现海洋电法视电阻率在浅水区对地下高阻薄层有明显的异常反应。相比海洋可控源电磁法,海洋直流电法不受空气波影响,在浅海区油气资源勘查有着较好的优势,而前人大多基于各向同性模型进行研究,本文实现一维海洋直流电阻率法各向异性正演模拟,算法计算精度高,能很好的为二维、三维正演模拟提供理论参考。更多还原

简介：标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。