简介：传统Robinson褶积模型主要受缚于三种不合理的假设,即白噪反射系数、最小相位地震子波与稳态假设,而现代反射系数反演方法（如稀疏约束反褶积等）均在前两个假设上寻求突破的同时却忽视了一个重要事实：实际地震信号具有典型的非稳态特征,这直接冲击着反射系数反演中地震子波不随时间变化的这一基础性假设。本文首先通过实际反射系数测试证实,非稳态效应造成重要信息无法得到有效展现,且对深层影响尤为严重。为校正非稳态影响,本文从描述非稳态方面具有普适性的非稳态褶积模型出发,借助对数域的衰减曲线指导检测非稳态影响并以此实现对非稳态均衡与校正。与常规不同,本文利用对数域Gabor反褶积仅移除非稳态影响,而将分离震源子波和反射系数的任务交给具有更符合实际条件的稀疏约束反褶积处理,因此结合两种反褶积技术即可有效解决非稳态特征影响,又能避免反射系数和地震子波理想化假设的不利影响。海上地震资料的应用实际表明,校正非稳态影响有助于恢复更丰富的反射系数信息,使得与地质沉积和构造相关的细节特征得到更加清晰的展现。

简介：非局部均值滤波是一种基于图像信息冗余的去噪方法，其认为图像自身的有效结构具有一定的重复性，而随机噪声则不具备这一特点，通过利用图像本身的自相似性来达到压制随机噪声的目的，是一种全局的去噪方法。本文把这一思想引入地震数据随机噪声压制中，针对传统非局部均值滤波计算量过大的问题，文章采用分块非局部均值的方式来减少计算量；针对滤波参数选取会影响非局部均值滤波效果的问题，提出一种简单的自适应滤波参数地震数据分块非局部均值算法。模型和实际数据处理结果表明：相对于传统的去噪算法（如f-x反褶积），该方法在压制随机噪声的同时对有效信号保护地更好，具有更高的保真度，更有利于后续的处理和解释工作。

简介：常规长排列非双曲动校正公式是在VTI介质中得到的,它不能满足任意空间取向TI（ATI）条件下的扩展.本文以VTI介质中非双曲动校正公式为基础,基于我们推导得出的ATI介质中精确四次时差系数解析解和NMO速度解析解,给出ATI介质中长排列优化的非双曲动校正公式.通过与各向异性射线追踪方法计算所得出的＂精确走时＂结果对比,研究表明优化后的非双曲动校正公式能精确地描述任意强弱、ATI介质中随测线方位变化的走时曲线,可以用来替代耗时、多偏移距、多方位的射线追踪方法正演拟合ATI介质中长偏移距反射走时,为利用非双曲时距的各向异性参数反演提供理论基础性认识。

简介：通过最新的全球地磁模型——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）处衰减较快,在地表处则趋于稳定。

简介：通过引入岩石骨架参数表达式到Gassmann近似流体方程中，实现了基于该方程的地震孔隙度反演。但是岩石骨架参数表达式有很多种，使得地震孔隙度反演公式也各不相同，实际应用不方便，反演结果也难以评价。作者在研究现有常用岩石骨架参数表达式如Esheby－Walsh，Pride，Geertsma，Nur，Keys－Xu以及Krief等的基础上，提出了含有调节参数的岩石骨架模型统一表示式，发展了基于Gassmann方程的地震孔隙度反演方法。该方法应用范围宽，参数调节方便、灵活。为验证所提公式和方法的实用性，作者们结合ZJ地区钻井岩心的岩石物理样品测试数据和测井数据，进行了地震孔隙度反演。由于该地区油气储层的存在与中等程度孔隙度的分布有关，因此地震孔隙度反演对预测和识别油气层、干层和含水层有重要作用。反演结果与工区内孔隙度测井瞌线吻合度高，说明本文所提公式与方法是有效的，可靠的。

简介：在储层发育带，通常地震波波形会变复杂，所以关联维会变小。但由于地震波是有一定带宽的信号，波形受这个带宽内所有频率成分的影响，如果直接进行关联维的计算，其结果往往不能很好地反映储层。鉴于这种情况，提出了经验模态分解（EMD）与关联维相结合的方法，利用EMD对地震波进行分解，对分解后的每一个平稳的固有模函数（IntrinsicModeFunction，IMF）分量进行关联维计算，将计算结果与有效的IMF分量进行储层预测。利用该方法对xx区三叠系中油组进行了试算，取得了较好的效果。

简介：本文是中国陆域航磁地质构造研究最新进展的综述性文章，对基于航磁数据的中国陆域地质构造编图最新成果和取得的新认识进行了论述。该项研究利用我国50多年的海量实测航空磁测数据，采用最新的数字化编图技术及先进实用的位场数据处理方法编制了1：5000000万和1：1000000万中国陆域航磁系列图件，以板块构造和大陆动力学理论为指导，以岩石物性研究与航磁异常特征分析为基础，定性、定量与综合解释相结合，编制了中国陆域磁性基底深度图及主要地质构造单元图件，并依据航磁数据对新划分的地质构造进行了论述。

