准东地区云质岩储层叠前反演及脆性预测

脆性是一种评价致密储层可压裂性的重要指标,对油气勘探与开发有着重要的意义。准东地区勘探程度低、探井少,二叠系中下部云质岩致密储层的精细描述和脆性的准确预测均有着很大的难度。为此,根据岩石物理实验结果,利用叠前地震角度域数据体,开展储层叠前反演及脆性预测。通过AVO反演,获取P(截距)+G(梯度)和P×G属性;通过Y(杨氏模量,E)P(泊松比,σ)D(密度,ρ)和L(拉梅常数,λ)M(μ)R(ρ)反演,获取弹性参数的变化率(反射系数);通过测井约束反演,获取弹性参数的实际值,并用多种表征公式计算脆性,分析其物理含义及应用效果。研究结果表明:①研究区横波速度与纵波速度有较好的拟合关系,杨氏模量、泊松比和密度对泥岩、非泥岩区分度明显;②云质岩具有高速特征,研究区储层呈双“甜点”结构;③用三个小角度数据体直接反演弹性参数反射系数,YPD反演方法相较LMR的储层识别效果更好、分辨率更高,密度反射系数数据体可见比较明显的扇体及河道冲积特征;④脆性表征结果与弹性参数特征有一定的区别,脆性表征的云质岩有利区呈条带状分布,而弹性参数预测的云质岩有利区基本呈连续、成片分布;⑤ρE/σ为研究区最佳的脆性表征公式,可为类似油田云质岩储层脆性预测提供借鉴。

卷积神经网络在AVA反演应用中影响因素研究

从基于卷积神经网络开展AVA(振幅随入射角变化)反演问题出发,建立网络训练所需的数据集,分析不同超参数调整对预测结果的影响,确定这些超参数的最佳设置值,并总结超参数调整规律。最后通过对比3种不同卷积神经网络模型的预测结果,寻找出适用于叠前参数反演问题的最佳网络。最终形成一套完备的基于叠前角度域道集数据应用卷积神经网络进行弹性参数预测工作流程。结果表明,基于本文构建的训练数据集,选取的卷积神经网络结构以及设置的超参数组合预测得到的弹性参数结果具有较高的精度。

深反射地震成像揭示的班公湖-怒江缝合带中段Moho断阶及其大地构造意义

揭示班公湖-怒江(班怒)缝合带Moho(莫霍面)结构对于认识中特提斯洋壳俯冲和南羌塘坳陷成因具有重要地球动力学意义。基于横跨班怒缝合带的深反射地震数据(88°30′E),本文采用了中长波长静校正、噪声压制、优化叠加和叠前深度偏移(PSDM)等地震处理技术,获得了深度域地震反射偏移剖面、层速度场和高分辨率Moho结构。由深度域剖面显示,班怒缝合带Moho位于地表以下65~80 km,呈不连续北向抬升趋势,指示在拉萨地块与南羌塘地块之间存在岩石圈上地幔断阶,最大阶步可达15 km。综合分析缝合带两侧的Moho形态认为,这些断阶受南侧拉萨地体的岩石圈上地幔以19.5°北倾俯冲与北侧南羌塘地块的上地壳抬升驱动,可能与深部存在局部熔融相关。班怒缝合带下的Moho结构表明,随着晚侏罗世—早白垩世中特提斯洋闭合,南羌塘地体由边缘海沉积向前陆盆地转换,形成南羌塘坳陷。

一种基于改进门控循环单元的叠前时变子波提取方法

子波的精确提取是地震勘探后续反演与成像的前提,针对传统时变子波提取方法受到的各类假设限制,且需分别提取子波振幅谱与相位谱的问题,本文提出了一种基于改进门控循环单元(GRU)网络的叠前时变地震子波提取方法.根据实际叠前地震数据分布特征与非平稳性质,本方法首先建立非平稳地震记录与添加随机噪声的时变子波训练数据集;为对提取出的时序特征进行拓展,提升传统GRU网络对长时序列的处理能力,本方法搭建起含多层GRU模块与全连接神经网络的改进门控循环单元网络模型;利用建立的训练数据集对网络模型进行训练使网络具备提取时变子波的能力;为提高训练效率与提取精度,本方法在训练的反向传播过程中应用自定义WaveLoss损失函数衡量误差,最终实现叠前时变子波的估计.经合成数据仿真实验与不同方法对比验证,本文提出的叠前时变子波提取方法具有更高的准确度;经对中国西部不同地区实际叠前地震资料处理与反褶积验证分析,该方法可有效提高目标区叠前地震剖面分辨率.

