一种适用于海相碳酸盐岩的快速古盐度计算方法与流程

背景技术：

2.海水盐度是海水中全部溶解固体与海水重量之比，通常以每千克海水中所含的克数表示，是人们用来表示海水中盐类物质的质量分数。现今世界大洋的平均盐度为35‰，最高为红海盐度，高达41‰，而最低为波罗的海盐度，为10‰。在地质历史时期中，海水的盐度是不断变化的，不同盐度下沉积形成的岩石类型变化很大，通过恢复沉积时期的盐度来判断岩石沉积环境，可以进一步分析可能的烃源岩或储层分布规律。3.恢复古盐度的方法很多，用地球化学方法(如b/ga、sr/ba、mg/ca、sr/ca和na/ca等)来重建不同沉积地层的古盐度在国内得到了较为广泛的应用。有研究利用对mn同位素值域、mn/sr值以及碳氧同位素交汇分析盐度。有研究利用地化mgo、cao、na2o、k2o等主量元素的含量和稀土元素的含量来判识古环境。有研究利用泥质岩的硼元素含量作为古盐度计，重建了湖泊沉积的古盐度变化曲线。另外，也有学者利用岩石脉体里的流体包裹体盐度测试恢复成岩过程中古流体盐度。4.目前与流体盐度的研究方法有关的研究仅有cn201711418008.x，其主要基于胶结物包裹体盐度，恢复后期流体盐度。其步骤如下：判断地层岩石胶结物次序，筛选发育多期胶结物且富存流体包裹体的样品；选择不同次序流体包裹体，进行均一温度和盐度测试，获得多个均一温度和盐度数据；结合埋藏史，确定地温演化与地质时间的对应关系，根据获得的多个流体包裹体的均一温度和盐度数据画出盐度随地质时间变化的趋势线，从而获得古流体盐度演化史。5.海相碳酸盐岩领域一直是深层油气勘探的主要阵地，我国近期突破的安岳、顺北、普光、元坝、川西等大型油气田都属于海相层系。而海相层系大油气田形成时的环境差异巨大，顺北奥陶系地层属于灰岩正常海水沉积，普光长兴组和安岳龙王庙组地层的形成环境为白云岩沉积，形成时海水略微咸化，靖边气田马家沟组和川西雷口坡组形成盐度更大。可见不同地层盐度差别较大，而在同一段地层中盐度变化也较大，与储层的关联性也不一样。通过盐度的分析，能够恢复地层沉积时海平面的变化过程，进而推测储层空间展布范围。而以往盐度分析需通过大量的岩石地球化学分析，时间和经济成本都较高。因此，目前亟待提出一种适用于海相碳酸盐岩的快速古盐度计算方法。

技术实现要素：

