张海峰, 黄旭日, 张俊, 任梦宇. 基于导航金字塔的强地震约束建模方法——以胜利油田沙二段河流相为例[J]. 石油地球物理勘探, 2024, 59(6): 1197-1205.
Haifeng ZHANG, Xuri HUANG, Jun ZHANG, Mengyu REN. Strong seismic constraint modeling based on steerable pyramid: A case study of fluvial facies in Sha-2 member in Shengli Oilfield[J]. Oil Geophysical Prospecting, 2024, 59(6): 1197-1205.
河流相储层横向变化快、非均质性强,传统建模方法已不能满足目前储层表征的需求。为此,首先利用导航金字塔技术,将地震数据进行分解与重构,重构后的数据更加凸显地质规律,有利于沉积相的刻画;然后采用贝叶斯—序贯高斯地震约束建模方法进行储层建模,建立测井数据与地震属性的约束关系,提高模型垂向分辨率。胜利油田沙二段河流相储层建模实践表明,经导航金字塔技术处理后的地震数据横向表征能力增强,砂体边界识别效果较好;与传统序贯高斯方法相比,贝叶斯—序贯高斯方法模拟模型垂向分辨率更高,砂体预测厚度与实际厚度吻合率达86%,更有利于砂体识别。该方法可为剩余油开采提供指导。
赵亮, 孙小东, 李振春, 秦宁, 王九拴, 杨静. 利用双路卷积神经网络的速度自动拾取方法[J]. 石油地球物理勘探, 2024, 59(6): 1206-1216.
Liang ZHAO, Xiaodong SUN, Zhenchun LI, Ning QIN, Jiushuan WANG, Jing YANG. Automatic velocity picking using dual-path convolutional neural network[J]. Oil Geophysical Prospecting, 2024, 59(6): 1206-1216.
人工拾取速度谱是地震资料常规处理中速度分析的主要手段,此方法既耗时、耗力,又限制了大规模三维地震资料处理的效率和准确性。为此,提出了一种利用双路卷积神经网络的速度谱自动拾取方法。首先,采用卷积神经网络结合注意力机制作为主网络,从速度谱数据中提取能量团的特征并实现速度的自动拾取;其次,训练主网络在输出时间—速度序列之前,通过特征融合和特征转换将速度与另一个卷积神经网络(辅网络)输入的未校正CMP道集的隐藏表征进行信息融合,重构成校正后的CMP道集;最后,通过辅网络模拟CMP道集动校正的过程,利用动校正优化速度拾取的精度。模型和实际资料测试结果表明,在加入辅助神经网络引入动校正信息后,文中方法比单一的卷积神经网络在速度拾取方面具有更高的精度。
武庭润, 高建虎, 常德宽, 王海龙, 陶辉飞, 李沐阳. 基于Transformer的地震数据断层识别[J]. 石油地球物理勘探, 2024, 59(6): 1217-1224.
Tingrun WU, Jianhu GAO, Dekuan CHANG, Hailong WANG, Huifei TAO, Muyang LI. Faults identification of seismic data based on Transformer[J]. Oil Geophysical Prospecting, 2024, 59(6): 1217-1224.
利用地震资料识别断层在油气勘探中有着重要的作用。目前,机器学习和深度学习技术提高了断层识别的精度和效率,但断裂预测结果仍难以满足生产需求。为此,提出基于Transformer的地震断层识别方法,即3D SwinTrans-U-Net。该网络由Swin Transformer模块、卷积模块组成。其中,Swin Transformer模块可以利用Transformer的注意力机制提取全局信息,并将计算全局注意力转变为计算窗口的注意力,从而比Transformer减少了计算复杂度;卷积模块具有归纳偏置的特性,避免了Swin Transformer存在弱归纳偏置的缺陷;最后,利用U-Net结构,结合Swin Transformer层与卷积层,融合深层与浅层的信息并提取相关特征,充分学习全局性和局部依赖性信息,在保证断层识别精度的基础上提高了计算效率,实现端到端的地震断层学习。模型数据和实际数据测试均表明,3D SwinTrans-U-Net网络能进一步提升断层识别精度。
李克文, 李文韬, 窦一民, 朱信源, 阳致煊. 基于AHRFaultSegNet深度学习网络的地震数据断层自动识别[J]. 石油地球物理勘探, 2024, 59(6): 1225-1234.
Kewen LI, Wentao LI, Yimin DOU, Xinyuan ZHU, Zhixuan YANG. Automatic fault recognition in seismic data based on AHRFaultSegNet deep learning network[J]. Oil Geophysical Prospecting, 2024, 59(6): 1225-1234.
断层识别是地震数据解释的重要环节之一。深度学习技术的发展有效提高了断层自动识别的效率和准确性。然而,目前在断层的自动识别任务中,如何准确捕捉断层细微结构并有效抵抗噪声干扰仍然是一个具有挑战性的问题。为此,在HRNet网络的基础上,构建了一种基于解耦自注意力机制的高分辨率断层识别网络模型AHRFaultSegNet。对于自注意力机制解耦,结合空间注意力和通道注意力,代替HRNet中并行传播的卷积层,在减少传统自注意力机制计算量的同时,模型可以在全局范围内计算输入特征的相关性,更准确地建模非局部特征;对解耦自注意力使用残差连接来保留原始特征,在加速模型训练的同时,使模型能够更好地保持细节信息。实验结果表明,所提出的网络模型在Dice、Fmeasure、IoU、Precision、Recall等性能评价指标上均优于其他常见的断层自动识别网络模型。通过对合成地震数据与实际地震数据等进行测试,证明了该方法对断层细微结构具有良好的识别效果并且具有良好的抗噪能力。
首皓, 曾庆才, 胡莲莲, 丁玲, 王彦春, 孙鲁平. 应用曲率谱和Siamese网络的叠前深度偏移速度建模[J]. 石油地球物理勘探, 2024, 59(6): 1235-1243.
Hao SHOU, Qingcai ZENG, Lianlian HU, Ling DING, Yanchun WANG, Luping SUN. Prestack depth migration velocity modeling with curvature spectrum and Siamese network[J]. Oil Geophysical Prospecting, 2024, 59(6): 1235-1243.
速度建模是叠前深度偏移的重要环节,通常需要在层位约束下对观测点的地层速度进行横向外推,然而在速度建模初期缺乏地震解释层位等格架信息。为此,文中提出了一种基于曲率谱横向相似性和改进循环结构Siamese网络的速度模型建立方法。Siamese网络是目前常用的基于深度学习的目标识别和追踪网络,可以快速进行目标图像的相似度对比,而且不需要人工制作标签。曲率谱可以看成反应地层特征和速度信息的二维图像,将速度建模作为横向特征相似性类比问题,通过类比曲率谱可以自动得到地层的格架和速度更新信息。首先,将叠前深度偏移后的道集转换为曲率谱;其次,确定待搜索曲率谱图像及其对应的目标追踪对象,并求取当前追踪对象与目标追踪对象的相似系数;然后,基于相似系数更新参考曲率谱图像和当前追踪对象;最后,在遍历完全部追踪对象时,基于各个追踪对象的层速度及深度建立速度模型。理论模型和实际数据试验结果表明,该方法能在没有解释资料的条件下快速生成符合地质构造和地层特征的速度模型。
邬世英. 时频域保幅面波压制方法[J]. 石油地球物理勘探, 2024, 59(6): 1244-1251.
Shiying WU. Amplitude-preserving surface wave suppression method in time-frequency domain[J]. Oil Geophysical Prospecting, 2024, 59(6): 1244-1251.
陆上勘探中面波是常见干扰波之一,面波压制的保幅效果对后续资料处理具有较大影响。传统的面波压制方法很难在有效保护反射波的情况下对面波进行精确压制。由于面波具有低速、低频、高振幅的特征,文中提出了一种基于Gabor变换的时频域保幅面波压制方法。首先,利用Gabor变换将地震道划分成若干个地震道片段,在地震道片段内根据面波和有效波振幅能量差异进行面波识别,从而实现时间域内面波和有效波精细分离;然后,在含有面波的地震道片段上采取频率域带通滤波,根据面波和有效波频率的差异进行面波压制;最后,利用Gabor反变换得到面波压制后的地震道。模型测试和实际数据处理表明,该方法在有效压制面波的同时,最大程度地保护了有效波的信息,提高了面波处理的保幅性。
吴迪, 文武, 门哲, 马一凡. 基于扩散模型的地震数据随机噪声压制方法[J]. 石油地球物理勘探, 2024, 59(6): 1252-1259.
Di WU, Wu WEN, Zhe MEN, Yifan MA. A random noise suppression method for seismic data based on diffusion model[J]. Oil Geophysical Prospecting, 2024, 59(6): 1252-1259.
地震数据中的随机噪声毫无规律,常规去噪方法难以达到理想的效果,影响后续的地震数据解释和分析。为此,提出一种基于扩散模型的地震信号去噪方法。该方法的前向扩散过程是通过对地震数据进行一定程度的加噪,将地震数据变成存在大量各向同性的高斯噪声的含噪地震数据,再利用训练后的扩散模型对含噪数据进行重建,提高地震数据的信噪比。预测网络部分是基于改进的U-Net网络,该网络中引入了注意力模块和ResNet模块,以提高网络对重要区域的关注度,避免深度网络中的梯度消失问题。理论数据和实际数据的应用结果均验证了文中方法的有效性。该方法去噪效果远超FX滤波、SVD等传统去噪方法,同时也比经典的深度学习网络CNN、GAN更加优秀,能够完整地保留有效信号,极大提升地震数据的质量。
周继继, 许辉群, 李彪, 杨梦琼, 黄燕洁. 利用动态时间规整的多层位追踪方法[J]. 石油地球物理勘探, 2024, 59(6): 1260-1268.
Jiji ZHOU, Huiqun XU, Biao LI, Mengqiong YANG, Yanjie HUANG. A multihorizon tracking method using dynamic time warping[J]. Oil Geophysical Prospecting, 2024, 59(6): 1260-1268.
层位追踪是地震资料处理、解释中非常重要的步骤,现有的层位自动追踪技术在遇到断层时追踪效果不佳。为此,根据相邻地震道之间波形相似性,提出一种利用动态时间规整(DTW)的纯数据驱动的多层位追踪方法。首先,提取地震道特征值,将地震道按波谷、波峰或过零点等特征划分,使所追层位点严格遵循以上特征,保证追踪结果的高精度;其次,利用DTW方法计算两道相邻地震道划分好的特征值序列以得到相似路径;最后,从相似路径中提取所有彼此相似的特征值点对,并根据参考层划分特征值点对得到各套层位。模型数据试算和实际资料测试结果表明,所提方法能快速追踪到目标区域的各套层位,并有效克服断层对层位追踪的阻碍,对不同地质条件的地震资料具有一定的适用性,有一定的应用价值。
张元鹏, 张强, 王玲谦, 周辉, 董春晖, 张明珠. 应用自适应方向TV正则化的多道反褶积方法[J]. 石油地球物理勘探, 2024, 59(6): 1269-1279.
Yuanpeng ZHANG, Qiang ZHANG, Lingqian WANG, Hui ZHOU, Chunhui DONG, Mingzhu ZHANG. Multichannel deconvolution based on adaptively directional TV regularization[J]. Oil Geophysical Prospecting, 2024, 59(6): 1269-1279.
作为提高地震数据纵向分辨率的有效手段,反褶积技术一直是地震勘探领域的研究热点。传统反褶积方法采取逐道稀疏反演的策略重构地下反射系数序列,忽视了相邻地震道之间的空间关系,反褶积结果空间连续性不佳。应用全变分(Total Variation,TV)约束的多道地震反褶积方法虽然能够改善由逐道反演引起的空间不连续,但由于未考虑地质构造信息,因而仅适用于“块状”构造。为了克服这一问题,文中提出一种自适应方向TV(Adaptively Directional TV,ADTV)约束的多道地震反褶积算法。该方法基于局部地质构造方向构建ADTV正则化项,不像传统TV约束仅考虑水平和垂直方向的信息。因此,ADTV可以更好地保护弱小反射信号,更清晰地刻画地质体的空间展布,在地质构造复杂的情况下也能保证反褶积结果的空间连续性。模型测试和实际资料处理结果均表明,该方法可以提供分辨率高、空间连续性强的反褶积结果。
李默言, 黄旭日, 史云燕, 陈永灿, 黎孟承. 基于P、z分量Radon变换和非稳态自适应相减的OBN横波噪声压制[J]. 石油地球物理勘探, 2024, 59(6): 1280-1285.
Moyan LI, Xuri HUANG, Yunyan SHI, Yongcan CHEN, Mengcheng LI. OBN shear wave noise suppression based on P and z component Radon transform and unsteady adaptive subtraction[J]. Oil Geophysical Prospecting, 2024, 59(6): 1280-1285.
在OBN地震采集中,检波器与海底耦合问题经常导致陆检z分量接收到横波噪声。横波噪声对OBN数据的处理、成像、解释、油藏监测均带来严重干扰,而且横波噪声的频带、曲率与有效信号非常相似,难以用常规方法有效压制。为此,提出了一种基于抛物线Radon变换和非稳态自适应相减的OBN横波噪声压制方法。首先,将OBN数据的P分量和z分量数据通过抛物线Radon变换从时空域变换到τ-q域,利用横波噪声与反射信号传播速度的差别,在τ-q域尽可能分离横波噪声和一次波;其次,通过非稳态自适应相减算法的高度灵活性,从包含横波噪声的z分量数据中去除P分量中的反射信号成分,提取横波噪声;最后,将提取的横波噪声分量变换回时空域,并从原始z分量数据中去除横波噪声,实现消除噪声的同时,尽可能保护地震数据中的反射信号。实际资料的处理结果证明,该方法能有效压制OBN数据的横波噪声,提高地震资料的信噪比。
马继涛, 赵凯歌, 廖震. 应用主频带约束高精度Radon变换的多次波压制方法[J]. 石油地球物理勘探, 2024, 59(6): 1286-1296.
Jitao MA, Kaige ZHAO, Zhen LIAO. Multiple attenuation by high precision Radon transform with dominant frequency band constraint[J]. Oil Geophysical Prospecting, 2024, 59(6): 1286-1296.
Radon变换是地震资料处理中常用的一种算法,通过特定路径求和实现地震数据插值、多次波压制和波场分离。然而,由于采集数据和算法限制,变换域的分辨率较低,提高Radon域分辨率一直是研究的热点。最常用的提高分辨率方法是多次迭代重加权算法,通过迭代更新权重,将加权值聚焦到地震数据曲率上,但很难将加权值聚焦到真实曲率位置。文中提出了一种优化提高变换域分辨率的新方法,即在地震数据的主频带内求取一个加权矩阵,并将加权值聚焦至真实地震曲率位置。首先,计算地震数据的主频,并以主频为中心频率,向低频和高频方向各取一个范围作为约束频带;对主频带内的地震数据采用低频约束策略,从低频到高频迭代计算加权矩阵;从主频带的最后一个加权矩阵得到最终的加权矩阵,并应用于所有其他频率的计算。主频带数据的信噪比高、振幅强,因此,该方法更加稳定,得到的加权矩阵可以显著增强变换域的分辨率。此外,与其他迭代方法相比,它避免了每个频率加权矩阵的迭代过程,计算效率更高。合成数据和实际数据的测试结果证明了该方法的在多次波压制方面的有效性和优势。
李斐, 许磊明, 李美, 陈娟, 黄黎刚, 吴朦朦. 提高黄土山地三维地震资料信噪比和分辨率的保真处理方法[J]. 石油地球物理勘探, 2024, 59(6): 1297-1304.
Fei LI, Leiming XU, Mei LI, Juan CHEN, Ligang HUANG, Mengmeng WU. Fidelity processing methods for improving signal-to-noise ratio and resolution of 3D seismic data in loess mountainous landforms[J]. Oil Geophysical Prospecting, 2024, 59(6): 1297-1304.
黄土山地地表条件复杂、干扰强,表层介质非均质性强,同时受表层巨厚黄土吸收衰减影响,出现地震资料信噪比低、主频低、频带窄等问题。经叠前保真噪声衰减处理和叠前时间偏移后的叠加数据仍存在随机噪声及散射干扰,制约了薄储层的预测精度。业界通常采用F-XY域四维去噪+叠后零相位反褶积方法提高信噪比、分辨率,笔者则采用F-X域Cadzow滤波实现三维叠后数据的散射干扰及随机噪声同步压制,在此基础上利用连续小波变换方法提高分辨率。各处理环节中,应用噪声剖面、F-K谱、振幅谱、时间切片、分频扫描、井震标定等多种质控手段确保参数的合理选择。该套处理方法具有良好的保真性,可为后续反演、属性分析等环节提供良好的数据基础。
江雨濛, 周昕, 郝越翔, 孙耀光, 陈雪, 李宜真. 基于重加权TL1范数的非稳态谱反演算法[J]. 石油地球物理勘探, 2024, 59(6): 1305-1314.
Yumeng JIANG, Xin ZHOU, Yuexiang HAO, Yaoguang SUN, Xue CHEN, Yizhen LI. Nonstationary spectral inversion algorithm based on reweighted TL1 norm[J]. Oil Geophysical Prospecting, 2024, 59(6): 1305-1314.
地震数据反射系数反演是联结地下储层和地震数据的桥梁,一直是研究的热点。目前反射系数反演大多基于L1范数约束的稀疏脉冲反褶积。近年来,奇偶分解算法的出现使得子波间调谐效应减弱,这使得基于L1范数约束的谱反演得到进一步应用。稀疏约束的能力关系到待求解反射系数的准确性,本文针对常用的L1范数和Lp范数稀疏能力约束度不足的问题,引入TL1范数(Transformed L1 Norm)稀疏约束,有利于获得更准确的反演反射系数;同时考虑反射系数较大的位置拟合能力需要增强,提出重加权TL1范数(Reweighted Transformed L1 Norm,RTL1)进一步提高稀疏约束能力。参数测试结果证明了重加权类范数重建能力强于未重加权类范数,RTL1在稀疏重建上的有效性。模型和实际数据处理结果证明RTL1相比于常用的稀疏约束项更有利于提高谱反演的反射系数精度。
马雄, 申天赐, 田永净, 韩刚, 桂志先, 彭晓波. 频率域稳定化的黏声波逆时偏移方法及其应用[J]. 石油地球物理勘探, 2024, 59(6): 1315-1329.
Xiong MA, Tianci SHEN, Yongjing TIAN, Gang HAN, Zhixian GUI, Xiaobo PENG. Stabilized viscoacoustic reverse time migration in the frequency domain and its application[J]. Oil Geophysical Prospecting, 2024, 59(6): 1315-1329.
地震波的吸收衰减和频散均与频率有关,因此,在频率域研究黏声波逆时偏移具有完备的理论基础和明确的物理含义。与时间域黏声波逆时偏移一样,频率域方法在补偿地层吸收衰减时也将不可避免的引起高频噪声放大,造成吸收衰减补偿的不稳定问题,影响偏移算法的稳定性和成像精度。为此,利用Kolsky-Futterman模型在频率域建立黏声波波动方程,并据此实现了一种稳定化的黏声波逆时偏移方法。其实现过程如下:首先,利用Kolsky-Futterman模型推导出频率域的黏声波波动方程,在该方程中振幅衰减和相速度频散是解耦的;然后,基于该方程的解耦特性,利用黏声波波场和仅含频散波场的稳定化比值,构建一种稳定化的吸收衰减补偿算子;最后,利用稳定化的吸收衰减补偿算子对震源正传波场和检波点反传波场分别进行吸收补偿处理,并应用互相关成像条件对地下结构进行黏声波偏移成像。模型实验和实际资料应用表明:该方法在补偿地层吸收的同时,较好地抑制了高频噪声的放大,保证了偏移算法的稳定性和成像精度,改善了成像剖面的质量。
杨春生, 焦艳丽, 蔡东梅, 杨莹彬, 边婧. 井—震结合曲流型储层构型研究及应用——以大庆长垣西部斜坡区为例[J]. 石油地球物理勘探, 2024, 59(6): 1330-1339.
Chunsheng YANG, Yanli JIAO, Dongmei CAI, Yingbin YANG, Jing BIAN. Meandering reservoir architecture based on well-seismic analysis and its application: Taking the western slope area of Daqing Oilfield as an example[J]. Oil Geophysical Prospecting, 2024, 59(6): 1330-1339.
大庆油田西部斜坡区萨尔图油层萨Ⅰ组储层以曲流型河道砂体沉积为主,由于点坝砂体边界及其内部废弃河道识别精度低,制约着开发水平井的部署。为此,在高密度地震资料基础上,采用基于动态Q补偿的偏移处理技术,大幅提高了河道砂体的成像精度,明确了地震剖面上曲流型河道砂体点坝及其内部废弃河道的波组特征;平面上为进一步突出曲流型河道砂体内部非均质特征,利用楔状地震正演模型,建立了不同频率地震数据体可识别的极限砂体厚度,根据研究区砂体厚度分布值域,优选基于S变换的3个单频体采用RGB三色融合技术,实现了不同厚度级别砂体精准预测;同时,建立了点坝砂体边界及其内部废弃河道的岩心相-测井相-地震相的响应模式,明确了点坝及废弃河道的含油性。研究结果表明:基于动态Q补偿的偏移处理技术,提高了河道砂体在地震剖面上的识别精度;基于S变换的分频体RGB属性融合技术,能够在平面上精准预测废弃河道、点坝边界及其内部非均质性;建立的岩心相—测井相—地震相的响应模式,能够准确预测不同类型砂体的含油性。应用研究成果指导部署井位304口,其中188口直井钻井成功率98.0%,116口水平井平均实钻水平段长度387.8 m,平均含油砂岩钻遇率93.4%。
李创, 程冰洁, 徐天吉, 王啟, 华豪. 基于OVT域HTI介质衰减各向异性的储层裂缝预测方法及应用[J]. 石油地球物理勘探, 2024, 59(6): 1340-1352.
Chuang LI, Bingjie CHENG, Tianji XU, Qi WANG, Hao HUA. Reservoir fracture prediction based on OVT domain HTI medium attenuation anisotropy and its application[J]. Oil Geophysical Prospecting, 2024, 59(6): 1340-1352.
品质因子Q是表征地震波衰减的一种重要参数,广泛应用于储层裂缝发育特征及含油气性识别。然而,在各向同性、各向异性、黏弹性等不同介质模型的约束下获取的Q值具有不同的储层表征能力。为了精确描述裂缝型、裂缝—孔隙型、孔隙—裂缝型等储层的各向异性及裂缝发育特征,基于HTI介质模型和地震波衰减各向异性理论,提出了一种基于OVT(Offset Vector Tile)域HTI介质衰减各向异性的储层裂缝预测方法。首先,利用广义S变换(GST)对OVT矢量片叠前道集进行高精度时频分析;然后,通过T-K(Teager-Kaiser)能量算子计算地震波的单频T-K能量;最后,计算抽样波长距离内的相对衰减量,获得衰减拟Q值进行椭圆拟合,并以椭圆扁率表征衰减各向异性强度和预测裂缝发育密度。该方法不仅充分发挥了OVT叠前道集较叠后数据包含了更丰富的储层流体和裂缝信息的优势,而且将衰减属性与储层裂缝表征联系在一起,提供了一种裂缝预测新思路。该方法在川西坳陷XC气田的应用中取得了地震资料预测结果与岩心分析、测井资料解释成果吻合度较高的裂缝识别效果。
丁富峰, 吴少波, 姚志纯, 马玉龙, 林鹏, 汪泳吉. 基于精确Zoeppritz方程的储层流体识别[J]. 石油地球物理勘探, 2024, 59(6): 1353-1361.
Fufeng DING, Shaobo WU, Zhichun YAO, Yulong MA, Peng LIN, Yongji WANG. Reservoir fluid identification based on exact Zoeppritz equation[J]. Oil Geophysical Prospecting, 2024, 59(6): 1353-1361.
流体识别对于储层的预测至关重要。在各向同性介质中,Russell流体因子的获取主要基于Zoeppritz方程近似公式的反演或者弹性参数间接计算,然而,基于近似公式的反演容易受到小入射角和低对比度等假设的干扰,影响反演精度。为此,提出了一种基于精确Zoeppritz方程的储层流体识别方法。首先,将精确Zoeppritz方程重新改写为包含Russell流体因子、泊松比和密度的形式;其次,在贝叶斯反演框架下构建了流体因子、泊松比和密度同时反演的目标函数;最后,通过迭代重加权最小二乘法(IRLS)对反演目标函数进行非线性求解。合成数据和实际工区数据测试结果均表明,所提方法能够有效识别流体。该方法可为储层流体预测、储层流体标定工作提供新的研究思路。
刘旭锋, 李凤磊, 刘华, 徐杰. 牛庄地区沙三中亚段重力流沉积储层预测技术与应用[J]. 石油地球物理勘探, 2024, 59(6): 1362-1377.
Xufeng LIU, Fenglei LI, Hua LIU, Jie XU. Prediction technologies and its application of reservoirs for gravity flow sedimentary in the middle Third Member of Shahejie Formation in Niuzhuang Area[J]. Oil Geophysical Prospecting, 2024, 59(6): 1362-1377.
断陷湖盆重力流沉积砂体相变快、规模相对较小、形成机制及分布规律复杂,决定了该类储层地震响应特征复杂。通过有针对性地分析现有三维地震资料品质,优选适合进行目标处理的地震数据体开展研究。利用沙四段油页岩地震响应特征、多井正负极性子波分别标定等手段,判断资料极性,为目标处理打下基础。针对牛庄油田沙三中亚段地球物理特征、浊积砂岩体发育特征,设计不同砂体组合的理论正演模型和实际砂体正演模型,确定合理的参数完成三参数小波高分辨率目标处理方法,效果显著。获得以下认识:①处理后地震数据体高频成分较老资料明显拓宽,且实现了复波分离,波组特征清楚,合成记录标定显示,处理后资料品质明显优于原纯波资料;②处理后资料波组连续性提高,能量关系清晰,有利于浊积砂体等的识别;③处理后资料利于标定,通过多井标定,建立精细构造模型与井约束模型完成地震反演,预测结果与实钻吻合度较高,在详细地层划分对比、地质研究基础上,采用地层切片技术,完成对沙三中亚段第5~第8砂层组主力含油砂岩体纵向发育、平面分布研究,确定砂体边界,为断陷湖盆重力流砂体储层精细刻画提供了新的思路。
刘智颖, 徐松, 罗兵, 张绪平, 许巍. 基于微扰法的裂缝发育层段的感应测井电磁场解析解[J]. 石油地球物理勘探, 2024, 59(6): 1378-1394.
Zhiying LIU, Song XU, Bing LUO, Xuping ZHANG, Wei XU. Analytical solution of induction logging electromagnetic field in fracture development formation based on perturbation method[J]. Oil Geophysical Prospecting, 2024, 59(6): 1378-1394.
根据微扰法的思想,将裂缝的影响看作原状地层电磁场的微扰,在径向分层介质的亥姆霍兹方程解的基础上,研究了柱外无穷小天线振子电磁场的解析解,从而提出了一种适用于任意裂缝产状,以及任意井径、任意泥浆电阻率的感应测井电磁场的解析解法。据此,以中国常用的HDIL 1515 EA/MA阵列感应测井仪为例,计算了4种裂缝性地层模型中接收线圈的感应电动势,并绘制了视电导率曲线(理论图版)。研究表明,在裂缝发育层段,阵列感应视电导率曲线并不会出现“大幅度尖峰”的现象。同时,在油基泥浆条件下,微裂缝密集发育的地层段中,7个子阵列的视电导率曲线呈现出幅度差异,且其特征明显不同于泥浆侵入层段,可以作为微裂缝发育段的识别标志。
钟志豪, 李予国. 基于GLC多项式谱元法的大地电磁场二维正演模拟[J]. 石油地球物理勘探, 2024, 59(6): 1395-1409.
Zhihao ZHONG, Yuguo LI. Two-dimensional forward modeling of magnetotelluric field based on GLC polynomial spectral element method[J]. Oil Geophysical Prospecting, 2024, 59(6): 1395-1409.
为了提高大地电磁场数值模拟的精度和效率,提出基于Gauss-Lobatto-Chebyshev(GLC)基函数谱元法的大地电磁场正演模拟方法。该方法首先从麦克斯韦方程组出发,推导了二维大地电磁场边值问题;然后基于Galerkin加权余量法将微分形式的边值问题转换成积分弱形式;最后采用GLC正交多项式插值基函数对全局问题进行离散化,利用Pardiso求解器求解大型稀疏线性方程组得到大地电磁场,实现了二维大地电磁场数值模拟。为了提高数值模拟效率,文中算法采用变密度规则网格剖分技术,即在电性复杂区域使用细网格,在电性均匀区域使用粗网格,采用OpenMP编程模式实现了多个频点的并行计算,以此达到缩短计算时间的目的。通过一维层状介质模型数值模拟结果验证了文中算法的正确性和精度,相比于Gauss-Lobatto-Legendre(GLL)多项式谱元法的数值模拟结果,GLC多项式谱元法的模拟结果精度更高。针对国际标准模型COMMEMI 2D-1和带地形模型,进行了基于GLC多项式谱元法、有限差分法和基于三角网格有限元法的正演模拟,三者数值结果对比表明,GLC多项式谱元法计算精度更高、网格依赖性更低。
秦策, 付文亮, 王绪本, 米晓利, 杨云见, 赵宁. 应用非一致性网格和高阶有限元的井地电磁三维各向异性正演算法[J]. 石油地球物理勘探, 2024, 59(6): 1410-1419.
Ce QIN, Wenliang FU, Xuben WANG, Xiaoli MI, Yunjian YANG, Ning ZHAO. Three-dimensional anisotropic forward modeling algorithm for borehole-to-surface electromagnetic method using non-conforming meshes and high-order finite element[J]. Oil Geophysical Prospecting, 2024, 59(6): 1410-1419.
地下介质通常表现为电各向异性,增加了井地电磁法响应的复杂性。因此,针对各向异性介质模型进行正演模拟并总结规律,对于正确解释电磁勘探的观测数据至关重要。目前,对复杂三维地质模型进行精确建模是正演中的难点之一。本文使用一种基于八叉树网格的高精度建模方法,通过将六面体单元划分为八个较小的单元来局部细化网格,能够使用较少数量单元的同时保证建模的精度,从而降低计算成本。此外,为提高数值解的精度,在计算时使用高阶形函数。结合八叉树网格和高阶形函数实现了井地电磁法三维各向异性正演方法。通过对比各向异性层状模型的数值解与解析解,验证了算法的正确性。使用八叉树网格实现复杂各向异性油气模型精确建模,将异常体设置不同的水平和垂直电阻率,计算结果表明垂直电阻率对响应的影响大于水平电阻率。研究成果对实际的野外施工具有一定的指导作用,为进一步进行各向异性介质的三维反演打下基础。
李曼, 张志勇, 刘佳成, 杨军, 鄢毛毛, 周小宝. 先验岩石物性约束的地球物理正则化反演[J]. 石油地球物理勘探, 2024, 59(6): 1420-1432.
Man LI, Zhiyong ZHANG, Jiacheng LIU, Jun YANG, Maomao YAN, Xiaobao ZHOU. Geophysical regularization inversion constrained by priori petrophysical properties[J]. Oil Geophysical Prospecting, 2024, 59(6): 1420-1432.
随着勘探程度的深入,可用于正则化反演的先验信息越来越丰富和精确。在分析参数变换函数、模糊C均值聚类模型约束、交替方向乘子优化等引入岩石先验物性约束方法的基础上,提出了一种基于反演单元物性互异原则的模型约束方法。基于单元互异原则的模型约束从地质单元数量有限、每个离散反演单元只能属于某一种岩性的本质出发,通过限定离散反演单元的物性取值实现物性引入。将上述先验物性引入与经典正则化反演相结合,构建了统一的目标函数并通过高斯—牛顿法进行优化求解;对四种方案实现岩石物性引入的作用机理和数学本质进行了分析,提出了同时利用多种方案提高反演效果的组合策略;对比了L1和L2范数最小结构模型约束与先验岩石物性引入策略结合的反演效果。模型和实测数据反演结果均表明,充分利用岩石物性测量的先验信息可有效提高反演效果。
钟连诚, 李伟勤, 刘红岐, 马琰祁. 基于深度学习的可控源电磁法GEMTIP激电参数预测[J]. 石油地球物理勘探, 2024, 59(6): 1433-1446.
Liancheng ZHONG, Weiqin LI, Hongqi LIU, Yanqi MA. Deep learning based controlled source electromagnetic method for GEMTIP induced polarization parameter extraction[J]. Oil Geophysical Prospecting, 2024, 59(6): 1433-1446.
广义有效介质极化理论(Generalized Effective Medium Theory for Induced Polarization,GEMTIP)提供了岩石物理参数与复电阻率(Complex Resistivity,CR)的频散关系,据此可基于观测到的激发极化(Induced Polarization,IP)数据反演岩矿石的激电参数。然而,传统的反演方法在非线性优化问题上存在局部最小值、计算量大和对初始模型依赖度高等问题,且含噪数据反演结果不稳定。此外,当前的激电参数反演研究主要集中在微观岩石孔隙表征和电化学机制领域,基于宏观地球物理观测数据直接进行反演和解释的相关研究不足。为此,提出了一种基于U-Net深度学习网络的方法,利用该方法可基于GEMTIP三维地电模型的地面IP差分数据直接提取激电参数。该方法将多个频率下的差分观测磁场振幅和相位作为网络输入,将异常区域的零频电阻率、体积分数、充电率、时间常数及弛豫常数作为输出标签。通过合成GEMTIP三维地电模型的可控源电磁样本数据训练深度神经网络,得到能够准确预测地下异常区域激电参数分布的网络模型。对包含GEMTIP激电参数的综合模型测试了该网络模型,并将测试结果与传统的正则化共轭梯度(Regularized Conjugate Gradients,RCG)反演结果进行比较,表明U-Net网络反演在耗时、求解精度和抗噪声能力方面均更具优势,能够从地面观测到的IP数据中直接、准确地预测GEMTIP激电参数。最后,利用深度学习方法对亚利桑那州南部North Sliver Bell地区的辉铜矿实际勘测数据进行训练,成功预测了该地区地下辉铜矿富集层分布,并与传统反演方法获得的地质解释成果进行对比,进一步证明了本文方法在实际应用中的可靠性和有效性。该方法可用于矿物组成和储层孔隙空间分布的预测,有望在宏观地球物理反演解释中得到广泛应用。
杨俊, 刘雪军, 王财富, 胡祖志, 何展翔, 吴鹏. 时频电磁法在深层页岩气有利区预测中的应用[J]. 石油地球物理勘探, 2024, 59(6): 1447-1457.
Jun YANG, Xuejun LIU, Caifu WANG, Zuzhi HU, Zhanxiang HE, Peng WU. Application of time-frequency electromagnetic method in predicting favorable areas for deep shale gas[J]. Oil Geophysical Prospecting, 2024, 59(6): 1447-1457.
时频电磁法基于时—频域、电—磁场一体化探测理念,实现了多信息联合,弥补了传统电磁方法抗干扰能力不足的缺点,集成了可控源电磁法的优点:探测深度大,精度高,抗干扰能力强。在四川泸州地区开展时频电磁勘探,结合岩石物理和测井数据,对宝藏—云锦向斜的志留系五峰组—奥陶系龙马溪组页岩层系进行页岩气有利区预测。对本区外露头信息和电阻率测井数据进行统计,结果表明页岩含气性与电阻率具有相关性,高电阻率页岩段普遍表现出较高的含气性,本区具有开展电法勘探的物性基础。研究区不同地层之间存在明显的电性差异,纵向可划分为高、低相间的6套电性层,其中页岩气主力产层五峰—龙马溪组与上、下围岩电阻率差异明显。从时频电磁勘探原理出发,利用一维有限元和有限差分方法进行数值模拟,研究方法的可行性,并根据正演结果设计最佳采集参数。在围绕目的层进行针对性数据采集的基础上,采用自主研发的重、磁、电处理解释系统GeoGME进行数据处理解释,通过弱异常信息提取和井震双控约束反演等关键技术,探明了龙一1亚段各小层的电阻率和极化率分布特征。综合考虑目标页岩层段电阻率异常和极化率异常,预测出页岩气勘探有利区4个,为下一步该区域页岩气钻井部署提供了技术支撑,表明时频电磁法是进行深层页岩气有利区预测的有效勘探手段。
宋莎, 董怀民, 曹晓萌. 南极南设得兰陆缘天然气水合物储层地球物理研究进展[J]. 石油地球物理勘探, 2024, 59(6): 1458-1468.
Sha SONG, Huaimin DONG, Xiaomeng CAO. Geophysical research progress on gas hydrate reservoirs in the South Shetland margin (Antarctica)[J]. Oil Geophysical Prospecting, 2024, 59(6): 1458-1468.
南极南设得兰陆缘天然气水合物分布广泛,因其在资源、地质灾害及全球气候变化等方面具有重要作用,成为国内外研究的热点。为此,从天然气水合物储层的地震识别特征、水合物相关地质构造、水合物储层定量评价以及水合物与气候变化关系等方面对南设得兰陆缘水合物储层的地球物理研究工作进行了梳理和总结,获得了几点思考与认识,以期为今后天然气水合物勘探开发及气候变化研究提供一定的参考借鉴。最后对该区域水合物储层地球物理研究方向进行了展望,认为综合利用海底地震仪数据的纵横波信息开展水合物储层三维表征研究、深入研究天然气水合物对全球变暖的响应机制及其反馈作用等是未来重要的发展方向。