CN115793053B - 一种基于地震、钻井与录井的寒武系膏盐岩综合预测方法 - Google Patents

Abstract

本发明提供了一种基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，是首先对地震资料进行保幅处理，再基于保幅处理后的地震资料开展叠前弹性参数反演获取纵波阻抗数据，预测寒武系膏盐岩目标层深度范围并制定针对性的钻探方案，然后开展寒武系膏盐岩钻探参数综合变化率分析，计算寒武系膏盐岩钻探参数综合变化率，依据寒武系膏盐岩钻探参数综合变化率进一步校准膏盐岩顶面位置，最后在地震、钻井预测寒武系膏盐岩深度范围开展录井参数岩屑观察和X射线元素分析，实现寒武系膏盐岩综合预测。本发明构思合理，过程简单，去除了不同入射角度子波的影响，实现了膏盐岩层综合预测，降低了钻探风险，提高了油气层保护程度。

Description

一种基于地震、钻井与录井的寒武系膏盐岩综合预测方法

技术领域

本发明涉及油气勘探技术领域，具体涉及一种基于地震、钻井与录井的寒武系膏盐岩综合预测方法。

背景技术

膏盐岩主要产于气候干旱的内陆盆地盐湖中或被砂坝与海隔绝的浅水泻湖和干涸的海湾中，是盆地内水体遭受蒸发、盐分逐渐浓缩以至发生沉淀而形成。Usiglio(1848)的实验表明当海水蒸发时，可溶盐是按溶解度由小至大的顺序依次沉淀形成蒸发矿物的。海水在浓缩过程中，溶解度最小的碳酸盐首先沉淀，海水继续蒸发浓缩，浓度达到15％-17％时，石膏类矿物开始析出，海水浓度为26％时，石盐开始结晶。

膏盐岩地层中的沉积岩除了盐岩和膏岩外，常还有与其共生的碳酸盐岩系，如灰岩、灰质白云岩、泥灰岩等。膏盐岩常较致密、塑性好，可作为封存油气的优质盖层(金之钧等，2010)。全球120个含油气盆地中约58％的油气田与膏盐岩有关，从前寒武纪到第四纪各时期的含油气盆地中均可见油-盐“共生”的现象，其中46％的油气层位于膏盐岩层之下，13％的油气层位于膏盐岩层之间。中国已发现的165个天然气田中，受膏盐岩控制的约占40％、天然气总储量约占36％(胡安平，2019)。

塔里木盆地早寒武世经历了短暂的深水盆地环快演变成以云坪为主的局限台地环境，发育下寒武统优质海相烃源岩和丘滩相碳酸盐岩储层；中寒武世海平面下降，气候炎热干燥，海水咸化程度高，发育了蒸发湖占主导的强镶边型台地，并形成中寒武统膏岩和盐岩相间、厚度不一、普遍分布的盖层(图8)。这套生储盖组合分布面积达约24万km²，油气潜力巨大，是油气增储上产的重要勘探对象。

图9为塔里木盆地西北部多口寒武系盐下膏盐岩及其背景碳酸盐岩的纵波阻抗对比。可以看出，云岩和灰岩的纵波阻抗值分别为2×10⁷kg/m³*m/s和1.74×10⁷kg/m³*m/s，相对较高，盐岩的纵波阻抗值仅为9.23×10⁶kg/m³*m/s，可显著区分于其他岩性，但硬石膏的纵波阻抗值为1.82×10⁷kg/m³*m/s，与背景碳酸盐岩存在显著的数据分布范围重叠，单一依据纵波阻抗数据难于识别膏盐岩。此外，常规叠后地震资料反演获得的波阻抗数据受限于AVO效应的影响，其精度较低。

尽管可作为优质油气盖层，但膏盐岩地层在钻井过程中经常发生塑性蠕变，盐岩易垮塌，膏岩、泥岩膨胀缩径，盐上、盐间和盐下地层压力差异大，常常导致复杂工程事故，井控风险大。针对盐膏层，钻探工程必须采用专打专封开次，钻探过程中需用特殊钻探方式配合高密度泥浆体系，钻后采用高强度套管予以固井封隔，因而对盐膏层的层位、深度、厚度预测等要尽可能准确。如果对盐膏顶卡层失误，低密度钻井液情况下误揭开盐膏层，可导致瞬间缩径卡钻，导致侧钻；如果盐膏底卡层失误，高密度钻井液情况下揭开下部低压地层，可导致瞬间失返性漏失，井筒泥浆压力降低，极可能造成盐膏层缩径卡钻、岩屑下沉卡钻、溢流等井控风险发生。2016-2022年间，塔里木盆地西北部针对寒武系盐下目标、因为膏盐岩层预测失准而最终回填侧钻的钻井高达8口，占比90％以上。

