简介:本文给出一种既能有效衰减地震噪音又可保护地层及构造的不连续性的新方法。构造约束保边平滑技术需要已知反射局部方位和边界信息,通常这些信息由全频率地震资料估算获得,但在资料信噪比很低的情况下,噪音往往会降低估算的可靠度。对于信噪比极低的地震资料,其主频成分相对非主频成分信噪比高,所以由主频资料获取的方位和边界信息比由其它频率成分获取的更可靠。方位和边界信息通常用倾角和相干值差异来描述。由于不同频率所引起的倾角和相干值差异的变化均比地震记录的变化缓慢,所以由主频资料获取的倾角及边界信息能够近似代表所有频率成分的倾角及边界信息。Ricker子波广泛用于地震勘探,Marr小波与Ricker子波在时间和频率域均具有相同的形态,所以选用Marrl小波变换将地震数据按照倍频程分为几个分频体。扫描主频分频体,用不等权二次曲面拟合并求解极大值来获取视倾角,通过比较9个滑动窗口的相干值来确定反射边界。将这些信息用构造约束保边平滑技术可选择性地(selectively)对主频、低频、高频分频体做平滑处理,最后将平滑后的各频段地震记录合成为滤波去噪后的地震记录。理论模型和实际资料处理效果表明该方法能有效压制噪音,保护边界,保护同相轴的连续性,且灵活地保留地震记录中的有用信息。
简介:Q因子估算是近地表吸收补偿提高地震记录分辨率最为重要的基础工作之一。我们采用了一种新的井地联合地震数据采集方式减小检波器耦合对Q因子估算的影响。在以激发井为中心的圆周上,按照井下检波器的设计深度,布置多口深度不同的接收井,将检波器直接安置在每口井的井底,以消除常规井间观测方式造成的检波器与井壁耦合对Q因子估算的影响。在此基础上,我们提出了一种不受激发影响的Q因子层析反演方法,利用模型数据就该方法的稳定性和可靠性进行了测试分析。使用两个不同的地震频带对大港油田实际近地表观测数据进行了吸收反演,反演结果表明,近地表的吸收系数远大于地下地层的吸收系数,近地表吸收补偿对于提高地震资料分辨率具有重要的现实意义。另外,两个频带反演得到了不同的Q因子,这在一定程度上支持了Q因子对频率依赖性的认识。
简介:相干属性是检测地下介质不连续性特征的主要参数之一,但倾斜地层会对相干属性产生干扰,故在计算相干属性时需要消除地层倾角的影响。传统的相干属性计算一般使用时域倾角扫描法来消除地层倾角的影响,但这种方法的精度受到采样间隔等因素的限制,并不能很好地消除地层倾角的影响。根据振幅谱具有不受地层倾角影响的特点,本文提出将预先定义的子数据体内的地震数据变换为振幅谱,再利用振幅谱构造协方差矩阵,协方差矩阵最大特征根与矩阵迹的比值作为子数据体中心点的相干值。利用振幅谱计算相干属性的方法能够较好地消除地层倾角的影响,而且可以使用不同频段的振幅谱计算相干体,具备多尺度相干体算法的特性。模型数据和实际地震数据的应用结果表明,振幅谱相干属性可以更好地消除地层倾角对相干属性的影响。
简介:我们改进了共偏移距-共反射面(COCRS)法,可用以衰减地滚波,即由于低速、低频/高振幅瑞雷波通常产生的相干噪声。COCRS算子是基于双曲线,因此它可以拟合双曲走时的同相轴,如叠前数据中的反射同相轴。相反,地滚波在共中点(CMP)和共炮点道集中是线性的并可以可以利用COCRS算子鉴别与压制。因此,我们在共偏移距剖面之前共炮道集中搜索反射倾斜和曲率。因为这对反射振幅的危害最小化是最理想的,我们只对在地滚波区多次覆盖的数据进行叠加。在CO剖面前搜索CS道聚集是对常规COCRS叠加的另一个改进。我们使用合成和真实数据集测试了所提出的方法,数据采自伊朗西部地区。我们将本方法压制地滚波的结果与f-k滤波和f-k滤波后常规COCRS叠加压制地滚波的结果进行了比较。结果表明,该方法对真伪滚压制效果优于F-K滤波与传统CRS叠加。然而,计算时间高于其他常规的方法,如f滤波。
简介:地震数据重构是地震数据处理的重要步骤之一,重构算法的精度、效率与抗噪性是地震数据重构技术的核心研究内容。研究针对傅里叶域凸集投影(POCS)算法,在定义的最优阈值评价标准基础上,提出了反比例阈值模型,该模型具有在大系数区间比指数模型更快下降速率、而在小系数区间比指数模型更慢下降速率,从而在保证弱反射信号重构精度的同时有效提高POCS地震数据重构算法计算效率。为提高反比例阈值对不同地震数据特点的适应性,在地震数据谱能量分布差异性特征分析基础上,研究提出了在反比例阈值模型分母上增加适应地震数据谱能量特征的因变参数,通过调节该因变参数获得适应不同地震数据特点的最佳阈值曲线,进一步提高算法的计算精度与计算效率。为了实现重构过程中随机噪音的自适应衰减,提高重构后地震数据信噪比,研究提出了数据驱动的加权回加系数计算策略,利用每次迭代对应数据驱动阈值占阈值区间的百分比获得加权回加系数。研究将新方法应用于模拟三维数据和实际三维地震数据,分析结果表明反比例阈值相对传统阈值在提高数据重构计算效率和精度方面具有明显的优越性,新提出的加权回加系数计算策略能有效提高重构数据的信噪比。
简介:目前在研究地球物理参数的空间变异特性时常常是基于几何各向异性假设,然而实际中几何各向异性假设常常不成立。本文提出了一种求解地球物理参数各向异性的新方法,首先对各方向上的变异值作归一化处理,然后用拟合的方法得到各方向单位步长上的变异性强度,进而实现对全局各向异性的求解。以上海地区合成孔径雷达干涉测量(InSAR)中差分大气延迟样本数据为例,验证了新方法求解各向异性的优越性:新方法的偏差只有6.4%,而基于几何各向异性假设方法的偏差达到了21.2%,而且新方法的运算效率有非常显著的提高。进一步,通过克里金插值交叉验证实验,证明了基于新算法得到的结构函数的插值效果最好,从而从另外一个角度验证了新算法求解的各向异性更准确,更好地表征了区域化变量的空间结构特征。因此,与几何各向异性假设的方法相比,新方法不仅能更准确的描述地球物理参数的各向异性,而且有更高的运算效率,为更准确地估计出所需的地球物理参数奠定了基础。
简介:目前对动电测井的研究主要集中于数值模拟和岩石物理实验。为研发出能下井测量的动电测井仪器,进行了一系列实验研究工作。文中介绍了一种融合了不同结构和材料要求的新型声电复合式动电测井探测器的结构及原理,对声、电信号接收器在仪器机械设计上无法布置在同一位置的难题也提出了解决方法。在水池中对其辐射声场特性进行了测试,计算了辐射声场声压及发射电压响应,分析了基于相控线阵的阵列式声波发射探头的指向性及实际应用效果。研究表明,在源距1500mm处产生的声压可达47.2kPa,且会随着激励声源主频增大而减小。随着相邻声波发射子阵激励信号延迟时间的增加,辐射声束主瓣会发生明显偏转,主瓣方向的能量也逐渐增强,可有效增强声电转换效率。本文为探测器研制及开展井下动电测井研究奠定了重要基础。
简介:致密砂岩储层的孔隙结构对其渗透性和电性影响显著,是此类复杂储层岩石物理研究的关键。针对仅从连通喉道半径评价渗透率的多解性以及储层孔隙结构与电性关系研究欠缺等不足,综合影响物性的主要因素,提出了一种同时考虑孔隙度、最大连通喉道半径及分选性三种因素的新型孔隙结构参数δ的计算公式。利用岩心及实测数据对比分析表明,δ值能够较连通喉道半径等传统方法更精确地刻画致密砂岩储层渗透性,同时它与储层电性具有密切关系,可用于估算地层因素F和胶结指数m。据此提出将孔隙结构对电阻率的影响进行归一化校正以及基于核磁共振测井预测储层完全含水电阻率R0的评价方法,从而突出储层流体性质变化引起的电性变化,并提供了一种新的致密砂岩储层流体识别思路,研究结果得到了实验资料和实际测井试油资料的验证。
简介:有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程,高阶显式有限差分方法通常用来提高求解精度,已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数,发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明:忽略这些小的差分系数能够近似维持有限差分的精度,但是显著减小了计算量。然后,引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性,采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后,给出了针对均匀和非均匀介质模型的弹性波模拟例子,表明了本文方法的优点。
简介:Inthispaper,wetheoreticallyandnumericallystudyacombinedmonopole–dipolemeasurementmodetoshowitscapabilitytoovercometheissuesencounteredinconventionalsingle-wellimaging,i.e.,thelowsignal-to-noiseratioofthereflectionsandazimuthambiguity.First,theazimuthambiguity,whichexistsextensivelyinconventionalsingle-wellimaging,issolvedwithanimprovedimagingprocedureusingcombinedmonopole–dipoleloggingdatainadditiontoconventionalloggingdata.Furthermore,wedemonstratethatthedirectwavespropagatingalongtheboreholeswithstrongenergy,canbeeffectivelyeliminatedwiththeproposedcombinedmonopole–dipolemeasurementmode.Thereflectionsarethereforepredominantinthecombinedmonopole–dipoledataevenbeforethesignalsarefiltered;thus,thereflections’arrivaltimesineachreceiverareidentified,whichmayhelpminimizethedifficultiesinfilteringconventionalloggingdata.Theoptimizedprocessingflowofthecombinedmeasurementmodeloggingimageisgiveninthispaper.Theproposedcombinedmonopole–dipolemeasurementmodemayimprovetheaccuracyofsingle-wellimaging.
简介:常规欧拉反褶积法中构造指数的选取以及分散解存在较多的问题,提出了基于联立垂向一阶导数与解析信号的欧拉齐次方程的RDAS-Euler反演方法。该方法可以更为精确的估计场源的范围及埋深,且不需考虑构造指数N的影响,避免了因构造指数不当而引起的反演误差。通过对单一地质体及组合地质体模型的实验证明本文方法能有效地完成目标体的反演工作,反演结果与理论值之间的误差小于10%,且相对于常规欧拉反褶积法更加稳定准确,能够更好的得到地质体边界及深度信息。将RDAS-Euler法应用于黑龙江省虎林盆地实测布格重力异常数据,获得了丰富的断裂信息,说明RDASEuler法增强了对断裂平面位置的识别能力。
简介:河流相储层是一种重要的储层类型,寻找河道砂体是油气勘探的一项重要工作。但是河道变迁比较频繁,叠置比较严重,加之地震资料本身信噪比和分辨率的限制,河流相储层的精细描述至今仍存在一定的困难。本文以胜利油田老河口三维工区为例,首先研究了本区河流相储层的基本特征,发现该类储层剖面上呈强相位、切片上连续性好和频谱上呈低频特征;通过颜色处理和地震精细解释,发现本区河流相储层存在类串珠状特征,为了验证这一观测结果,我们用正演模拟分析了其形成机理,为在类似工区识别此类古河道找到了新的途径;应用多属性融合和RGB显示技术,河道下切特征更加明显,很好地彰显了河道结构特征,提高了河道识别能力;最后,我们研究并开发了多子波检测技术,检测出更多河流相弱反射信息。
简介:探索地震前兆是一件具有挑战性的研究。本文采用高通滤波的方法,处理分析了来自世界不同位置4个超导重力台站5次大地震前的秒级固体潮重力数据,发现在这些地震发生之前数天内,能够接收到稳定的重力高频扰动信号,其中2008年5月12日的汶川地震前出现一个主频为0.1-0.3Hz的重力扰动,而其他4次地震前出现两个主频为0.12-0.17Hz和0.06-0.085Hz的重力扰动。另一个发现是当地震发生在大陆与海洋板块断裂带时这种重力扰动信号具有远场效应,即重力扰动异常出现在距离震中较远的超导重力仪上,而较近的台站异常较小或没有异常。以上研究结果表明,这种重力扰动信号与地震发生具有较密切的相关性,可作为一种有潜力的临震预报信号;重力扰动信号的远场效应也许能够揭示某些地震发生前地球内部板块之间的相互作用机制。但由于地球上重力固体潮台站分布的非常不均匀,上述结果还需要今后研究更多的样本数据加以证实。
简介:边界识别是重磁数据解释中的常用方法之一,依据其结果可划分出地质体的水平范围。边界识别结果受地质体埋深及导数计算误差的影响所识别边界与真实边界之间存在一定的差距,且边界识别法无法直观地给出地质体的深度信息。为了获得异常体的水平位置和深度信息,本文提出空间归一化边界识别方法,其对不同深度的边界识别函数进行归一化计算,空间归一化边界识别法的最大值对应于异常体的水平位置和深度。常规边界识别结果的误差随理深的减小而减小,而空间归一化边界识别法是通过最大值来判断地质体的位置,最大值是在地质体处获得,因此归一化边界识别方法所获得的结果是准确的。通过理论模型试验证明归一化边界识别方法能有效地完成异常体的水平位置和深度的计算,所获得的水平位置和深度信息与理论值相一致,为下一步的勘探计划提供了更加可靠的依据。将其应用于实际航磁数据的解释,获得了断裂的具体分布形式。更多还原
简介:由于地震数据中包含的噪声在不同频率或者频带数据中的分布强度存在差异,使得全频带数据上进行的噪声衰减处理改变了地震反射波信号的动力学特征,干扰后期的地震资料解释、储层预测、油气检测等问题,提出边界和振幅特性保持自适应噪声衰减方法。首先应用小波包变换对全频带地震数据进行多频段划分,然后对分频段数据进行非线性各向异性倾角导向边界保持自适应滤波处理。在该方法中,由结构张量计算的扩散张量实现自适应地确定平滑滤波方向,加入的不连续结构置信度量和不连续性算子自适应地控制不连续结构特征的保持程度,引入的去相关滤波迭代停止准则自适应地确定滤波迭代次数。这些参数的引入具有减少处理人员的干预和人为的主观性,且执行简单的特点。对合成地震记录和实际地震记录处理结果表明,提议的方法能够自适应地衰减地震数据中噪声,同时既能保持地震反射波中有效的不连续性信息,也能有效地保持有效信号的频率分布规律。能够为后期的地震资料解释和分析提供高品质的基础数据。
简介:本文基于弹性波动方程,从其弱形式出发,利用Galerkin变分原理,通过对方程进行空间和时间上的离散,在空间域中引入预条件共轭梯度的逐元算法,在时间域中引入时间积分的交错网格预处理/多次校正算法,发展了弹性波模拟的Chebyshev谱元算法。针对均匀固体介质和具有倾斜分层的分区均匀固体介质模型,通过与有限差分算法结果相比较验证其精度的可信性,同时利用该算法模拟了弹性波在具有水平分层的任意起伏自由表面模型中的传播,并分析了其传播特点。研究表明,我们提出的交错网格预处理/多次校正算法的Chebyshev谱元算法,保留了有限元法的优势,并且采用了具有最优张量乘积技术的元到元的算法,能够处理带有起伏自由表面的复杂介质模型,它具有比有限元法收敛快,计算效率较高等优点,特别适合于复杂结构和复杂介质中的弹性波传播的数值模拟。