简介:用基于张量格林函数的体积分方程法对三维异常体进行瞬变电磁响应的正演模拟,首先在频率域内计算电磁场分量的频率域响应,然后利用快速数字滤波技术将计算结果转换到时间域。设计和计算了水平电偶极子源激发下层状水平地层模型背景下的常见地形如山谷、山峰地形的模型,并考察分别把源和接收器放于这些地形中的瞬变电磁场响应,详细分析了这些地形对长偏移距瞬变电磁测深(LOTEM)的影响。结果表明,山谷和山峰地形对LOTEM的结果均有不同程度的影响。当电偶极子源放在山谷谷底时,地形对观测异常场的畸变非常严重;当接收器放在山谷中时,接收器处地形的影响强烈但该影响在空间和时间上只是局部的。总体来讲,不论山峰地形位于何处,其对LOTEM的影响相对较小。当地形处于发射源与接收器之间时,地形对LOTEM的影响非常小,表明在进行LOTEM勘探时,选择发射源的放置比接收器的位置更加重要,野外勘探是尽量把发射源选择在开阔的平坦位置。
简介:裂缝识别对碳酸盐岩储层测井评价至关重要,传统的测井仪器由于探测深度浅(小于3米)而无法对井壁外围(大于3米)裂缝发育情况进行评价,远探测声波测井仪器采用相控阵发射、同时加大源距,有效提升了测井仪器的探测深度。但由于缺少正演模拟研究,目前对于裂缝的解释往往是基于经验而缺乏理论依据,导致很多现象难以解释。本文利用高阶有限差分方法对远探测声反射波测井裂缝识别进行了正演模拟及叠前逆时偏移成像研究,首先在理论研究的基础上构建正演模型,重点研究不同裂缝离井壁距离、裂缝张开度和倾角的响应特征;其次在单因素变化基础上提取成像区域的能量强度,分析确定出在实际地层速度有波动变化时远探测声反射波测井方法能够识别裂缝的条件;最后通过对裂缝识别的影响因素定量化分析,确定裂缝识别的最大距离、最小张开度和最小倾角,降低了裂缝识别与评价中的多解性。研究成果对远探测声反射波测井仪器的发展、数据处理方法的改进,以及后续的测井解释工作都有一定的借鉴意义。
简介:ThevariationsintheEarth’srotationareimportanttospacedynamictheoryandnaturaldisastersbecausetheyaffectthelength-of-day(LOD)andconsequentlyhumanlives.Weusemaximumentropymethod(MEM),Lombmethod(LOMB),andphasedispersionminimization(PDM)todeterminethenaturalperiodsofequallyspacedLODtimeseries.WetransformtheobservationalmonthlyLODtimeseries(LODM)tounequallysampledseries(LODMD)byremovingeveryfourth,third,andhalfofthetotalsamples.WealsoapplysplineinterpolationtoLODMDtoyieldequallyspacedtimeseries(LODMDN).Theresultssuggestthatregardlessofthetimeseries,theMEMfrequencyis0.1660month^-1and0.0840month^-1,whereasLOMBandPDMyield0.166month^-1and0.083month^-1,respectively.Furthermore,missingdatathatarelessthanhalfofthetotaldataorsplineinterpolationdonotaffecttheanalysis.Fortheamplitude,neithermissingdatanorsplineinterpolationaffecttheanalysis.
简介:在复杂地表条件的区域,地震数据的采集和处理是一项极大的挑战。虽然可以通过静校正来消除起伏地表的影响,然而当地表高程以及近地表速度剧烈变化时,简单的垂直时移对地震波场造成的畸变会严重降低偏移成像的质量。基于射线的偏移方法可以直接在起伏地表面进行波场的延拓成像,是解决上述问题的有效手段。本文针对复杂地表条件下的高斯波束叠前深度偏移进行研究,对倾斜叠加公式进行修改,使之包含地表高程以及速度的信息,通过直接在复杂地表面进行平面波的合成,得到了一种具有更高成像精度的改进方法。首先简单介绍常规高斯波束偏移的基本原理和计算流程,并以此为基础,给出复杂地表条件下高斯波束偏移原有的实现方法以及本文的改进方法,最后通过模型和实际资料的试算验证本文方法的有效性。
简介:有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。
简介:复杂地层曲面构建可广泛应用于石油勘探、地质建模、地质构造分析等领域中,也是这些领域数据可视化和可视分析的重要基础。现有方法在多种类型断层处理、曲面光滑性等方面存在不足。为解决地层数据分布不规则性带来的曲面光滑性问题,本文引入几何偏微分方程的曲面造型方法开展地层曲面构建;针对多种类型复杂断层存在情况下的曲面构建问题,提出三维空间与二维平面的投影互换算法,给出了基于几何偏微分方程的多种形式断层的统一处理方法;针对复杂地层曲面构建问题,建立了相应几何偏微分方程,给出演化求解算法。实际地层数据的空间曲面构建验证了本文方法计算效率高、能处理不规则的数据分布,同时能对存在断层尤其是逆掩断层的地层面进行曲面恢复。
简介:初至拾取是起伏地表地震资料处理的关键环节,直接影响近地表建模的精度及静校正效果,但起伏地表资料极低的信噪比使传统的自动初至拾取算法几乎失效,而手动拾取要耗费巨大的人工和时间成本。本文研究了改进超虚干涉法,结合多道多域初至质量监控技术,实现起伏地表资料初至波自动拾取。改进超虚干涉法首次将近地表散射波纳入干涉法提高信噪比的范围,并通过折射波和散射波的线性组合使干涉法适用于起伏地表条件下任意复杂射线路径的初至波类型;提出反向干涉和多域干涉的概念,显著增强了所估计的虚源信号;使用波形反褶积滤波器较好地抑制了干涉导致的"假事件"的形成;采用多道多域初至质量监控技术,实现错误"假事件"初至的自动归位,提高了拾取初至的稳定性。本文研究的初至波自动拾取理论与技术具有突出的鲁棒性和稳定批量处理大量实际三维地震数据的能力,在中国西部山区三维地震勘探数据处理的应用效果显著,质量优于某常用商业软件。
简介:频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和Hessian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。
简介:地震波场数值模拟中不可避免地会出现边界反射,一般采用吸收边界条件以压制人工边界反射。目前常用的分裂式完全匹配层(PML)边界条件需要在边界处进行特殊处理,尤其是在三维情况下需要将变量分裂为三个分量,增加了数值模拟的计算时间和内存占用量。与分裂式PML吸收边界条件相比,混合吸收边界条件(HABC)具有易于实现、计算量小和吸收效果好等优点,可以提高三维波动方程数值模拟的计算效率。本文将基于一阶Higdon单程波方程的混合吸收边界条件从二维计算域发展到三维,提出了适用于三维弹性波数值模拟的混合吸收边界条件。均匀模型以及复杂模型的三维数值模拟结果表明,混合吸收边界条件与传统的完全匹配层边界条件相比,具有效率高、吸收效果好的优势。
简介:基于双相介质理论的AVO正演技术是储层性质描述和流体预测的有效技术手段之一,但是输入参数中基质矿物模量的准确性和双相介质模型的的合理性极大地影响双相介质AVO正演效果。因此,本文采用基于流体因子的基质矿物模量反演方法,自适应反演基质矿物体积模量。引入具有岩石物理意义的多约束条件,使得流体替换技术制作的双相介质模型具有岩石物理意义。保证获得的双相介质AVO特征反映实际地层响应,真实可靠。通过不同岩性岩样的对比分析,说明反演方法的优越性和准确性。同时LH地区实际资料应用,获得孔隙度和流体饱和度等重要岩性参数变化时双相介质AVO特征,特别是不同储层孔隙度在同一入射角对应快纵波和横波反射系数幅值的大小差异和突变角差异是分辨储层孔隙度大小的依据。
简介:High-qualityseismicgeometryisthekeytoobtainhigh-qualityseismicdata,andcanaffecttheaccuracyofdataprocessingandimaging.Basedontheanalysisoftherelationshipbetweenthequalityofthegeometryandthefouracquisitionparameters(thenumberoftraces,shotlinespacing,andthespaceandnumberofreceiverlines),aqualityevaluationmethodofthegeometrybasedoncomprehensivequalityfactor(CQF)isproposed,andtherelationshipbetweenthegeometryqualityandthefourparametersisgiven.WeusefielddatacollectedinanoilfieldinWesternChinawithcomplexgeology:Firstweuseawideazimuthgeometry.Then,wecalculatetherelationshipcurvebetweengeometryanddataqualitybyvaryingeachparameterwhilekeepingtherestfixed.andtheanalysisresultsaregivenbyusingtheCQFevaluationmethod.Theresultsshowthattheshot-linespacinghasthegreatesteffectonthequalityofthegeometry,andtheincreaseofthereceiverlinespacingcanappropriatelyimprovethequalityofthegeometry,andtheincreaseofthenumberofreceivingtracescanimprovethegeometryquality.Thedifferentacquisitionparametershavedifferenteffectsontheimagingqualityofshallowanddeepevents.Themodelforwardandprestackdepthmigrationareusedtogenerateprestackdepthmigrationprofileswithdifferentacquisitionparameters.Theimagingresultsareconsistentwiththeabovecalculatedresults.Accordingtothedepthofthetargetlayer,thequalityfactorevaluationmethodisappliedtoguidethedesignofthegeometryandoptimizetheacquisitionparameterstoimprovetheimagingaccuracyofseismicdata.
简介:为了解决海底起伏变化对地震波场的影响问题,本文提出将(x-z)域中的曲网格映射成(ξ,η)域中的矩形网格,推导出(ξ,η)域中的二维标量声波方程,根据推导出来的波动方程采用逆时有限差分法将海面上采集到的地震波场在(ξ,η)域中向下延拓至海底面,延拓时采用海水的速度,然后采用顺时有限差分法将延拓后的地震波场再反延拓到海面上,延拓时采用海底面以下地层的速度,从而消除了海底起伏带来的负面影响。模型及实际地震资料的计算分析表明该方法不但能够校正由于海底起伏所引起的海底面下地层反射波场的不连续性还能够校正由于海底起伏所引起的地震波的动力学特征的变化。对延拓前后的地震波场进行速度反演,延拓后反演的地层速度比延拓前反演的地层速度的精度提高很多,延拓前后地震波场的叠加剖面对比表明该延拓方法能够明显提高地震波场的成像质量。
简介:在多个区块的测井评价工作发现,孔隙结构直接影响储集层的品质和油气层的电阻率,是测井准确评价流体性质的关键。岩石物理资料表明不同区块内影响储层孔隙结构的因素不同,但效果是一致的,即孔隙结构的复杂程度控制着储层的储集能力和渗透能力。孔隙结构的复杂程度影响了储层中导电流体的分布和含量,从而控制了储层的电阻率。储层出现低阻油气层的内因均为复杂的孔隙结构(骨架导电及工程原因除外)。测井储层评价在分析控制储层孔隙结构复杂程度的地质因素及储层分类的基础上,针对不同类型储层采用不同的模型、参数和标准,可以有效的认识储层品质和识别不同类型储层的流体性质。
简介:在水库现场试验了RISK2型探地雷达探测水库冰厚度的能力,试验时所用天线频率为600MHz;同步钻孔测量雷达探测处的冰厚度;以及在一个点上取样测试分析冰晶体、冰内气泡和冰密度。试验时冰面积雪厚度0.03-0.05m,冰层上部有0.24m粒状冰,其下均为柱状冰;冰内气泡含量呈表层高底层低分布;冰密度随气泡含量变化;冰厚度在平面内不均一。通过探测厚度和实测厚度的对比分析以及气泡含量对介电系数影响的理论分析,建立了积雪、粒状冰和柱状冰三层介质模型,获取雷达波在冰内的理论传递时间。结果发现:能够利用等效介电常数或等效传播速度评价雷达波传递时间,结冰期冰层1/3深度处的对应介电常数或传递速度可以作为等效值;另外因冰内大气泡造成的理论传递时间大于雷达探测时间,其差值随理论传递时间或冰厚的增加呈非线性增加。
简介:Tikhonov正则化(TR)方法在重磁数据处理中发挥了重要的作用,本文在研究如何利用Tikhonov正则化方法方法解决重力数据3D反演的同时,深入讨论了可进一步提高拟合误差的ExtrapolationTikhonov正则化方法(EXTR)的原理,并就其参数选择方法及各参数对拟合误差、迭代次数及反演结果的影响进行研究。常密度及变密度组合模型试算结果表明,与TR方法相比,EXTR方法不仅可以达到解释人员设定的先验拟合误差水平,在计算时间及迭代次数相应增加的前提下有更高的拟合精度;同时其反演结果也更加紧致,进一步改善了TR反演结果的发散性;并且其反演数据范围更贴近预设模型参数范围,模型特征与预设模型密度分布吻合较好。更多还原