塔里木盆地中寒武统膏盐岩层发育段无固定沉积模式，膏盐层顶底均为白云岩，缺乏标志层，且横向上单层和累计厚度变化较大，预测十分困难。

由于膏盐岩具有“低温脆变、高温塑变”的特征，特别是盐岩岩石特征表现为低速(约220us/m)、低密度(2030kg/m³)，与灰岩(约156us/m和2710kg/m³)和云岩(约143us/m和2870kg/m³)有较大差异，因而常用地震波阻抗反演结果对盐岩发育对予以预测。但膏岩的速度和密度(硬石膏约164us/m和2980kg/m³；石膏约174us/m和2350kg/m³)常与背景碳酸盐岩有重叠，难于依照波阻抗反演予以区分(图2)。此外，地震资料分辨率不高，传统叠后波阻抗反演受限于AVO效应影响，其预测精度有限(Zhang et al.,2011)。

膏盐岩层钻探施工中，针对录井岩屑开展特殊的化学试剂实验、元素分析(X射线荧光分析、XRF)等也可对岩性进行预测，但XRF实验分析结果受岩屑选择影响，而常用欠饱和泥浆体系下盐岩段岩屑难于获取。此外，XRF分析耗时较长，全井段分析成本较大。

发明内容

针对上述背景技术中存在的技术问题，本发明提出了一种构思合理，过程简单，去除了不同入射角度子波的影响，构建了膏盐层钻探参数综合变化率计算公式，实现了膏盐岩层综合预测，降低了钻探风险，提高了油气层保护程度的基于地震、钻井与录井的寒武系膏盐岩综合预测方法。

为解决上述技术问题，本发明提供的一种基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，是首先对地震资料进行保幅处理，再基于保幅处理后的地震资料开展叠前弹性参数反演获取纵波阻抗数据，预测寒武系膏盐岩目标层深度范围并制定针对性的钻探方案，然后开展寒武系膏盐岩钻探参数综合变化率分析，计算寒武系膏盐岩钻探参数综合变化率，依据寒武系膏盐岩钻探参数综合变化率进一步校准膏盐岩顶面位置，最后在地震、钻井预测寒武系膏盐岩深度范围开展录井参数岩屑观察和X射线元素分析，实现寒武系膏盐岩综合预测。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中，所述纵波阻抗数据的计算公式为：

其中，Rpp为地震反射系数，反射界面上下介质的纵横波速度和密度分别为v_p1、v_p2、v_s1、v_s2、ρ₁、ρ₂，θ₁、θ₂分别是反射角和透射角，则v_p＝(v_p1+v_p2)/2，

ρ＝(ρ₁+ρ₂)/2，Δv_p＝v_p2-v_p1，Δv_s＝v_s2-v_s1，Δρ＝ρ₂-ρ₁，θ＝(θ₁+θ₂)/2，θ₂＝arcsin(sinθ₁/v_p1*v_p2)，v_p*ρ为纵波阻抗。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中：所述叠前弹性参数反演获取纵波阻抗数据中纵波阻抗数值小于1.45×10⁷kg/m³*m/s的范围为预测的膏盐岩发育层段。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中：所述寒武系膏盐岩钻探参数综合变化率的计算公式为：

其中，Dc(t₀)为t₀深度的钻探参数综合变化率；DT(t₀)、Dp(t₀)、Or(t₀)分别为t₀深度处的钻时、钻压和出口电导率，DT(t₀-1)、Dp(t₀-1)、Or(t₀-1)分别为t₀前一深度处的钻时、钻压和出口电导率，a、b和c为权重系数，且a+b+c＝1；Dc(0)＝DT(0)＝Dp(0)＝Or(0)＝0。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中：所述权重系数a、b和c中默认a＝b＝c＝1/3，在纵波阻抗预测可疑深度范围，钻探参数变化率首次大于3的深度视为膏盐岩顶面位置。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中，所述寒武系膏盐岩综合预测的策略为：(1)挑选欠饱和盐水钻井液体系下岩屑，若无岩屑返出，直接判断为盐岩；(2)依据欠饱和盐水钻井液体系下挑选的岩屑开展X射线元素分析。

