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

简介：基于Aki-Richards公式和贝叶斯原理,本文发展了利用叠前PP波和PS波资料联合反演P波速度比、S波速度比和密度比的方法。该方法假设参数之间满足正态分布,引入参数协方差矩阵来描述反演参数之间的相关性以提高反演过程的稳定性,并同时使反演的参数序列服从Cauchy分布,引入矩阵Q来描述参数序列的稀疏性以提高反演结果的分辨率。采用本文提出的方法对模型数据和实际多波资料进行反演,结果表明：本文方法正确有效;与传统的单一PP波反演相比,PP波和PS波AVO联合反演具有稳定性更好和反演精度更高等优点。

简介：储层物性参数的改变伴随着弹性参数某种程度的改变，而这种改变并非简单的线性关系。加上观测信息的缺失、重叠、噪声破坏及模型理想化等原因，致使这种改变关系的获取具有较大的不确定性。针对传统统计岩石物理物性参数反演方法的不足，以利用弹性参数反演储层物性参数为目的，依据贝叶斯反演框架，我们建立了新的储层物性参数目标反演函数。首先，采用兼具确定性与随机性特点的统计岩石物理模型，考虑到不同弹性参数间的精度存在差异，引入权重系数，建立起储层物性参数与弹性参数间的加权统计关系。其次，基于这种加权统计关系，结合马尔科夫链蒙特卡洛随机模拟技术产生储层物性参数、弹性参数随机联合样本空间作为目标函数求解样本空间。最后，建立解的快速求解准则，求取最大后验概率密度对应的储层物性参数取值作为最终解。实际应用表明，该方法具有较高的反演效率，应用前景较好。

简介：本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法,该方法综合利用岩石物理和地质统计先验信息,在贝叶斯理论框架下,首先通过变差结构分析得到合理的变差函数,进而利用快速傅里叶滑动平均模拟算法（FastFourierTransformMovingAverage,FFT-MA）和逐渐变形算法（GradualDeformationMethod,GDM）得到基于地质统计学的储层物性参数先验信息,然后根据统计岩石物理模型建立弹性参数与储层物性参数之间的关系,构建似然函数,最终利用Metropolis算法实现后验概率密度的抽样,得到物性参数反演结果。并将此方法处理了中国陆上探区的一块实际资料,本方法的反演结果具有较高的分辨率,与测井数据吻合度较高;由于可以直接反演储层物性参数,避免了误差的累积,大大减少了不确定性的传递,且计算效率较高。