简介：Gabor变换和S变换是常用的时频分析工具。根据测不准原理,它们的时频分解结果无法在时间域和频率域同时具有很高的分辨率。为了提高非平稳信号时频分解结果的分辨率,本文提出瞬时频率分布函数（IFDF）并利用它表达非平稳信号。当非平稳信号时频成分的分布满足测不准原理对信号可分辨的要求时,瞬时频率分布函数的支集和短时Fourier变换的小波脊支集是同一个集合。利用IFDF的该特征,本文提出一种迭代算法（Sparse-STFT）实现了信号的稀疏时频分解。该算法在每次迭代过程中利用残留信号的短时Fourier变换结果的脊支集更新信号的时频成分,每次迭代得到的时频成分的叠加结果即为最终的稀疏时频分解结果。文中的数值实验证明了Sparse-STFT可以有效地提高非平稳信号时频分解结果的分辨率。最后,本文将该方法应用于地震数据面波的压制中,取得了理想的处理结果。

简介：本文通过地震正演模型和面波振幅数据对比分析了天然地震、爆炸,与核试验数据,以探索改进辨别天然地震与爆炸的方法。所提方法是基于波数域和频率域的双积分变换解。用于研究的地震事件都发生在北朝鲜,其中包括2001年6月26日（39.212°N,125.383°E）和2001年10月30日（38.748°N,125.267°E）所记录的爆炸数据,2006年10月9日（41.275°N,129.095°E）所记录的一次核试验数据、以及2002年4月14日（39.207°N,125.686°E）和2002年6月7日（38.703°N,125.638°E）所记录的两次天然地震。基于这些数据通过地震波分析与理论模拟表明爆炸引起的地震与天然地震有着不同的波型特征。爆炸引起的信号特征是P波的能量较S波强,在爆炸记录上0.05-0.5Hz频率之间Rg波清晰呈现,而在然地震记录上没有。这是由于爆炸地震记录上P波是优势波,与SH成分发生了耦合。

简介：数据挖掘是从大量的、不完全的、有噪声的、模糊的数据中提取隐含的、但又是潜在有用的信息和知识的过程，由于数据挖掘具有出色的非线性建模能力和自组织学习能力，因此可以在复杂储层的测井解释中发挥作用。本文用数据挖掘方法识别复杂储层的岩性。将岩性识别作为一种分类任务建立数据挖掘流程，包括特征提取、特征选择和建立模型等步骤。本文用独立成分分析法从测井曲线中提取信息；然后使用分支定界算法寻找最佳的特征子集，并消除冗佘信息；最后采用C5．0决策树算法建立分类模型的测井曲线。模型和实际测井数据吻合较好，表明在复杂油藏的研究中数据挖掘方法是有效的。

简介：常规地震观测系统设计方法基于地下水平层状介质的假设,通常不能适应复杂构造情况。我们从控制照明的思想出发,提出了一种面向目标成像的地震观测系统设计方法,该方法需要一个由初步地震解释得到的速度模型。利用单程傅立叶有限差分波场传播算子将目标层的平面源延拓到地表,通过分析从目标层延拓到地表的波场能量的分布,可以确定目标层成像所需要的炮点或者检波点的位置。利用SEG-EAGE盐丘模型数值试算结果表明,该方法用于设计面向目标成像的特定地震采集系统。

简介：本文通过人工变换T2分布和建立管-球模型模拟法研究含水合物地层渗透率与水合物含量之间的关系。首先,在渗透率的模拟试验中,我们改变了束缚水与可动水的比例、总孔隙度以及与之关联的T2分布。试验结果表明,相对渗透率与水合物含量之间的关系受到这些因素的制约。随后,我们用管-球模型表示水合物生长的孔隙空间,并把水合物的生长过程看成是向孔隙空间随机扔小球的过程。在此过程中,采用两种方法计算渗透率,一是Schlumberger＇sT2公式（即SDR模型）,二是Darcy定律与Poiseuille流动方程相结合的方法。前人的实验研究表明,在一定的水合物含量范围内,渗透率基本保持不变。以此为参考,我们将计算结果与之进行比较。我们发现,采用SDR模型时,渗透率的数值模拟曲线与Masuda模型N=15时的结果相近。而采用Darcy定律时,渗透率模拟值较高,但与实验结果的趋势相一致,都会出现渗透率的平直阶段。尤其,当水合物晶体在孔隙体内优先生成时,优先的概率越高,渗透率的平直范围越大。

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