所述基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其中，所述步骤(2)的具体方法为：首先确定该区目标地层泥岩含量，然后按粘土、白云石、黄铁矿、石膏、方解石、硅质的顺序进行计算定量矿物含量，确立岩性组成。

采用上述技术方案，本发明具有如下有益效果：

本发明基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法构思合理，过程简单，首先基于保幅地震资料开展叠前弹性参数反演获得纵波阻抗，去除AVO效应对传统叠后波阻抗反演的影响，并依据反演的纵波阻抗数据预测膏盐岩层深度段大致范围，指导钻探方案编制；而后在随钻过程中开展钻探参数综合变化率分析，进一步确定膏盐岩层顶界面，并开展针对性的录井技术分析；再后在膏盐岩发育层段开展岩屑观察和X射线元素分析，实现了膏盐岩层综合预测，降低了钻探风险，提高了油气层保护程度。

本发明提出的钻探参数综合变化率可以快速反映地下地层可钻性、复杂性和电离子含量变化信息，支撑膏盐岩顶界面确立和录井分析技术优选，建立的基于岩屑X射线元素分析的膏盐岩层预测策略可以提升盐岩判断效率，并大幅提升膏盐岩等不同矿物含量预测准确程度。

本发明采用保幅处理地震资料进行叠前弹性参数反演获得纵波阻抗，对求取的纵波阻抗，相比较于传统叠后波阻抗反演去除了不同入射角度子波的影响。

附图说明

为了更清楚地说明本发明具体实施方式或现有技术中的技术方案下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法的流程图；

图2为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法中过新苏参1井(设计井)的地震剖面图；

图3为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法中图3所示剖面钻前地震波阻抗反演结果示意图；

图4为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法中新苏参1井所在地震道反演预测膏盐岩段(e)及其随钻过程中出口电导率(a)、钻时(b)、钻压(c)和钻探参数综合变化率(d)对比示意图；

图5为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法的实施例中钻探参数变化率(a)、地震反演预测(b)、X射线元素录井参数(c-g)和基于地震、钻井和录井综合预测的寒武系岩性分布(h)示意图；

图6为传统地震预测膏盐岩(a)、钻后测井解释膏盐岩(b)、两者误差对比(c-d)密度反演结果对比示意图；

图7为本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法中综合预测膏盐岩(a)、钻后测井解释膏盐岩(b)、两者误差对比(c-d)密度反演结果对比示意图；

图8为塔里木盆地中下寒武统地层柱状图；

图9为塔里木盆地中下寒武统典型岩性纵波阻抗分布图。

具体实施方式

下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

下面结合具体的实施方式对本发明做进一步的解释说明。

如图1所示，本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法,是首先对地震资料进行保幅处理，抽取高品质共反射点道集，采用叠前弹性参数反演获取纵波阻抗数据；基于叠前反演波阻抗数据，预测中寒武统膏盐岩发育层段深度范围并制定针对性的钻探方案；然后开展出口电导率、钻压、钻时等钻探参数综合变化率分析，计算钻探参数综合变化率，进一步校准中寒武统膏盐岩顶面位置；最后在地震、钻井预测膏盐岩深度范围开展录井参数岩屑观察及X射线元素分析，实现寒武系膏盐层综合预测。

其中,采用保幅地震资料进行叠前弹性参数反演获得纵波阻抗，相比较于传统叠后波阻抗反演去除了不同入射角度子波的影响。叠前弹性参数反演中纵波阻抗数值小于1.45×10⁷kg/m³*m/s的范围为预测的膏盐岩发育层段；对求取的纵波阻抗相比较于传统叠后波阻抗反演去除了不同入射角度子波的影响。

如图1所示，本发明基于地震、钻井与录井的寒武系膏盐岩综合预测方法,具体流程为：

S110、叠前弹性参数反演

Zoeppritz方程是进行AVO叠前弹性反演的理论基础，但是由于方程的过于复杂，通常人们使用其近似方程用于各种计算。以波阻抗反射系数作为参数反演，可以使反演问题更稳定为起点，按照波阻抗反射系数对Aki和Richards方程进行重新整理，给出了以波阻抗反射系数表示的近似方程，即：

