简介:在文献中,Kirchhoff型反偏移场的稳相分析主要是在下列两个条件下进行的:(1)等时面和目标反射面相切;(2)深度偏移像场信号的长度接近于零。对于与目标反射面不相切的等时面和长度远大于零的深度偏移像场子波,已有的结果将不再成立。为了在等时面和目标反射面不相切和深度偏移像场的子波长度远大于零的条件(一般条件)下对Kirchhoff型反偏移场进行稳相分析,我推导了出现在二维稳相分析公式中的诸因子的计算公式,并从中发现:(1)对于不同的等时面,距离差函数的稳相点具有不同的水平坐标;(2)Kirchhoff型真振幅反偏移的输出场由两部分(真振幅反偏移信号与振幅畸变因子)的乘积组成。由此得到下列两个结论:(1)一个给定的反偏移信号由多个深度偏移信号上的采样点组装而成,反偏移信号上的采样点个数等于对于这种组装有贡献的偏移信号的个数。(2)振幅畸变效应是Kirchhoff型反偏移中的固有效应,靠反偏移本身无法消除。如果一定要消除这种振幅畸变效应,必须对反偏移结果进行振幅校正。
简介:本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。更多还原
简介:塔里木盆地碳酸盐岩岩溶型储层广泛分布于塔北、塔中、巴楚等地区,是塔里木盆地碳酸盐岩勘探的重点,也是塔里木油田增储上产的主要层系之一.由于埋藏深度大、反射信号饺酢⑿旁氡冉系?更由于岩溶型储层非均质性严重,使得储层预测的难度增大.通过多年的攻关研究,从碳酸盐岩风化壳储层的地球物理响应特征出发,逐步探索并建立了技术、岩溶型储层测井评价技术、岩溶型储层地震描述技术、应变量分析技术为手段的溶型储层综合评价技术.在塔北、以地质概念模型为基础,以基于层序格架的潜山古地貌恢复塔中等地区应用中取得了良好的地质效果,相继获得了工业油气流.本文介绍了该技术在塔北轮古西地区的具体实施和应用效果.
简介:为了探明介质的非完全弹性对瑞利型槽波传播的影响以及利用瑞利型槽波品质因子QR反演煤层碳化程度、裂隙、煤层厚度等煤层属性,在本文中,我们将复速度引入瑞利型槽波频散方程,计算了煤层瑞利型槽波的品质因子QR曲线;采用控制变量法,分析了瑞利型槽波品质因子QR随煤层厚度、煤层纵横波品质因子以及围岩纵横波品质因子的变化规律。研究表明:瑞利型槽波品质因子QR曲线与群速度曲线变化趋势一致;埃里相(Airy-phase)频率附近瑞利型槽波的品质因子QR最小,且埃里相频率随煤层厚度增加而减小;瑞利型槽波品质因子QR随煤层横波品质因子QS2的增大而增大。利用瑞利型槽波品质因子QR曲线可以实现对瑞利型槽波的吸收衰减进行补偿、对煤层厚度进行预测以及岩性成像。更多还原
简介:有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层(PerfectMatchedLayer,PML)吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的"反射波",以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。
简介:Q因子估算是近地表吸收补偿提高地震记录分辨率最为重要的基础工作之一。我们采用了一种新的井地联合地震数据采集方式减小检波器耦合对Q因子估算的影响。在以激发井为中心的圆周上,按照井下检波器的设计深度,布置多口深度不同的接收井,将检波器直接安置在每口井的井底,以消除常规井间观测方式造成的检波器与井壁耦合对Q因子估算的影响。在此基础上,我们提出了一种不受激发影响的Q因子层析反演方法,利用模型数据就该方法的稳定性和可靠性进行了测试分析。使用两个不同的地震频带对大港油田实际近地表观测数据进行了吸收反演,反演结果表明,近地表的吸收系数远大于地下地层的吸收系数,近地表吸收补偿对于提高地震资料分辨率具有重要的现实意义。另外,两个频带反演得到了不同的Q因子,这在一定程度上支持了Q因子对频率依赖性的认识。
简介:通过最新的全球地磁模型——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)处衰减较快,在地表处则趋于稳定。
简介:常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。
简介:3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点,但快速推进法在震源附近区域的计算精度不高,降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时,在剩余区域使用快速推进法计算走时,由于模型中绝大多数网格节点走时是通过快速推进法计算的,故新方法保留了快速推进法高效的特点,同时由于震源附近网格节点走时精度的提高,整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。
简介:为了有效处理人工源的影响,本文开发了带源的CSAMT二维正反演算法,可用于全区(近区、过渡区和远区)资料的反演。引入正则化因子完成磁法二维反演,并且将模型参数调整为磁化率的对数,保证反演过程中磁化率始终为正值。本文基于交叉梯度原理,将CSAMT和磁法进行联合反演,通过搜索交叉梯度项权重的方法,避免了不同异常源引起的两种异常相互干扰的问题。理论模型算例表明基于交叉梯度的联合反演方法优于单独反演。本文开发的带源CSAMT二维正反演算法,有效处理了人工源的影响,保证了最终联合反演算法的可靠性。
简介:火山沉积型硼矿是一种重要的硼矿床成因类型,主要发育在第三纪湖相蒸发岩沉积层和火山岩互层的二元结构单元中。针对硼矿具有低密度的特点,本文对西藏雄巴地区进行了1:5万高精度重力勘探,获得了研究区布格重力异常和剩余重力异常。根据其异常特征进一步分析了区域构造特征、断裂体系、沉积单元以及火成岩的分布与火山沉积型硼矿的形成关系;通过对所获得重力资料的处理,清楚地揭示了局部异常的变化和研究区的断裂构造特征;通过对重力异常进行优选延拓和小波变换处理,并结合研究区的其他地质资料的综合分析,对火山沉积型硼矿的分布范围进行了预测,为该地区火山沉积型硼矿的钻探工程部署提供了地球物理依据。
简介:Tilt—depth法可以用于快速反演磁源的上顶埋深,但其无法反演磁源底部深度,同时该方法的反演解过于单一。针对于此,本文在有限厚度台阶磁异常导数公式基础上,推导出了可以同时反演地质体上、下界面埋深的改进型Tilt—depth法计算公式,并采用在Tilt梯度图上选取多特征点进行联合反演的计算模式来提高反演解的可靠性。二维及三维组合模型试验均证实了改进型Tilt—depth法可以有效地反演出磁性体的上顶与下底深度,且上顶埋深的反演精度明显高于常规方法的。最后将改进型Tilt—depth法用于松辽盆地长岭断陷航磁数据反演,其上顶埋深的反演结果与钻孔钻遇火山岩的深度基本一致,证实改进型Tilt—depth法相对于常规方法的反演结果更加准确。
简介:本文提出了一种基于Kirchhoff积分偏移和逆时偏移的联合速度分析及成像方法,采用剩余曲率分析及层剥离策略进行偏移速度建模。本文方法改善了Kirchhoff积分偏移在复杂构造时计算精度不高和逆时偏移计算效率慢的缺点,兼有计算效率高、成像精度高的优点,并将其应用在反射波法隧道超前预报中。通过模型试算,发现隧道中使用逆时偏移的成像结果在多方面优于Kirchhoff积分偏移的成像结果;通过对实测数据的处理,验证本方法计算效率较高,建立的速度模型合理,成像剖面清晰,结合地质调绘资料可对隧道开挖前方地质构造做出较准确预报。
简介: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.
简介:本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的P和S波被同一台站记录)来联合反演纵波速度(Vp)和纵、横波速度比(Vp/Vs),然后单独反演横波速度Vs,在反演过程中同时对地震参数进行定位。该方法不需要假设P和S波的射线路径一致,它是沿着P和S波射线路径计算相对慢度扰动值。该方法直接把Vp/Vs作为一个模型参数,由此能获得比采用从独立反演获得的Vp和Vs计算出Vp/Vs的方法更精确的速度比值。该新方法被应用到反演日本东北地区的壳幔速度及波速比结构的研究中,获得了较好的效果。反演结果表明,在日本东北地区,太平洋俯冲板块为一高Vp,高Vs和低Vp/Vs异常区,而在活火山下方的浅部地幔楔以及背弧深部地区为低Vp,低VS和高Vp/VS异常。虽然这些特征在前人的研究中已经报道过,但与前人的研究结果相比,本次研究所获得的Vp/Vs的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。
简介:JointPP–PSinversionoffersbetteraccuracyandresolutionthanconventionalP-waveinversion.P-andS-waveelasticmodulideterminedthroughdatainversionsarekeyparametersforreservoirevaluationandfluidcharacterization.Inthispaper,startingwiththeexactZoeppritzequationthatrelatesP-andS-wavemoduli,acoefficientthatdescribesthereflectionsofP-andconvertedwavesisestablished.Thismethodeffectivelyavoidserrorintroducedbyapproximationsorindirectcalculations,thusimprovingtheaccuracyoftheinversionresults.Consideringthattheinversionproblemisill-posedandthattheforwardoperatorisnonlinear,priorconstraintsonthemodelparametersandmodifiedlow-frequencyconstraintsarealsointroducedtotheobjectivefunctiontomaketheproblemmoretractable.Thismodifiedobjectivefunctionissolvedovermanyiterationstocontinuouslyoptimizethebackgroundvaluesofthevelocityratio,whichincreasesthestabilityoftheinversionprocess.TestsofvariousmodelsshowthatthemethodeffectivelyimprovestheaccuracyandstabilityofextractingPandS-wavemodulifromunderdetermineddata.Thismethodcanbeappliedtoprovideinferencesforreservoirexplorationandfluidextraction.