6.本发明的目的是针对现有技术的不足，提出一种适用于海相碳酸盐岩的快速古盐度计算方法。本发明的方法可以较好地应用于海相碳酸盐岩地层中。通过测井曲线计算碳酸盐岩地层沉积时的海水盐度，快速简便，不一定需要岩石地球化学分析，节约了时间和经济成本，为进一步分析烃源岩和储层分布规律奠定重要基础。在鄂尔多斯盆地已得到了较好的应用。7.为了实现上述目的，本发明提供一种适用于海相碳酸盐岩的快速古盐度计算方法，该方法包括如下步骤：8.s1：确定碳酸盐岩地层的多种岩性组合；9.s2：根据每种岩性组合对应的体积模型和测井曲线建立每种岩性组合中的各种岩性的体积百分含量计算公式；10.s3：利用各种岩性沉积时的海水盐度赋值和步骤s2建立的每种岩性组合中的各种岩性的体积百分含量计算公式建立每种岩性组合地层沉积时海水盐度计算公式，进而对碳酸盐岩地层沉积时的海水盐度进行计算。11.本发明的技术方案的有益效果如下：12.(1)本发明提供了一种适用于海相碳酸盐岩识别古海水盐度的方法，为进一步分析烃源岩和储层分布规律奠定重要基础。13.(2)本发明形成一套适用于海相碳酸盐岩的利用测井曲线快速定量识别古海水盐度的方法，解决了仅能依靠地球化学分析识别古海水盐度的问题，节约了时间和经济成本。14.(3)本发明的方法可操作性强，识别直观、清晰、简便，具有很高的应用价值，其方法步骤可以很方便地推广应用到海相碳酸盐岩地层的勘探开发中。利用本发明的方法还可以分段获得纵向上的盐度变化，进而推测不同时间段的海水环境变化。15.本发明的其它特征和优点将在随后具体实施方式部分予以详细说明。具体实施方式16.下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本发明更加透彻和完整，并且能够将本发明的范围完整地传达给本领域的技术人员。17.本发明提供一种适用于海相碳酸盐岩的快速古盐度计算方法，该方法包括如下步骤：18.s1：确定碳酸盐岩地层的多种岩性组合；19.s2：根据每种岩性组合对应的体积模型和测井曲线建立每种岩性组合中的各种岩性的体积百分含量计算公式；20.s3：利用各种岩性沉积时的海水盐度赋值和步骤s2建立的每种岩性组合中的各种岩性的体积百分含量计算公式建立每种岩性组合地层沉积时海水盐度计算公式，进而对碳酸盐岩地层沉积时的海水盐度进行计算。21.在本发明中，碳酸盐岩地层中的孔隙在地层中所占体积较小，因此忽略孔隙体积的影响。22.根据本发明，优选地，在步骤s1中，对所述碳酸盐岩地层的沉积微相和成岩作用进行分析，进而确定所述碳酸盐岩地层的多种岩性组合。23.根据本发明，优选地，所述多种岩性组合包括白云岩+灰岩+泥质组合、白云岩+灰岩+硬石膏+泥质组合和白云岩+硬石膏+盐岩+泥质组合。24.根据本发明，优选地，所述白云岩+灰岩+泥质组合在单位体积地层的矿物组成包括白云石、方解石和泥质。25.根据本发明，优选地，所述白云岩+灰岩+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0026]vdol1＝2.242*ρb1-0.142*u1-4.077ꢀꢀꢀꢀꢀꢀꢀ(1)；[0027]vlim1＝1.132*ρb1+0.261*u1-5.753ꢀꢀꢀꢀꢀꢀꢀ(2)；[0028]vsh1＝-3.375*ρb1-0.119*u1+10.830ꢀꢀꢀꢀꢀꢀꢀ(3)；[0029]其中，vdol1、vlim1、vsh1分别为所述白云岩+灰岩+泥质组合中的白云石、方解石和泥质在单位体积地层中的体积含量；[0030]ρb1为密度测井值；[0031]u1为光电吸收截面测井值；[0032]所述白云岩+灰岩+泥质组合对应的体积模型为vdol1+vlim1+vsh1＝1；[0033]公式(1)、公式(2)和公式(3)中各项前的系数是通过根据所述所述白云岩+灰岩+泥质组合对应的体积模型建立的所述白云岩+灰岩+泥质组合对应的密度测井曲线和光电吸收截面测井曲线与所述白云岩+灰岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0034]在本发明中，所述白云岩+灰岩+泥质组合对应的体积模型为单位体积地层中所述白云岩+灰岩+泥质组合的各矿物的体积含量之和为1，即：vdol1+vlim1+vsh1＝1。通过该体积模型建立的“密度测井曲线和光电吸收截面测井曲线”与“所述白云岩+灰岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量”之间的关系为：所述密度测井值ρb1和所述光电吸收截面测井值u1分别为上述三种矿物的密度测井值和光电吸收截面测井值与各对应矿物体积含量的乘积之和，再通过利用所述白云岩+灰岩+泥质组合对应的体积模型，进而求得公式(1)、公式(2)和公式(3)中各项前的系数。[0035]根据本发明，优选地，所述白云岩+灰岩+硬石膏+泥质组合在单位体积地层的矿物组成包括白云石、方解石、硬石膏和泥质。[0036]根据本发明，优选地，所述白云岩+灰岩+硬石膏+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0037]vdol2＝0.715*ρb2-0.188*u2-0.003*dt1+1.235ꢀꢀꢀꢀꢀꢀꢀ(4)；[0038]vlim2＝-3.609*ρb2+0.116*u2-0.010*dt1+10.742ꢀꢀꢀꢀꢀꢀꢀ(5)；[0039]vanhy1＝3.120*ρb2+0.095*u2+0.007*dt1-10.855ꢀꢀꢀꢀꢀꢀꢀ(6)；[0040]vsh2＝-0.226*ρb2-0.023*u2+0.007*dt1-0.122ꢀꢀꢀꢀꢀꢀꢀ(7)；[0041]其中，vdol2、vlim2、vanhy1、vsh2分别为所述白云岩+灰岩+硬石膏+泥质组合中的白云石、方解石、硬石膏和泥质在单位体积地层中的体积含量；[0042]ρb2为密度测井值；[0043]u2为光电吸收截面测井值；[0044]dt1为声波测井值；[0045]所述白云岩+灰岩+硬石膏+泥质组合对应的体积模型为vdol2+vlim2+vanhy1+vsh2＝1；[0046]公式(4)、公式(5)、公式(6)和公式(7)中各项前的系数是通过根据所述白云岩+灰岩+硬石膏+泥质组合对应的体积模型建立的所述白云岩+灰岩+硬石膏+泥质组合对应的密度测井曲线、光电吸收截面测井曲线和声波测井曲线与所述白云岩+灰岩+硬石膏+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0047]在本发明中，所述白云岩+灰岩+硬石膏+泥质组合对应的体积模型为单位体积地层中所包含各矿物的体积含量之和为1，即：vdol2+vlim2+vanhy1+vsh2＝1；通过该体积模型建立的“密度测井曲线、光电吸收截面测井曲线和声波测井曲线”与“所述白云岩+灰岩+硬石膏+泥质组合中的各个矿物组成在单位体积地层中的体积含量”之间的关系为：密度测井值、光电吸收截面测井值和声波测井值分别为上述四种矿物的密度测井值、光电吸收截面测井值、声波测井值与各对应矿物体积含量的乘积之和，再通过利用所述白云岩+灰岩+硬石膏+泥质组合对应的体积模型，进而求得公式(4)、公式(5)、公式(6)和公式(7)中各项前的系数。[0048]根据本发明，优选地，所述白云岩+硬石膏+盐岩+泥质组合在单位体积地层的矿物组成包括白云石、硬石膏、石盐和泥质。[0049]根据本发明，优选地，所述白云岩+硬石膏+盐岩+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0050]vdol3＝0.946*ρb3-0.196*u3-0.003*dt2+0.547ꢀꢀꢀꢀꢀꢀꢀ(8)；[0051]vanhy2＝-0.265*ρb3+0.204*u3-0.003*dt2-0.779ꢀꢀꢀꢀꢀꢀꢀ(9)；[0052]vhal＝-1.647*ρb3+0.053*u3-0.005*dt2+4.903ꢀꢀꢀꢀꢀꢀꢀ(10)；[0053]vsh3＝0.966*ρb3-0.061*u3+0.010*dt2-3.672ꢀꢀꢀꢀꢀꢀꢀ(11)；[0054]其中，vdol3、vanhy2、vhal、vsh3分别为所述白云岩+硬石膏+盐岩+泥质组合中的白云石、硬石膏、石盐和泥质在单位体积地层中的体积含量；[0055]ρb3为密度测井值；[0056]u3为光电吸收截面测井值；[0057]dt2为声波测井值；[0058]所述白云岩+硬石膏+盐岩+泥质组合对应的体积模型为vdol3+vanhy2+vhal+vsh3＝1；[0059]公式(8)、公式(9)、公式(10)和公式(11)中各项前的系数是通过根据所述白云岩+硬石膏+盐岩+泥质组合对应的体积模型建立的所述白云岩+硬石膏+盐岩+泥质组合对应的密度测井曲线、光电吸收截面测井曲线和声波测井曲线与所述白云岩+硬石膏+盐岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0060]在本发明中，所述白云岩+硬石膏+盐岩+泥质组合对应的体积模型为单位体积地层中所包含各矿物的体积含量之和为1，即：vdol3+vanhy2+vhal+vsh3＝1；通过该体积模型建立的“密度测井曲线、光电吸收截面测井曲线和声波测井曲线”与“所述白云岩+硬石膏+盐岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量”之间的关系为：密度测井值ρb3、光电吸收截面测井值u3和声波测井值dt2分别为上述四种矿物的密度测井值、光电吸收截面测井值、声波测井值与各对应矿物体积含量的乘积之和，再通过利用所述白云岩+硬石膏+盐岩+泥质组合对应的体积模型，进而求得公式(8)、公式(9)、公式(10)和公式(11)中各项前的系数。[0061]根据本发明，优选地，[0062]所述白云岩+灰岩+泥质组合地层沉积时海水盐度计算公式为：[0063]s1＝(k*vdol1+l*vlim1)/(vdol1+vlim1)ꢀꢀꢀꢀꢀꢀ(12)；[0064]所述白云岩+灰岩+硬石膏+泥质组合地层沉积时海水盐度计算公式为：[0065]s2＝(k*vdol2+l*vlim2+m*vanhy1)/(vdol2+vanhy1+vlim2)ꢀꢀꢀꢀꢀꢀ(13)；[0066]所述白云岩+硬石膏+盐岩+泥质组合地层沉积时海水盐度计算公式为：[0067]s3＝(k*vdol3+n*vhal+m*vanhy2)/(vdol3+vanhy2+vhal)ꢀꢀꢀꢀꢀꢀ(14)；[0068]其中，vdol1、vlim1分别为所述白云岩+灰岩+泥质组合中的白云石、方解石在单位体积地层中的体积含量；[0069]vdol2、vlim2、vanhy1分别为所述白云岩+灰岩+硬石膏+泥质组合中的白云石、方解石、硬石膏在单位体积地层中的体积含量；[0070]vdol3、vanhy2、vhal分别为所述白云岩+硬石膏+盐岩+泥质组合中的白云石、硬石膏、石盐在单位体积地层中的体积含量；[0071]k为白云岩沉积时的海水盐度赋值，依据现今白云岩沉积环境确定，在本发明中，白云岩沉淀的德克萨斯和科威特的潟湖和潮坪中盐度为42-44‰，因此，k取值优选为43‰。[0072]l为灰岩沉积时的海水盐度赋值，依据现今灰岩沉积环境确定，在本发明中，灰岩沉积时盐度与现今海水大致相近，因此，l取值优选为35‰。[0073]m为硬石膏沉积时的海水盐度赋值，依据现今硬石膏沉积环境确定，在本发明中，硬石膏发生沉淀时的盐度要明显高于灰岩和白云岩沉积时的盐度，因此，m取值优选为100‰。[0074]n为盐岩沉积时的海水盐度赋值，依据现今盐岩沉积环境确定，在本发明中，盐岩沉淀时的盐度可高达260‰，因此，n取值优选为260‰。[0075]s1为碳酸盐岩地层的白云岩+灰岩+泥质组合沉积时的海水盐度；[0076]s2为碳酸盐岩地层的白云岩+灰岩+硬石膏+泥质组合沉积时的海水盐度；[0077]s3为碳酸盐岩地层的白云岩+硬石膏+盐岩+泥质组合沉积时的海水盐度。[0078]以下通过实施例具体说明本发明。[0079]实施例1[0080]本实施例提供一种适用于海相碳酸盐岩的快速古盐度计算方法，计算鄂尔多斯盆地某区块的某井深3056.6m的碳酸盐岩地层沉积时的盐度。[0081]该方法包括如下步骤：[0082]s1：确定该井的碳酸盐岩地层的岩性组合为白云岩+灰岩+泥质组合；所述白云岩+灰岩+泥质组合在单位体积地层的矿物组成包括白云石、方解石和泥质。[0083]s2：根据白云岩+灰岩+泥质组合中的白云岩、灰岩、泥质岩性对应的体积模型和测井曲线建立所述白云岩+灰岩+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0084]vdol1＝2.242*ρb1-0.142*u1-4.077ꢀꢀꢀꢀꢀꢀꢀ(1)；[0085]vlim1＝1.132*ρb1+0.261*u1-5.753ꢀꢀꢀꢀꢀꢀꢀ(2)；[0086]vsh1＝-3.375*ρb1-0.119*u1+10.830ꢀꢀꢀꢀꢀꢀꢀ(3)；[0087]其中，vdol1、vlim1、vsh1分别为所述白云岩+灰岩+泥质组合中的白云石、方解石和泥质在单位体积地层中的体积含量；[0088]ρb1为该井在这一深度的密度测井值为2.810g/cm3；[0089]u1为该井在这一深度的光电吸收截面测井值为11.286b/cm3；[0090]所述白云岩+灰岩+泥质组合对应的体积模型为vdol1+vlim1+vsh1＝1；[0091]公式(1)、公式(2)和公式(3)中各项前的系数是通过根据所述各种岩性对应的体积模型建立的密度测井曲线和光电吸收截面测井曲线与所述白云岩+灰岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0092]由此可得：vdol1为62.6％、vlim1为31.2％、vsh1为0.2％；[0093]s3：本实施例中，利用白云岩沉积时的海水盐度赋值、灰岩沉积时的海水盐度赋值和步骤s2建立的白云岩+灰岩+泥质组合中的各种岩性的体积百分含量计算公式(1)、公式(2)和公式(3)建立白云岩+灰岩+泥质组合地层沉积时海水盐度计算公式为：[0094]s1＝(k*vdol1+l*vlim1)/(vdol1+vlim1)ꢀꢀꢀꢀꢀꢀꢀ(12)；[0095]其中，vdol1为62.6％，vlim1为31.2％；[0096]k为白云岩沉积时的海水盐度赋值，依据现今白云岩沉积环境确定为43‰；[0097]l为灰岩沉积时的海水盐度赋值，依据现今灰岩沉积环境确定为35‰；[0098]进而计算得到碳酸盐岩地层的白云岩+灰岩+泥质组合沉积时的海水盐度为40.4‰。[0099]实施例2[0100]本实施例提供一种适用于海相碳酸盐岩的快速古盐度计算方法，计算鄂尔多斯盆地某区块的某井深3097.1m的碳酸盐岩地层沉积时的盐度。[0101]该方法包括如下步骤：[0102]s1：确定该井的碳酸盐岩地层的岩性组合为白云岩+灰岩+硬石膏+泥质组合；所述白云岩+灰岩+硬石膏+泥质组合在单位体积地层的矿物组成包括白云石、方解石、硬石膏和泥质。[0103]s2：根据白云岩+灰岩+硬石膏+泥质组合中的白云岩、灰岩、硬石膏、泥质岩性对应的体积模型和测井曲线建立所述白云岩+灰岩+硬石膏+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0104]vdol2＝0.715*ρb2-0.188*u2-0.003*dt1+1.235ꢀꢀꢀꢀꢀꢀꢀ(4)；[0105]vlim2＝-3.609*ρb2+0.116*u2-0.010*dt1+10.742ꢀꢀꢀꢀꢀꢀꢀ(5)；[0106]vanhy1＝3.120*ρb2+0.095*u2+0.007*dt1-10.855ꢀꢀꢀꢀꢀꢀꢀ(6)；[0107]vsh2＝-0.226*ρb2-0.023*u2+0.007*dt1-0.122ꢀꢀꢀꢀꢀꢀꢀ(7)；[0108]其中，vdol2、vlim2、vanhy1、vsh2分别为所述白云岩+灰岩+硬石膏+泥质组合中的白云石、方解石、硬石膏和泥质在单位体积地层中的体积含量；[0109]ρb2为该井在这一深度的密度测井值为2.825g/cm3；[0110]u2为该井在这一深度的光电吸收截面测井值为10223b/cm3；[0111]dt1为该井在这一深度的声波测井值为154.768μs/m；[0112]所述白云岩+灰岩+硬石膏+泥质组合对应的体积模型为vdol2+vlim2+vanhy1+vsh2＝1；[0113]公式(4)、公式(5)、公式(6)和公式(7)中各项前的系数是通过根据所述各种岩性对应的体积模型建立的密度测井曲线、光电吸收截面测井曲线和声波测井曲线与所述白云岩+灰岩+硬石膏+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0114]由此可得：vdol2为81.0％、vlim2为10.8％、vanhy1为0.2％、vsh2为8.0％；[0115]s3：本实施例中，利用白云岩沉积时的海水盐度赋值k、灰岩沉积时的海水盐度赋值l、硬石膏沉积时的海水盐度赋值m和步骤s2建立的白云岩+灰岩+硬石膏+泥质组合中的各种岩性的体积百分含量计算公式(4)、公式(5)、公式(6)和公式(7)建立白云岩+灰岩+硬石膏+泥质组合地层沉积时海水盐度计算公式为：[0116]s2＝(k*vdol2+l*vlim2+m*vanhy1)/(vdol2+vanhy1+vlim2)ꢀꢀꢀꢀꢀꢀ(13)；[0117]其中，vdol2为81.0％、vlim2为10.8％、vanhy1为0.2％；[0118]k为白云岩沉积时的海水盐度赋值，依据现今白云岩沉积环境确定为43‰；[0119]l为灰岩沉积时的海水盐度赋值，依据现今灰岩沉积环境确定为35‰；[0120]m为硬石膏沉积时的海水盐度赋值，依据现今硬石膏沉积环境确定为100‰，[0121]进而计算得到碳酸盐岩地层的白云岩+灰岩+硬石膏+泥质组合沉积时的海水盐度为42.2‰。[0122]实施例3[0123]本实施例提供一种适用于海相碳酸盐岩的快速古盐度计算方法，计算鄂尔多斯盆地某区块的某井深3661.9m的碳酸盐岩地层沉积时的盐度。[0124]该方法包括如下步骤：[0125]s1：确定该井的碳酸盐岩地层的岩性组合为白云岩+硬石膏+盐岩+泥质组合；所述白云岩+硬石膏+盐岩+泥质组合在单位体积地层的矿物组成包括白云石、硬石膏、石盐和泥质。[0126]s2：根据白云岩+硬石膏+盐岩+泥质组合中的白云岩、硬石膏、盐岩、泥质岩性对应的体积模型和测井曲线建立所述白云岩+硬石膏+盐岩+泥质组合中的各种岩性的体积百分含量计算公式分别为：[0127]vdol3＝0.946*ρb3-0.196*u3-0.003*dt2+0.547ꢀꢀꢀꢀꢀꢀꢀ(8)；[0128]vanhy2＝-0.265*ρb3+0.204*u3-0.003*dt2-0.779ꢀꢀꢀꢀꢀꢀꢀ(9)；[0129]vhal＝-1.647*ρb3+0.053*u3-0.005*dt2+4.903ꢀꢀꢀꢀꢀꢀꢀ(10)；[0130]vsh3＝0.966*ρb3-0.061*u3+0.010*dt2-3.672ꢀꢀꢀꢀꢀꢀꢀ(11)；[0131]其中，vdol3、vanhy2、vhal、vsh3分别为所述白云岩+硬石膏+盐岩+泥质组合中的白云石、硬石膏、石盐和泥质在单位体积地层中的体积含量；[0132]ρb3为该井在这一深度的密度测井值为2.730g/cm3；[0133]u3为该井在这一深度的光电吸收截面测井值为11.187b/cm3；[0134]dt2为该井在这一深度的声波测井值为194.093μs/m；[0135]所述白云岩+硬石膏+盐岩+泥质组合对应的体积模型为vdol3+vanhy2+vhal+vsh3＝1；[0136]公式(8)、公式(9)、公式(10)和公式(11)中各项前的系数是通过根据所述各种岩性对应的体积模型建立的密度测井曲线、光电吸收截面测井曲线和声波测井曲线与所述白云岩+硬石膏+盐岩+泥质组合中的各个矿物组成在单位体积地层中的体积含量之间的关系，求解得出的。[0137]由此可得：vdol3为41.8％、vanhy2为21.0％、vhal为6.9％、vsh3为30.3％；[0138]s3：本实施例中，利用白云岩沉积时的海水盐度赋值、硬石膏沉积时的海水盐度赋值、盐岩沉积时的海水盐度赋值和步骤s2建立的白云岩+硬石膏+盐岩+泥质组合中的各种岩性的体积百分含量计算公式(8)、公式(9)、公式(10)和公式(11)建立白云岩+硬石膏+盐岩+泥质组合地层沉积时海水盐度计算公式为：[0139]s3＝(k*vdol3+n*vhal+m*vanhy2)/(vdol3+vanhy2+vhal)ꢀꢀꢀꢀꢀꢀ(14)；[0140]其中，：vdol3为41.8％、vanhy2为21.0％、vhal为6.9％；[0141]k为白云岩沉积时的海水盐度赋值，依据现今白云岩沉积环境确定为43‰；[0142]m为硬石膏沉积时的海水盐度赋值，依据现今硬石膏沉积环境确定为100‰；[0143]n为盐岩沉积时的海水盐度赋值，依据现今盐岩沉积环境确定为260‰。[0144]进而计算得到碳酸盐岩地层的白云岩+硬石膏+盐岩+泥质组合沉积时的海水盐度为81.7‰。[0145]通过实施例1-3可知：利用本发明的方法还可以分段获得纵向上的盐度变化，进而推测不同时间段的海水环境变化。[0146]以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。