其中，Rpp为地震反射系数，反射界面上下介质的纵横波速度和密度分别为v_p1、v_p2、v_s1、v_s2、ρ₁、ρ₂，θ₁，θ₂分别是反射角和透射角，则v_p＝(v_p1+v_p2)/2，

ρ＝(ρ₁+ρ₂)/2，Δv_p＝v_p2-v_p1，Δv_s＝v_s2-v_s1，Δρ＝ρ₂-ρ₁，θ＝(θ₁+θ₂)/2，θ₂＝arcsin(sinθ₁/v_p1*v_p2)；

分别代表纵横波反射系数和密度梯度；A(θ)＝sec²θ，/>

采用对反射系数的近似方法即下公式(2)，可将反射系数的表达式(即公式(2)左侧第一个等式)直接转化为地层的弹性参数，实现弹性参数的离散化：

根据上述公式(2)，可以将连续多界面的纵波AVO反射系数方程近似表示为下式：

通过式(3)建立线性反演方程组，加入子波矩阵W以考虑地震资料的带限特征及调谐效应，综合噪音的影响，则可得到如下K个界面N个入射角度的单道反演公式，

即，d_KN*1＝G_KN*3Km_3K*1+n (4)

其中，d为KN行实际地震观测数据；m为3K行待求反演的参数(纵波阻抗，横波阻抗和密度)；G为KN行3K列的正演算子，由子波矩阵W，一阶差分矩阵D和系数矩阵组成；n为KN行地震数据包含的噪音；L_p＝Ln(I_p)，L_s＝Ln(I_s)，L_D＝Ln(I_D)，

鉴于此，针对保幅处理(即强调在资料处理过程中保持有效信号相对振幅和相位信息不受破坏，尽可能地保持地震波振幅与地层界面反射系数成合理的比例关系，它既包括垂向上不同地层以及横向上同一地层的相对振幅关系，也包括同一反射点振幅随炮检距的振幅关系)后的地震道集(d)，估算子波(w)后，即可依据上式(5)进行叠前弹性参数反演求解获得纵波阻抗(Ip)数据，并依据反演获得的阻抗数据在考虑频散效应影响下进行膏盐岩层段预测。

图2为过塔里木盆地西北部柯坪地区过新苏参1井(设计井)的全叠加地震资料，图3为统一剖面位置处叠前弹性参数反演获得的纵波阻抗数据，井旁曲线为反演结果提取的拟阻抗曲线。可以看出，叠前弹性参数反演获得的纵波阻抗直接指示了中寒武统膏盐岩层段发育的大体区域(纵波阻抗<1.45×10⁷kg/m³*m/s)，以此可以为钻探工程的方案编制提供深度直接的依据。

S120、膏盐层钻探参数综合变化率识别

依据叠前弹性参数反演获知的膏盐岩层段发育深度段，可以制定针对性的钻探方案，包括井身结构、钻机、钻探参数优选、泥浆体系的配置等。一般针对膏盐层发育段，首先应简化钻具(扩孔钻头)，避免钻进中蠕变造成卡钻；其次应在膏盐层段采用双心钻头(比如8 1/2″-9 1/2″)，保证井径扩大率，保障套管顺利下入；第三钻探采用短起下钻和“进一退二”措施，加强短起及划眼，加强坐岗；第四应尽量饱和盐水体系控制盐层溶解，并选择较高的泥浆比重(1.85g/cm³-2.0g/cm³)，延缓膏盐层蠕变，若盐层因结晶卡钻，立即泵入低氯钻井液解卡。

一般钻头在碳酸盐岩硬地层背景下进入膏盐岩层段时，其钻时普遍相对降低。膏盐岩具有较大的蠕变特性、低渗透性、损伤恢复等特性，钻探过程中扭矩会明显波动，返出钻井液电导率增大，钻压变小。为了在叠前弹性参数反演获知纵波阻抗预测深度基础上进一步弄准膏盐岩层顶界面，本发明提出基于钻探参数综合变化率识别策略。具体依照纵波阻抗预测膏盐岩层深度周缘一定范围加强钻探参数监控与分析，包括钻压、钻时和出口电导率等参数，钻探参数综合变化率计算公式为：

