简介:Tilt—depth法可以用于快速反演磁源的上顶埋深,但其无法反演磁源底部深度,同时该方法的反演解过于单一。针对于此,本文在有限厚度台阶磁异常导数公式基础上,推导出了可以同时反演地质体上、下界面埋深的改进型Tilt—depth法计算公式,并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt—depth法可以有效地反演出磁性体的上顶与下底深度,且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tilt—depth法用于松辽盆地长岭断陷航磁数据反演,其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致,证实改进型Tilt—depth法相对于常规方法的反演结果更加准确。
简介:无网格法形函数构造不依赖预定义的单元,具有计算精度高、处理复杂模型便利等优点。本文介绍了无单元Galerkin法(EFGM)、点插值法(PIM)与径向基点插值法(RPIM)三种全域弱式无网格法的近似原理及特点;以二维泊松方程为例研究了支持域无量纲尺寸、场节点与背景网格设置对无网格法计算精度的影响。将RPIM与EFGM应用于频率域线源二维正演,给出了RPIM形状参数的推荐值;分析了均匀介质模型大地电磁(MT)二维正演无网格法边界条件直接加载与罚函数法加载的精度差异,结合PIM与RPIM边界条件加载便利及EFGM计算复杂模型精度高的优势,提出了EFG—PIM及EFG.RPIM耦合算法,数值计算结果验证了耦合算法的有效性。研究发现:无网格法及其耦合方法适用于电磁法数值模拟;支持域无量纲尺寸取1.0时无网格法精度与效率高,场节点与背景网格重合时计算效果佳;泊松方程求解PIM及RPIM精度较EFGM低,计算均匀介质MT响应精度较EFGM高;RPIM改善了PIM计算涉及的奇异性问题,对应支持域无量纲尺寸选择空间大。
简介:震源定位是微地震监测关键技术之一。本文提出用于微地震定位的弹性波和加权弹性波(WEW)干涉成像方法。该方法在保证定位精度的同时,还可避免震源假象。通过各向同性水平层介质状模型的数值试验,初步表明该方法可适应低信噪比微震信号、速度随机扰动、较稀疏的检波器分布等情况,并在速度模型存在一定的系统误差时也仍保持较高的定位精度。由于干涉成像方法不需要进行初至拾取,定位效率相对传统走时方法也得到了提高。采用二维断层模型试算Nnumericalresultsofusingatwo-dimensionalfaultmodel,表明方法还能实现多震源定位,且比逆时成像有更高的定位精度。
简介:在文献中,Kirchhoff型反偏移场的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像场信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像场的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移场进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。
简介:本文介绍了在B隧道勘查中应用高分辨电磁成像技术探测隐伏构造的一个成功实例.数据采集采用高频电磁成像系统(STRATAGEMEH4,频率范围从1Hz到90kHz).野外采集正交的电磁场分量,从地磁成像剖面中提取目标体有关电磁学信息.为获取高质量的采集数据,野外采用包含天然场源和全张量可控源的混合场源.B隧道坐落在位于中国中部的湖北省的西部,隧道埋深不到200m,但区域地质作用导致它的地质机构十分复杂.第一次勘查过程中,勘探人员误把一个脱落体的露头当成基岩面,施工过程中出现了冒顶.第二次勘探时采用高分辨电磁法和折射地震法,这次勘探找到了隐伏的基岩面和一个隐伏断裂.勘探结果与后来隧道挖掘揭示的构造吻合.
简介:塔里木盆地碳酸盐岩岩溶型储层广泛分布于塔北、塔中、巴楚等地区,是塔里木盆地碳酸盐岩勘探的重点,也是塔里木油田增储上产的主要层系之一.由于埋藏深度大、反射信号饺酢⑿旁氡冉系?更由于岩溶型储层非均质性严重,使得储层预测的难度增大.通过多年的攻关研究,从碳酸盐岩风化壳储层的地球物理响应特征出发,逐步探索并建立了技术、岩溶型储层测井评价技术、岩溶型储层地震描述技术、应变量分析技术为手段的溶型储层综合评价技术.在塔北、以地质概念模型为基础,以基于层序格架的潜山古地貌恢复塔中等地区应用中取得了良好的地质效果,相继获得了工业油气流.本文介绍了该技术在塔北轮古西地区的具体实施和应用效果.
简介:为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原
简介:随着煤矿的深部开采,多层积水采空区探测成为煤田水文物探工作的新内容。中心回线瞬变电磁(简称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,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:在直流电法有限元数值模拟中,针对常规结构化网格源点附近网格节点数值精度低的问题,并考虑计算效率和反演成像的需求,本文在二维点源问题中提出一种新的网格加密一收缩方法。其核心思想是在结构化网格计算区域中先后引入网格节点加密和网格节点收缩两个环节。通过在计算区域水平方向上均匀加密网格节点密度,降低源点奇异性的影响范围,并提升对地形的模拟能力;通过在计算区域垂向方向上将多个网格收缩到一个网格中,降低网格节点的规模,进而提高数值计算的效率。理论模型检验表明,网格加密-收缩方法能有效地解决源的奇异性问题,与均匀加密网格相比,网格节点总数降低了约80%。
简介:本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。
简介:受沉积环境影响,海底地层层理发育,各向异性介质更接近于真实海洋地质条件。对海洋各向异性研究可更好认识海底构造特征,有效地进行海底资源勘查。本文从各向异性麦克斯韦电磁方程出发,根据电磁场的无源特征引入标量位函数;利用电场和磁场的连续性分别向海底深部和海水中延拓,并将延拓后的电磁位函数在海底耦合到发射源上,从而实现海底电场和磁场的递推求解。我们首先探索如何利用海洋电阻率法识别和求解海底电各向异性特征。研究发现海底各向异性可从视电阻率测深曲线和海底视电阻率极性图进行求解。进而,我们通过海底各向异性地层中高阻体(油气藏)模型进行正演模拟,发现海洋电法视电阻率在浅水区对地下高阻薄层有明显的异常反应。相比海洋可控源电磁法,海洋直流电法不受空气波影响,在浅海区油气资源勘查有着较好的优势,而前人大多基于各向同性模型进行研究,本文实现一维海洋直流电阻率法各向异性正演模拟,算法计算精度高,能很好的为二维、三维正演模拟提供理论参考。更多还原
简介:标量CSAMT只适合一维及测量方向与构造方向垂直的二维情况,对于复杂的三维地电结构,CSAMT需采用张量测量。本文试图采用矢量有限元法实现三维张量CSAMT的正演模拟。为了验证算法的正确性,本文在层状介质中计算了三维CSAMT远区的电场,磁场及阻抗张量,并且与层状介质中的理论解进行了比较,接着还模拟了均匀半空间中含有三维异常体的模型,并且分析了四个阻抗张量、视电阻率及阻抗相位的响应特征。得出如下结论:采用矢量有限元法来模拟三维张量CSAMT,其电磁场及阻抗张量的实虚部计算精度都比较高,并且该方法本身满足电场法向不连续,不用进行散度校正。
简介:MostcurrentprestackAVAjointinversionmethodsarebasedontheexactZoeppritzequationanditsvariousapproximations.However,theseequationsonlyreflecttherelationbetweenreflectioncoefficients,incidenceangles,andelasticparametersoneithersideoftheinterface,whichmeansthatwave-propagationeffects,suchassphericalspreading,attenuation,transmissionloss,multiples,andeventmismatchingofP-andS-waves,arenotconsideredandcannotaccuratelydescribethetruepropagationcharacteristicsofseismicwaves.ConventionalAVAinversionmethodsrequirethatthesewave-propagationeffectshavebeenfullycorrectedorattenuatedbeforeinversionbuttheserequirementscanhardlybesatisfiedinpractice.Usingaone-dimensional(1D)earthmodel,thereflectivitymethodcansimulatethefullwavefieldresponseofseismicwaves.Therefore,weproposeanonlinearmulticomponentprestackAVAjointinversionmethodbasedonthevectorizedreflectivitymethod,whichusesafastnondominatedsortinggeneticalgorithm(NSGAII)tooptimizethenonlinearmultiobjectivefunctiontoestimatemultipleparameters,suchasP-wavevelocity,S-wavevelocity,anddensity.Thisapproachisrobustbecauseitcansimultaneouslycopewithmorethanoneobjectivefunctionwithoutintroducingweightcoefficients.Modeltestsprovetheeffectivenessoftheproposedinversionmethod.Basedontheinversionresults,wefindthatthenonlinearprestackAVAjointinversionusingthereflectivitymethodyieldsmoreaccurateinversionresultsthantheinversionbyusingtheexactZoeppritzequationwhenthewave-propagationeffectsoftransmissionlossandinternalmultiplesarenotcompletelycorrected.
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。
简介:火山沉积型硼矿是一种重要的硼矿床成因类型,主要发育在第三纪湖相蒸发岩沉积层和火山岩互层的二元结构单元中。针对硼矿具有低密度的特点,本文对西藏雄巴地区进行了1:5万高精度重力勘探,获得了研究区布格重力异常和剩余重力异常。根据其异常特征进一步分析了区域构造特征、断裂体系、沉积单元以及火成岩的分布与火山沉积型硼矿的形成关系;通过对所获得重力资料的处理,清楚地揭示了局部异常的变化和研究区的断裂构造特征;通过对重力异常进行优选延拓和小波变换处理,并结合研究区的其他地质资料的综合分析,对火山沉积型硼矿的分布范围进行了预测,为该地区火山沉积型硼矿的钻探工程部署提供了地球物理依据。