基于重加权L1的ATpV正则化叠前反演方法

地震叠前反演能够准确获取地下储层介质的各类参数,是油气的勘探与开发中重要技术之一。然而,地震反演是典型的病态问题,为了克服此问题,通常使用正则化约束目标函数,来减轻反演问题的病态性。但是正则化约束忽略了地层边界的振幅信息,使用重加权方法可以很好地克服这一问题,更好地恢复稀疏性。提出了一种基于重加权L1的ATpV正则化叠前三参数反演方法(ATpV-L1方法),首次将重加权L1方法与ATpV方法结合,并引入到叠前反演中。采用交替方向乘子算法(ADMM)建立反演框架,对目标函数进行分块优化,有效提高了收敛速度。首先,介绍ATpV-L1方法,建立了基于ATpV-L1的叠前反演目标函数;然后,应用理论模拟数据对比新方法和ATpV方法反演结果,验证了方法的效果;最后,使用实际数据进行实验分析,进一步验证了ATpV-L1方法的反演精度及可行性。实验结果表明,提出的ATpV-L1方法可以有效恢复反演结果的稀疏性,提高反演精度。

物理、数据先验认识融合的叠前解耦分步反演

AVA三参数反演在地层弹性参数预测中发挥着重要作用。由于AVA理论公式(即物理先验认识)中的角度不易确定,加之大型稀疏矩阵的病态性,导致常规叠前反演过程不稳定。为此,提出物理、数据先验认识融合的叠前解耦分步反演方法。首先,基于物理先验认识构建非稀疏正演框架,以增加参数反演的稳定性,为解耦分步反演奠定基础;然后,以井资料为数据先验认识,将物理、数据先验认识融合,对叠前地震数据进行解耦,以得到更准确的叠前地震属性数据;最后,对解耦后的叠前地震属性进行反演得到地层弹性参数。该方法通过井数据先验认识修正反演过程,可以避免因物理先验认识中角度不准带来的误差。实际数据测试结果表明,相比于业界三参数AVA反演方法,本方法的反演结果具有更高的精度,其中拉梅参数、剪切模量和密度的精度分别提高14.1%、13.6%和11.9%。

基于地震反演技术的地应力建模方法在江汉新沟嘴组页岩油中的应用

江汉盆地新沟嘴组页岩油地质评价资源量为3.6×108 t,具有良好的勘探潜力,其中地质评价Ⅰ类有利区陈沱口凹陷有多口井在新沟嘴组新下段Ⅱ油组泥质白云岩中获得油流。由于页岩油井通常采用“水平井+分段压裂”工艺实现经济开采,因而地应力研究对该区水平井储层压裂改造具有重要的意义。三维地应力建模研究的核心问题在于提升井间地应力预测效果,传统的有限元方法由于方法本身的限制,井间预测效果较差,有必要利用地震资料反演技术来提升井间地应力预测的效果。针对该区页岩层矿物组分多、局部发育硫酸盐岩、地应力条件复杂的问题,首先对新沟嘴组页岩层开展基于实验室测量值约束的精细测井评价,其结果与岩心全岩X衍射测量结果吻合;在此基础上,选择自相容近似地震岩石物理建模流程,对目标层的弹性参数进行了正演模拟,正演结果与实测曲线基本吻合。然后,基于地层压力、实验室测量的岩石力学性能参数和测井解释结果,依据不同岩性的动静态关系,采用有效应力比值和构造应变法计算井的地应力参数。最后运用基于地震叠前反演技术的三维地应力建模方法,使用与一维地质力学建模相同的相关性和参数约束,建立了新沟嘴组页岩油的三维地应力模型。与一维地质力学研究成果相比,该三维地应力模型提取的伪曲线与工区大部分井一维地应力结果基本吻合,验证了该方法在新沟嘴页岩油地应力评价中的适用性和准确性。

基于精确Zoeppritz方程的贝叶斯叠前地震随机反演