其中，Dc(t₀)为t₀深度的钻探参数综合变化率；DT(t₀)、Dp(t₀)、Or(t₀)分别为t₀深度处的钻时、钻压和出口电导率，DT(t₀-1)、Dp(t₀-1)、Or(t₀-1)分别为t₀前一深度处的钻时、钻压和出口电导率，a、b和c为权重系数，且a+b+c＝1；Dc(0)＝DT(0)＝Dp(0)＝Or(0)＝0。

图4d为依据新苏参1井钻时、钻压和出口电导率计算的钻探参数综合变化率结果，其中权系数a＝b＝c＝1/3。显然，在深度约3935m位置处，出口电导率剧烈增加，钻时与钻压均有显著降低并持续抖动，综合变化率出现一个从背景值约0.2％突然变化至6.8％的异常高峰，指示进入疑似膏盐岩发育段，反演纵波阻抗预测的结果(图4e)也有指示。鉴于此，从3935m深度上浮50m起开展X射线元素录井分析。

S130、膏盐层底面及其综合预测

常规录井一般包括岩屑岩心、荧光和综合录井等，膏盐岩层段可增加相对昂贵的X射线元素录井予以预测岩性(参考表1)。

膏盐岩地层通常采用欠饱和盐水钻井液，当从碳酸盐岩地层钻进至膏盐岩地层时，盐岩被溶解，井底返至地面的岩屑量有明显的减少。钻进膏岩(CaSO₄)时，其岩屑样品常呈灰白色，用手可碾碎为粉末状，不与稀盐酸反应，与热盐酸溶解后，加氯化钡溶液会产生硫酸钡白色沉淀，S元素会明显增加。而钻进盐岩(NaCl)时，欠饱和盐水钻井液体系下常常难有岩屑反出，且钻井液中氯离子含量显著上升。

中寒武统膏盐岩层段之下为吾松格尔和肖尔布拉克组云岩储层，是重要的油气富集层段。当钻头穿过膏盐岩层段进入下伏碳酸盐岩地层，录井气测值会上升，甚至发生井漏或溢流，岩屑成分逐渐变为碳酸盐岩，Mg或Ca元素会显著上升。

表1碳酸盐岩剖面中常见矿物的主要元素含量表

鉴于此，本发明提出欠饱和盐水钻井液体系下，目标层段基于岩屑X射线元素分析的寒武系膏盐岩综合预测的策略为：

S131、挑选欠饱和盐水钻井液体系下岩屑，若无岩屑返出，直接判断为盐岩。

S132、依据欠饱和盐水钻井液体系下挑选的岩屑开展X射线元素分析，首先确定该区目标地层泥岩(泥质)含量，然后按粘土、白云石、黄铁矿、石膏、方解石、硅质的顺序进行计算定量矿物含量，确立岩性组成。各种矿物的定量计算公式(参考表1)如下：

1)粘土含量计算公式：

X_iclay≈Al_i/9.961×100％

其中，X_iclay为i深度处粘土含量，Al_i为i深度处Al元素含量。

2)白云石含量计算公式：

X_idolo≈(Mg_i–1.212*X_iclay)/13％

其中，X_idolo为i深度处白云石含量，Mg_i为i深度处Mg元素含量。

3)石膏含量计算公式：

X_igyps≈[Si–53.3％(Fe_i–5.3155*Al_i/9.961)/46.7％]/23.5％

其中，X_igyps为i深度处石膏含量，Si为i深度处S元素含量，Fe_i为i深度处Fe元素含量。

4)方解石含量计算公式：

X_ilime≈[Ca_i–21.7％Xidolo–29.4％Xigyps]/40％

其中，X_ilime为i深度处方解石含量，Ca_i为i深度处Ca元素含量。

图5c-g为新苏参1井反演纵波阻抗和钻探参数确立可疑深度范围元素录井结果，以及最终解释的岩性柱状图(图5h)对比。显然，综合预测的岩性种类(6种)和较纵波阻抗(仅1种)更加丰富，且深度范围的分辨率显著提高。在3935m深度附近，反演纵波阻抗与钻探参数综合变化率均有异常，综合解释依据元素分析，进一步确切其矿物含量，属于含膏泥质云岩。在4300-4710m深度范围内，反演纵波阻抗数值显著降低，且综合参数变化率异常丰富，综合解释其为膏盐岩集中发育段，进一步确立了8套盐岩。