简介：在地震频段内,中观尺度局域流是引起孔隙岩层速度频散与衰减的主要原因,研究中观尺度孔隙中局域流的过渡频率及其对应的衰减（逆品质因子）峰值频率之间的关系,有助于深入分析孔隙介质内部构造的细节。基于Biot孔隙弹性理论得出时间域内局域流流体通量,引入傅立叶变换得到频域上的流体通量。为避免不同单元体之间局域流的相互影响,建立双相孔隙介质单重单元体模型,并且在流体通量基础上选择含水层的过渡频率来等价局域流的过渡频率。此外,针对不同孔隙介质参数,分析了局域流整体过渡频率及相应的衰减峰值频率的变化特征。为了了解边界条件对计算结果造成的影响,建立了多重单元体模型,引出局部局域流过渡频率及相应的衰减峰值频率。数值模拟结果表明,在低饱和度状态下,两种频率的变化趋势相同,但变化幅度存在较大差异;在高饱和度状态下,两种频率的变化趋势与变化幅度都很接近。

简介：多波勘探是提高致密、低渗透等复杂油气藏勘探开发精度的有效手段,但转换波属于典型的低信噪比、低分辨率资料,如何最大程度提高转换波纵向分辨率已成为多波处理的一个难点。谱反演技术能有效解决常规反褶积提频技术受频带限制影响,提高分辨率能力有限的问题,能最大程度提高数据分辨率,便于识别薄层,但该技术的难点是如何通过稳定的反演算法得到高精度反射系数,如何把反射系数恢复为宽频的地震数据？本文在前人研究的基础上利用全变差作为先验信息有助于求解欠定问题的优势,提出了一种基于全变差约束的最小二乘反演算法,提高反演的精度和稳定性,并利用高斯拟合振幅谱模拟得到宽频子波数据,通过在蓬莱地区的高分辨率数据恢复处理,得到分辨率更高的转换波。理论试验和实际资料应用证明了该方法能较大幅度地提高转换波资料分辨率,为后续速度反演、储层反射信息提取提供更准确的数据。

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

简介：致密储层与围岩波阻抗差异小，常规储层反演方法难以有效预测出有利区。中国西部致密储层非均质性强，常规储层反演预测难度大，本文提出了基于Xu-White模型的叠前储层预测方法。首先修正Xu-White模型的试验参数，直至该模型计算出来的纵、横波速度及密度与实际测井得到的数值基本吻合，然后用修正的模型计算不同储层物性、岩性、含油性下的纵、横波速度、密度及泊松比，进而建立不同的砂泥岩互层地质模型。其次，采用Zoeppritz方程近似式进行叠前正演，分析不同储层物性、岩性、含油性时的叠前道集响应特征，在对中国西部Z区实际叠前道集优化处理的基础上，选择其中1种对岩性、含油性、物性最敏感的弹性参数，对其有利储层分布区进行预测，其预测结果具有较高的精度。

简介：目前对动电测井的研究主要集中于数值模拟和岩石物理实验。为研发出能下井测量的动电测井仪器,进行了一系列实验研究工作。文中介绍了一种融合了不同结构和材料要求的新型声电复合式动电测井探测器的结构及原理,对声、电信号接收器在仪器机械设计上无法布置在同一位置的难题也提出了解决方法。在水池中对其辐射声场特性进行了测试,计算了辐射声场声压及发射电压响应,分析了基于相控线阵的阵列式声波发射探头的指向性及实际应用效果。研究表明,在源距1500mm处产生的声压可达47.2kPa,且会随着激励声源主频增大而减小。随着相邻声波发射子阵激励信号延迟时间的增加,辐射声束主瓣会发生明显偏转,主瓣方向的能量也逐渐增强,可有效增强声电转换效率。本文为探测器研制及开展井下动电测井研究奠定了重要基础。

简介：常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。

简介：通过地震数据获取裂缝储藏中流体的性质并对流体类型进行识别,是地震勘探岩性反演的重要问题之一。由于地震波的速度、储层的密度等弹性参数对某些流体不具有很强的敏感性,使只依赖振幅信息进行流体识别的传统AVO方法面临困境。作为传统叠前振幅反演的一个拓展,频变AVO（FDAVO）技术进一步考虑了振幅对频率的依赖关系,将这种依赖关系与地下裂缝结构、流体填充对应起来,能带来更丰富的流体信息。利用该技术,本文提出了一种基于地震数据参数化Chapman模型的贝叶斯反演新方法（BIDCMP）,它包含两步算法,即,FDAVO反演储层的非弹性属性和贝叶斯框架下的流体识别。首先,通过匹配观测数据和模型数据,构造差函数反演裂缝储层非弹性参数。随后,在贝叶斯框架下,使用马尔科夫随机场（MRF）作为先验模型,联合多参数场识别流体。本方法在计算过程中,除综合考虑了弹性参数场、测井资料等常规信息外,还特别地加人了第一步中反演得的非弹性参数的约束,从而充分利用了流体粘性差异,最后在最大后验概率（MAP）准则下输出最佳岩性一流体识别结果。分别对合成地震记录和模拟岩性—流体剖面验证本文方法的有效性,结果证明本文方法获得的流体识别结果准确可信。