基于精确Zoeppritz方程的叠前地震反演方法在面向低信噪比地震资料的应用时仍然存在较大挑战。马尔科夫链蒙特卡洛(Markov chain Monte Carlo, MCMC)模拟是一种启发式的全局寻优算法,是实现叠前弹性参数非线性反演的有效途径。常规基于MCMC算法的叠前反演采用高斯分布来刻画弹性参数的统计特征,在应用于复杂岩性储层时有较大的局限性。同时,由于地下模型参数空间巨大以及地震数据中噪声等因素的影响,MCMC对弹性参数后验概率分布的搜索过程极易受到局部极值的影响,这使得基于MCMC的叠前反演较难获得稳定、准确的结果。本文针对实际复杂储层及低信噪比地震资料条件下基于精确Zoeppritz方程的叠前反演问题,提出了一种改进的MCMC弹性参数反演方法。该方法首先利用低频模型约束,将待反演参数转换为模型参数的扰动量,从而降低后验概率分布的复杂度;其次,通过对似然函数取对数,并利用低频模型来约束地震正演过程;最后,利用基于随机子空间采样的多链算法对叠前非线性反演问题进行全局寻优,以避免采样过程过早地收敛到局部极值。低信噪比模拟数据和实际数据的测试表明,本文所提方法能够获得更加准确、稳定的弹性参数反演结果,并且能够对反演结果给出可信、定量的不确定性估计。

Estimation of Q and inverse Q filtering for prestack reflected PP-and converted PS-waves

Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resolution is one of the key problems. The main factor affecting converted wave resolution is the absorption of seismic waves in overlying strata. In order to remove the effect of absorption on converted waves, inverse Q filtering is used to improve the resolution. In this paper, we present a method to estimate the S-wave Q values from prestack converted wave gathers. Furthermore, we extend a stable and effective poststack inverse Q filtering method to prestack data which uses wave field continuation along the ray path to compensate for attenuation in prestack common shot PP- and PS-waves. The results of theoretical modeling prove that the method of estimating the S-wave Q values has high precision. The results from synthetic and real data prove that the stable inverse Q filtering method can effectively improve the resolution of prestack PP- and PS-waves.

Zoeppritz equation-based prestack inversion and its application in fluid identification

Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in these approximations increases with increasing angle of incidence and variation of the elastic parameters, which increases the number of inversion solutions and minimizes the inversion accuracy. In this study, we explore a method for solving the reflection coefficients by using the Zoeppritz equations. To increase the accuracy of prestack inversion, the simultaneous inversion of P-wave velocity, S-wave velocity, and density by using prestack large-angle seismic data is proposed based on generalized linear inversion theory. Moreover, we reduce the ill posedness and increase the convergence of prestack inversion by using the regularization constraint damping factor and the conjugate gradient algorithm. The proposed prestack inversion method uses prestack large-angle seismic data to obtain accurate seismic elastic parameters that conform to prestack seismic data and are consistent with logging data from wells.

Study of prestack elastic parameter consistency inversion methods

The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.

Prestack inversion based on anisotropic Markov random field-maximum posterior probability inversion and its application to identify shale gas sweet spots

Economic shale gas production requires hydraulic fracture stimulation to increase the formation permeability. Hydraulic fracturing strongly depends on geomechanical parameters such as Young＇s modulus and Poisson＇s ratio. Fracture-prone sweet spots can be predicted by prestack inversion, which is an ill-posed problem; thus, regularization is needed to obtain unique and stable solutions. To characterize gas-bearing shale sedimentary bodies, elastic parameter variations are regarded as an anisotropic Markov random field. Bayesian statistics are adopted for transforming prestack inversion to the maximum posterior probability. Two energy functions for the lateral and vertical directions are used to describe the distribution, and the expectation-maximization algorithm is used to estimate the hyperparameters of the prior probability of elastic parameters. Finally, the inversion yields clear geological boundaries, high vertical resolution, and reasonable lateral continuity using the conjugate gradient method to minimize the objective function. Antinoise and imaging ability of the method were tested using synthetic and real data.

Detection of gas hydrate sediments using prestack seismic AVA inversion

Bottom-simulating reflectors (BSRs) in seismic profile always indicate the bottom of gas hydrate stability zone, but is difficult to determine the distribution and features of gas hydrate sediments (GHS). In this study, based on AVA forward modeling and angle-domain common-image gathers we use prestack AVA parameters consistency inversion in predicting gas hydrate sediments in the Shenhu area at northern slope of South China Sea, and obtain the vertical and lateral features and saturation of GHS.

VTI介质流体因子及各向异性参数直接反演方法