S140膏盐层预测结果及其误差对比

图6为传统地震预测膏盐岩(a)、钻后测井解释膏盐岩(b)、两者误差对比(c-d)密度反演结果对比示意图。显然，叠前弹性参数反演的波阻抗分辨能力相对有限，只能依据相对低的纵波阻抗分布确立膏盐岩层段发育大体位置，且无法具体区分岩性，与钻后测井解释结果存在较大误差。

图7为本发明的综合预测膏盐岩(a)、钻后测井解释膏盐岩(b)、两者误差对比(c-d)密度反演结果对比示意图。可以看出，本发明综合预测的膏盐岩与最终钻后测井解释的结果形态基本一致，存在的由于岩屑掉块、捞取不准的较小误差，整体可靠，可以支撑钻探方案优化，降低钻探工程风险，实现油气层保护。

本发明提出的钻探参数综合变化率可以快速反映地下地层可钻性、复杂性和电离子含量变化信息，支撑膏盐岩顶界面确立和录井分析技术优选，建立的基于岩屑X射线元素分析的膏盐岩层预测策略可以提升盐岩判断效率，并大幅提升膏盐岩等不同矿物含量预测准确程度，从而降低钻探风险，保幅油气层。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (5)

1.一种基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其特征在于：首先对地震资料进行保幅处理，再基于保幅处理后的地震资料开展叠前弹性参数反演获取纵波阻抗数据，预测寒武系膏盐岩目标层深度范围并制定针对性的钻探方案，然后开展寒武系膏盐岩钻探参数综合变化率分析，计算寒武系膏盐岩钻探参数综合变化率，依据寒武系膏盐岩钻探参数综合变化率进一步校准膏盐岩顶面位置，最后在地震、钻井预测寒武系膏盐岩深度范围开展录井参数岩屑观察和X射线元素分析，实现寒武系膏盐岩综合预测；

所述寒武系膏盐岩钻探参数综合变化率的计算公式为：

其中上式(6)中Dc(t₀)为t₀深度的钻探参数综合变化率；DT(t₀)、Dp(t₀)、Or(t₀)分别为t₀深度处的钻时、钻压和出口电导率，DT(t₀-1)、Dp(t₀-1)、Or(t₀-1)分别为t₀前一深度处的钻时、钻压和出口电导率，a、b和c为权重系数，且a+b+c＝1；Dc(0)＝DT(0)＝Dp(0)＝Or(0)＝0；

所述权重系数a、b和c中默认a＝b＝c＝1/3，在纵波阻抗预测可疑深度范围，钻探参数变化率首次大于3的深度视为膏盐岩顶面位置。

2.如权利要求1所述的基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其特征在于，所述纵波阻抗数据的计算公式为：

其中，上式(1)中Rpp为地震反射系数，反射界面上下介质的纵横波速度和密度分别为v_p1、v_p2、v_s1、v_s2、ρ₁、ρ₂，θ₁、θ₂分别是反射角和透射角，则v_p＝(v_p1+v_p2)/2，v_s＝(v_s1+v_s2)/2，ρ＝(ρ₁+ρ₂)/2，Δv_p＝v_p2-v_p1，Δv_s＝v_s2-v_s1，Δρ＝ρ₂-ρ₁，θ＝(θ₁+θ₂)/2，θ₂＝arcsin(sinθ₁/v_p1*v_p2)，v_p*ρ为纵波阻抗。

3.如权利要求1所述的基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其特征在于：所述叠前弹性参数反演获取纵波阻抗数据中纵波阻抗数值小于1.45×10⁷kg/m³*m/s的范围为预测的膏盐岩发育层段。

4.如权利要求1所述的基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其特征在于，所述寒武系膏盐岩综合预测的策略为：

(1)挑选欠饱和盐水钻井液体系下岩屑，若无岩屑返出，直接判断为盐岩；

(2)依据欠饱和盐水钻井液体系下挑选的岩屑开展X射线元素分析。

5.如权利要求4所述的基于地震、钻井与录井参数的寒武系膏盐岩综合预测方法，其特征在于，所述步骤(2)的具体方法为：首先确定该区目标地层泥岩含量，然后按粘土、白云石、黄铁矿、石膏、方解石、硅质的顺序进行计算定量矿物含量，确立岩性组成。

Images