简介：本文提出了一种新的反演方法:通过采用纵、横波走时数据对(从相同的震源产生的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的空间分布具有较小的分散性,同时,它的分布特征能较好的与地震波速度结构相吻合。

简介：本文将相干成像用于WVSP数据的一次反射波和下行地表多次波的成像中,通过对叠加次数分布的分析,发现一次反射波相干成像的成像范围小但在井位置处的叠加次数高;下行地表多次波相干成像的成像范围大而在井位置处的叠加次数少。基于此,本文尝试利用一种WVSP的联合相干成像方法来解决这一问题。该方法利用WVSP的直达波与下行地表多次波和直达波与一次反射波来进行联合相干成像。这种联合相干成像方法在有效拓宽成像面积的同时增加了覆盖次数(特别是在井旁),使干涉记录中的虚假同相轴得到了有效地压制,提高了成像剖面的信噪比,使成像结果更加可靠。正演模拟得到的WVSP合成记录和野外实际采集到的WVSP数据处理结果均表明这种方法的有效性。更多还原

简介：常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。

简介：由中国科学院地质与地球物理研究所与中国石油勘探开发研究院联合主办的2009北京岩石物理研讨会才2009年彤月4—6日在北京文津国际酒店召开。采自英国帝国理工、波斯坦地学中心、卡廷大学、

简介：多分量勘探相比常规纵波勘探而言有许多的优势，PP／PS联合AVO分析和反演是强而有效的储层识别方法。在本文中，我们推导了更精确的转换波AVO反射系数公式；并在入射角小于30度时进行了进一步简化，其对于转换波AVO分析和反演而言更加简单有效。基于该近似，我们进行了PP／PS联合AVO反演。实际资料的实例表明，反演得到的纵横波速度比可很好地识别岩性及油气。反演得到的流体因子和泊松比等其他属性在储层处也显示出明显异常，效果显著。

简介：3D地震波走时计算是偏移、反偏移、层析等诸多地震勘探技术中的重要中间步骤。快速推进法计算3D地震波走时具有高效率、稳定性及适应能力强的特点，但快速推进法在震源附近区域的计算精度不高，降低了整个走时算法的计算精度。本文提出了一种联合3D走时计算方法来解决这一问题。该方法在震源附近小范围内使用计算精度较高的波前构建法计算走时，在剩余区域使用快速推进法计算走时，由于模型中绝大多数网格节点走时是通过快速推进法计算的，故新方法保留了快速推进法高效的特点，同时由于震源附近网格节点走时精度的提高，整个新算法的计算精度相对于快速推进法而言有了较大的改善。文中通过数值分析对上述结论进行了验证并使用三维岩丘模型验证了新方法的稳定性和适应能力。

简介：本文提出了一种基于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.

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

简介：Q因子估算是近地表吸收补偿提高地震记录分辨率最为重要的基础工作之一。我们采用了一种新的井地联合地震数据采集方式减小检波器耦合对Q因子估算的影响。在以激发井为中心的圆周上,按照井下检波器的设计深度,布置多口深度不同的接收井,将检波器直接安置在每口井的井底,以消除常规井间观测方式造成的检波器与井壁耦合对Q因子估算的影响。在此基础上,我们提出了一种不受激发影响的Q因子层析反演方法,利用模型数据就该方法的稳定性和可靠性进行了测试分析。使用两个不同的地震频带对大港油田实际近地表观测数据进行了吸收反演,反演结果表明,近地表的吸收系数远大于地下地层的吸收系数,近地表吸收补偿对于提高地震资料分辨率具有重要的现实意义。另外,两个频带反演得到了不同的Q因子,这在一定程度上支持了Q因子对频率依赖性的认识。

简介：JointPP–PSinversionoffersbetteraccuracyandresolutionthanconventionalP-waveinversion.P-andS-waveelasticmodulideterminedthroughdatainversionsarekeyparametersforreservoirevaluationandfluidcharacterization.Inthispaper,startingwiththeexactZoeppritzequationthatrelatesP-andS-wavemoduli,acoefficientthatdescribesthereflectionsofP-andconvertedwavesisestablished.Thismethodeffectivelyavoidserrorintroducedbyapproximationsorindirectcalculations,thusimprovingtheaccuracyoftheinversionresults.Consideringthattheinversionproblemisill-posedandthattheforwardoperatorisnonlinear,priorconstraintsonthemodelparametersandmodifiedlow-frequencyconstraintsarealsointroducedtotheobjectivefunctiontomaketheproblemmoretractable.Thismodifiedobjectivefunctionissolvedovermanyiterationstocontinuouslyoptimizethebackgroundvaluesofthevelocityratio,whichincreasesthestabilityoftheinversionprocess.TestsofvariousmodelsshowthatthemethodeffectivelyimprovestheaccuracyandstabilityofextractingPandS-wavemodulifromunderdetermineddata.Thismethodcanbeappliedtoprovideinferencesforreservoirexplorationandfluidextraction.

简介：探索地震前兆是一件具有挑战性的研究。本文采用高通滤波的方法，处理分析了来自世界不同位置4个超导重力台站5次大地震前的秒级固体潮重力数据，发现在这些地震发生之前数天内，能够接收到稳定的重力高频扰动信号，其中2008年5月12日的汶川地震前出现一个主频为0.1-0.3Hz的重力扰动，而其他4次地震前出现两个主频为0.12-0.17Hz和0.06-0.085Hz的重力扰动。另一个发现是当地震发生在大陆与海洋板块断裂带时这种重力扰动信号具有远场效应，即重力扰动异常出现在距离震中较远的超导重力仪上，而较近的台站异常较小或没有异常。以上研究结果表明，这种重力扰动信号与地震发生具有较密切的相关性，可作为一种有潜力的临震预报信号；重力扰动信号的远场效应也许能够揭示某些地震发生前地球内部板块之间的相互作用机制。但由于地球上重力固体潮台站分布的非常不均匀，上述结果还需要今后研究更多的样本数据加以证实。

简介：本文提出了基于测井、VSP和地震数据拟合的子波估计方法,从输入输出都包含随机噪声的统计模型出发,采用相关性拟合技术来提取子波。拟合度和误差分析为整个过程提供了定量的质量控制手段,可以评估数据拟合和子波估计的可靠性。实际数据试算表明,该方法在含有噪声的实际数据中稳定而有效,在地震频带内的子波估计和数据拟合是可靠的。该方法无需对子波相位和振幅谱进行任何假设,其主要优点在于确定相位的能力。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。

简介：渗透率评价是储层评价中的一个复杂问题，传统测井方法难以对储层渗透率参数进行直接和有效的评估。储层渗透率对钻井过程中的泥浆滤液侵入有着较大的影响，因此本文提出一种利用泥浆侵入效应来评价储层渗透率的数学模型和方法。首先构造含泥饼的泥浆侵入数值模型，然后以达西渗流理论为基础导出泥浆侵人深度和储层渗透率的近似数学关系，以此可以利用侵入深度测量值来估算储层渗透率。对孔隙度、渗透率和含水饱和度三个主要储层参数的敏感性分析，发现该方法适用于受到泥浆侵入影响的低孔隙度、低渗透率的油层或油水同层。采用现场测井和取芯数据进行数值模拟计算，结果表明估算出来的渗透率曲线和预设的渗透率曲线吻合较好，证明该方法具有一定的可行性。

简介：由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。