储层流体识别是油气勘探领域的重要研究方向。页岩储层具有较强的垂直横向各向同性(Transverse Isotropy with a Vertical Axis of Symmetry,VTI)介质特征,在流体因子估测中带来不容忽视的影响。为此,根据Rüger纵波反射系数方程,提出了一种基于VTI介质弹性阻抗的流体因子、杨氏模量及各向异性参数反演方法。首先,推导以Russell流体因子、杨氏模量和等效各向异性参数表示的反射系数近似公式和弹性阻抗方程;然后,通过弹性阻抗反演估计流体因子、杨氏模量和各向异性参数,并比较了各向异性参数对反演结果的影响。模型测试和实际工区数据应用结果表明,该方法可以合理、准确地预测流体因子、杨氏模量和各向异性参数,为页岩储层流体识别、脆性参数和各向异性发育程度预测提供了一种新的方法。

Attenuation compensation in multicomponent Gaussian beam prestack depth migration

Gaussian beam prestack depth migration is an accurate imaging method of subsurface media. Prestack depth migration of multicomponent seismic data improves the accuracy of imaging subsurface complex geological structures. Viscoelastic prestack depth migration is of practical significance because it considers the viscosity of the subsurface media. We use Gaussian beam migration to compensate for the attenuation in multicomponent seismic data. First, we use the Gaussian beam method to simulate the wave propagation in a viscoelastic medium and introduce the complex velocity Q-related and exact viscoelastic Zoeppritz equation. Second, we discuss PP- and PS-wave Gaussian beam prestack depth migration algorithms for common-shot gathers to derive expressions for the attenuation and compensation. The algorithms correct the amplitude attenuation and phase distortion caused by Q, and realize multicomponent Gaussian beam prestack depth migration based on the attenuation compensation and account for the effect of inaccurate Q on migration. Numerical modeling suggests that the imaging resolution of viscoelastic Gaussian beam prestack depth migration is high when the viscosity of the subsurface is considered.

High angle prestack depth migration with absorption compensation

The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the absorption effect and enhance the resolution. In this paper, we derive a one-way wave equation with an attenuation term based on the time- space domain high angle one-way wave equation. A complicated geological model is then designed and synthetic shot gathers are simulated with acoustic wave equations without and with an absorbing term. The derived one-way wave equation is applied to the migration of the synthetic gathers without and with attenuation compensation for the simulated shot gathers. Three migration profiles are obtained. The first and second profiles are from the shot gathers without and with attenuation using the migration method without compensation, the third one is from the shot gathers with attenuation using the migration method with compensation. The first and third profiles are almost the same, and the second profile is different from the others below the absorptive layers. The amplitudes of the interfaces below the absorptive layers are weak because of their absorption. This method is also applied to field data. It is concluded from the migration examples that the migration method discussed in this paper is feasible.

Anisotropic converted wave amplitude-preserving prestack time migration by the pseudo-offset method

In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.

基于无监督神经网络匹配算法的叠前表面多次波压制方法

为了有效地压制海上地震勘探数据中的表面多次波,实现勘探目标的正确成像,提出一种基于无监督神经网络的叠前表面多次波匹配算法,将神经网络方法与表面相关多次波压制方法相结合,通过设定学习率不断下降,用无监督神经网络取代匹配滤波算子,对叠前地震数据进行表面多次波的压制,既不需要传统的匹配算法,也不需要在标签数据集上进行训练。在简单合成数据、Sigsbee模型数据和实际数据上的应用结果验证了该方法对表面多次波压制的有效性。

相控叠前地质统计学反演在储层预测中的应用——以四川盆地川西坳陷须家河组二段致密砂岩储层为例

四川盆地川西坳陷新场构造带高庙子气田须家河组二段致密砂岩气藏多套砂体叠置,相变明显,表现为“整体富砂、局部含气”的特征,但有效储层发育段与泥岩的阻抗相当,常规叠后反演方法难以准确预测砂体边界及优质储层分布。在须二段地质特征研究基础上,首先开展有利相带边界识别及预测研究,并结合岩石物理敏感参数分析,开展相控约束下叠前地质统计学反演研究,实现对有效储层与泥岩层劈分,进而完成了相控有效储层预测。预测结果不仅提高了砂体识别精度,同时也提高了有效储层预测的可靠性,预测结果与实际钻井吻合度较高。该预测技术有效支撑了合兴场气田储量提交和探井井位部署,已产生了一定的油气效益。

准北地区侏罗系煤层脆性叠前预测

准噶尔盆地煤层气资源量超6万亿方,中石化探区1.3万亿方,资源潜力大;白家海地区彩探1H井,西山窑组煤层埋深2385 m,通过压裂日产气5.7万方,无水,一年累产气796万方,煤层气勘探取得了重大突破。准北地区石北凹陷煤层主要发育在西山窑组和八道湾组,西山窑组为主力层组.