CN105510993A - 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 - Google Patents

Abstract

本发明提供了一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法。该方法包括：根据露头数据、钻井数据和测井数据中的岩性垂向变化，保留因沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，建立原始岩性序列，通过该序列恢复研究区的古地貌模型，并恢复现今构造模型；开展三维重磁电震联合反演，依靠原始岩性序列和现今构造模型，在地层框架内通过电法、地震相和密度异常，开展钻井周围地区膏盐岩层的岩性岩相解释，确定岩性序列平面分布规律。本发明提供的方法通过地质建模搞清前陆盆地膏盐岩层变形特征的基础上，融合三维重磁电震等地球物理方法开展联合反演，能够准确确定变形后膏盐岩层的分布特点。

Description

前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法

技术领域

本发明涉及一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法，属于石油开采与地质勘探技术领域。

背景技术

膏盐岩是水体遭受蒸发、盐分逐渐浓缩后沉淀形成的化学成因沉积岩，又称“蒸发岩”。在我国沉积盆地内，寒武纪、奥陶纪、三叠纪、白垩纪和古近纪遍布了大量的膏盐岩，聚集较多时可形成多种沉积矿产(如钾盐矿)，可广泛应用于化工、医药、农业、玻璃、炸药等领域，是一种紧缺的矿物资源，对国家经济建设和资源安全具有重要意义。另外，膏盐岩一般是非渗透层，许多大型、特大型油气田往往都是以巨厚膏盐岩层作盖层，故膏盐岩层对油气藏的聚集和保存有重要意义。但是，在油气勘探和开发的钻井过程中，由于膏盐岩地层易溶易变，经常发生井壁垮塌、井漏和卡钻事故，因此，准确识别和预测膏盐岩层对保障井场施工安全也很重要。

现有的膏盐岩分布预测方法主要利用地表工作、岩心钻探、水化学分析、物探测井等技术，基于已提出的膏盐岩沉积理论模型，进行盐岩识别和分布预测。

根据现有的膏盐岩沉积理论模型，膏盐岩层往往与碳酸盐岩共存，且具有明显的平面分带性，常呈牛眼状分布，从盆地边缘向中心的岩性依次为石灰岩、白云岩、膏岩、石盐及钾镁盐岩等；膏盐岩偶与碎屑岩共存，其岩性变化为砂砾岩、泥岩、泥灰岩、膏盐岩。这与盐水浓缩过程中矿物结晶顺序密切相关，但仅能用于识别和解释蒸发岩所在岩石层序的岩性序列，对于含盐盆地内部膏盐岩的分布位置有指导作用，但也仅能根据岩性序列指出其分布方向，无法对其分布做出准确预测。

地表工作方法主要通过调查盐泉或者膏盐岩露头点，寻找亲盐岩石的指示矿物如天青石、菱镁矿、自形石英及其他卤素沉积物；寻找盐溶产物(例如盐溶角砾岩、岩溶假晶)；详细测量膏盐岩地层的野外剖面，分析其沉积相特征及变化，确定可能的膏盐岩层位。岩心钻探方法主要通过钻井，直接收集各井点的膏盐岩岩石标本，确定膏盐岩层的岩性和厚度。水化学分析方法是通过分析地下水的基本要素(矿化度、水质类型)及其与成矿有关的元素富集程度(钾、溴、氯、镁等)，来判断该地区是否发育膏盐岩层的一种化学分析方法。这些地质类方法可以准确识别膏盐岩层的岩石类型和岩性序列，但只是从单个地质点出发，无法从平面的角度精细预测膏盐岩的分布范围。

物探测井技术包括电法、重力法、地震法及放射性测井法等几种技术。电法主要分析卤水的电阻值判断岩性，具体包括电剖面法和电测物探法，前者可以识别岩性，后者主要用来预测范围。重力法是根据膏岩或盐岩地层密度特点，重力测量结果的异常来确定盐层的埋深、厚度和分布范围。地震法主要利用膏盐岩层的地震反射特征来探测埋藏较浅、范围不大的盐丘构造。这几种技术对于膏盐岩的识别和预测也具有各自的局限性，例如地震法对于埋深较大的膏盐岩层识别准确率较低，重力法受到地层密度的影响较大，电法勘探受到卤水浓度的影响，在对范围不大、岩性简单和岩层变形不强烈的膏盐岩层发育地区，这两种方法比较实用。但是如果遇到原始沉积是岩性序列就较为复杂、后期在前陆构造背景下膏盐岩受构造挤压发生了蠕变、流动、揉皱、重复和叠置时，例如前陆盆地背景下的膏盐岩层，上述这些物探测井方法就无能为力了，无法做到精细识别和预测。

目前而言，现有的地质理论模型对膏盐岩的识别比较有效，因为通过地表工作和岩心钻探可以直接获得膏盐岩层的岩石类型、岩性垂向变化序列和膏盐岩层厚度，从而预测膏盐岩的分布边界。例如通过露头区的岩性序列和盆地覆盖区的岩性序列进行对比，判断井点或露头点的岩性分布特征，根据膏盐岩环带分布特征加以判断该井点或露头点所处在环带哪个位置，从而判断膏盐岩的分布范围。但是这种方法仅能判断出哪些范围内有膏盐岩，但是具体厚度无法刻画，且仅针对于简单的膏盐岩层或者岩性较为单一的膏盐岩丘，对于前陆盆地背景下的挤压型复杂膏盐岩层而言，这种方法显得力不从心。

这种方法的缺点主要表现在以下几点。首先，地表工作依靠露头区的出露岩石，无法将岩性直接应用于盆地覆盖区，需要根据岩性序列进行推断；这种推断在理想模型下(即环带分布模式)是可以的，但是在强烈构造挤压后膏盐岩层岩性序列发生重复或反转的前陆盆地是行不通的，因为膏盐岩的厚度已经发生了变化，井点的变化可以通过已钻井来揭示，但井间的这种变化是无法判断的。

其次，受到钻井取心成本的控制，岩心钻探在盆地覆盖区得到的只是少量井点、部分井段的岩性序列，也不足以反映全井段或者全盆地的岩性分布。例如在塔里木盆地北部某坳陷总面积约2.8×10⁴km²，仅有探井40余口，且分布不均，根本无法控制全区膏盐岩的平面分布，更何况受到强烈构造挤压后，膏盐岩的厚度变化在平面上并非线性趋势变化，因此现有技术无法满足膏盐岩层的精细刻画和预测。

现有的物探测井方法对膏盐岩的分布预测有两个方面的技术，一个是通过测井资料进行解释和识别膏盐岩层，第二个是根据膏盐岩的岩石物理特征，通过地震资料或者电法资料对膏盐岩层平面分布特征进行反演预测，例如采用测井约束稀疏脉冲反演方法，即通过岩石物理建模计算岩石物理弹性参数，确定井上的膏岩、盐岩波阻抗范围，利用Jason软件计算三维空间的波阻抗体，根据膏岩、盐岩的波阻抗范围圈定膏盐岩的展布，其中最重要的参数就是依靠地震速度谱(参见CN101556337A“一种确定地下深层特殊岩性体的方法”)。

上述方法在膏盐岩层变化较小的地区可以适用，例如东部拉张背景的断陷盆地、中部稳定的克拉通盆地内的膏盐岩变形往往没有那么强烈，而且有变形往往是正断层的缺失或者是简单的褶皱穹窿，而对于西部前陆盆地挤压型的膏盐岩层来说，这种简单的地球物理反演方法则有些捉襟见肘。首先，前陆盆地构造挤压往往形成多个冲断带，膏盐岩层会随着冲断带的叠加而发生垂向逆冲叠加。因此，这种叠加导致垂向上可以出现多个相互间隔的膏盐岩层，利用上述现有的物探测井方法是无法准确预测的。其次，山前地震资料信噪比低，所依靠的地震速度谱资料失真，而电法探测的精确度较低，没有很好的层的概念，难以反映深层岩性的准确变化。

因此，研发出一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法，仍是本领域亟待解决的问题之一。

发明内容

为解决上述技术问题，本发明的目的在于提供一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法。该方法能够实现对前陆盆地背景下形成的受到强烈构造挤压发生变形的复杂膏盐岩层内部岩性序列的古环境恢复以及变形后的现今膏盐岩层的精细刻画和预测。

为达到上述目的，本发明提供了一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法，其包括以下步骤：

(1)获取前陆盆地深埋挤压型复杂膏盐岩层的露头数据和/或测井数据和/或钻井数据以及三维重磁电震数据；

(2)通过对比露头数据的岩性垂向变化规律、和/或对钻井数据和/或测井数据中的岩性垂向变化进行判别，保留因为沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，最终建立研究区膏盐岩的原始岩性序列(即原始沉积序列)，通过该序列恢复研究区的古地貌模型，并恢复现今的构造模型；

(3)用钻井数据标定相同深度的三维重磁电震数据，开展三维重磁电震联合反演，依靠步骤(2)中得到的膏盐岩的原始岩性序列和现今的构造模型，在地层框架内通过电法和/或地震相和/或密度异常，开展钻井周围地区膏盐岩层的岩性岩相解释，识别岩性，圈定各岩性的边界，确定岩性序列平面分布规律(包括厚度及边界)，完成前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测。

在本发明中，主要涉及的技术术语的定义如下所述。前陆盆地：位于造山带前缘与相邻克拉通之间的沉积盆地；由岩石圈挤压挠曲沉降作用而形成的沉积盆地，发育于聚合、碰撞或挤压环境；一般为位于造山带与克拉通之间的山前拗陷；具有与挤压有关的、从造山带-前渊-前隆的不对称的构造-地层格架。膏岩：蒸发岩的一种，主要由石膏和硬石膏组成的岩石。盐岩：蒸发岩的一种，由可溶盐类组成的化学成因沉积岩的统称，主要由石盐组成。深埋挤压型复杂膏盐岩层：前陆盆地背景下，盆地容易受到挤压应力场的影响而发生的一种构造变形，同时会影响到盆内的膏盐岩层，因为膏盐岩层在深埋地下时(4000m～6000m)具有极强的塑性，可以发生蠕变揉皱，故构造挤压可以使其岩性序列发生叠置或反转，以及膏盐岩层的各种厚度变化(增厚或减薄)。

在上述方法中，优选地，在步骤(1)中，获取露头数据的步骤包括：在研究区内膏盐岩层露头出露地区，利用工具(包括皮尺、放大镜、锤头等)测量及描述膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化。

在上述方法中，优选地，在步骤(1)中，获取钻井数据的步骤包括：根据钻井过程中返回的岩屑确定该井点处膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化。

在上述方法中，优选地，在步骤(1)中，获取测井数据的步骤包括：放置测井仪在井筒内开展自上而下的测试，收集反映岩性及其垂向变化的测井数据。

在上述方法中，优选地，在步骤(1)中，获取三维重磁电震数据的步骤包括：设计网格以足够覆盖研究区，然后在网格节点处放置震源和电源激发器和接收器，通过激发器释放震源和电源，在接收器同时接收地下返回的信息。这些信息通过计算机程序(如中国石油集团公司的Geoeast软件)处理后，可以直接加载在Landmark工作站上使用。

在上述方法中，优选地，在步骤(1)之后，还可以包括步骤(1)-1：从地震剖面(属于地震数据)选择研究区没有构造变形或构造变形较弱地区的钻井，根据该钻井的测井数据和/或钻井岩性特征，通过对比相同深度的伽马曲线，电阻率曲线，声波曲线，钻时曲线和井径曲线中的一种或几种的数值和形态，识别研究区膏盐岩的岩性，并划分岩性段(包括亚段)，作为最终结果的一部分。上述的地震剖面、伽马曲线，电阻率曲线，声波曲线，钻时曲线和井径曲线均可以由本领域常规的方法得到。其中，地震剖面可以是通过以下步骤获得的：将地震数据加载至Landmark工作站，通过SeisWork模块提取常规地震剖面；所述地震数据可以是通过以下步骤获得的：设计网格以足够覆盖研究区，然后在网格节点处放置震源激发器和接收器，通过激发器释放震源，在接收器同时接收地下返回的信息。

在上述方法中，优选地，步骤(2)具体可以包括：通过对比露头数据的岩性垂向变化规律、和/或对钻井数据和/或测井数据中的岩性垂向变化进行判别，并且对比钻井处地震数据(该地震数据是通过本领域的常规方法获得，例如可以参考文献：梁向豪等，2011，库车大北构造带三维叠前深度偏移处理解释技术，中国石油勘探，5：2-8；江民等，2015，库车山前带高陡盐下构造三维采集技术及效果，天然气勘探与开发，38(3)：35-38；杨红霞，2003，地震数据采集技术进展，勘探地球物理进展，26(5)：6-11；陈小二，2010，复杂山地地震采集技术在库车坳陷的应用，地质勘探，30(9)：25-27)中的地震相，判断此处是否有断裂或者地层变形；结合利用2D-Move软件开展的构造建模，计算出构造变形中挤压增厚或者断层重复增厚引起的岩性变化和厚度变化及其形成的岩性序列剖面的变化；保留因为沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，最终建立研究区膏盐岩的原始岩性序列(即原始沉积序列)和原始沉积剖面，通过该序列和剖面恢复研究区的古地貌模型，并恢复现今的构造模型。其中，恢复古地貌模型可以是通过拉平膏盐岩层原始岩性序列及沉积剖面的顶界面，以恢复古地貌模型。恢复现今的构造模型可以是通过对现今膏盐岩层的变形特征、断裂和褶皱类型进行抽象化，保留变形机制，去除重复变形，最终建立构造模型(利用2D-Move软件)。

在上述方法中，优选地，在步骤(3)中，开展三维重磁电震联合反演是利用Landmark工作站中的Seiswork模块进行的。

在上述方法中，优选地，在步骤(3)中，在地层框架内通过电法和/或地震相和/或密度异常开展钻井周围地区膏盐岩层的岩性岩相解释具体可以包括：寻找重磁密度异常和/或电法电阻率异常高的数据范围，以此作为依据圈定膏盐岩层的范围和厚度，和/或比对不同构造部位(依据步骤(2)中得到的现今的构造模型)的地震相反射特征，以开展钻井周围地区膏盐岩层的岩性岩相解释。其中，电法测得的电阻率，以及地震法得到的地震相(地震数据)均是本领域的常规技术手段，例如可以参考文献：郑晓月，2009，电法勘探中数据采集设备和计算机通信的实现，陕西理工学院学报(自然科学版)；陈儒军等，2003，高精度多频电法采集系统，物探与化探；何继善，1997，电法勘探的发展和展望，地球物理学报。

本发明提供了一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法。该方法完善了传统膏盐岩地质理论模型，并融合了地质与地球物理技术，综合地质建模、构造建模和重磁电联合反演，实现了对受到强烈构造挤压发生变形的复杂膏盐岩层内部岩性序列的古环境恢复以及变形后的现今膏盐岩层的精细刻画和预测。

与现有技术相比，本发明的技术方案主要具有以下优点：

(1)相比前人单独利用地质(地表工作、水化学分析、岩心钻探)或地球物理资料(如地震相、速度谱)进行膏盐岩层预测，本发明的方法融合了地质建模(沉积和构造模型)与地球物理技术，开展综合分析，提高了适用于前陆盆地复杂膏盐岩层的预测精度。

(2)本发明的方法弥补了地震资料、电法和重力资料在各自解决问题时的不足，融合形成三维重磁电震联合反演技术。

(3)本发明的方法增加了对膏盐岩层受到挤压后发生变形的形变分析，提出了一套如何消除构造挤压与断层重复引起膏盐岩变形的方法，恢复了原始沉积序列和古地貌。

综上所述，本发明提供的方法通过地质建模搞清前陆盆地膏盐岩层变形特征的基础上，融合三维重磁电震等地球物理方法开展联合反演，能够准确确定变形后膏盐岩层的分布特点，包括古环境特征、岩性序列组合与分布特征、以及现今膏盐岩厚度分布规律。该方法主要用于前陆盆地深埋挤压型复杂膏盐岩层的识别和预测，对膏盐岩层开展精细刻画，确定前陆盆地膏盐岩层的分布范围、厚度变化，预测各部位的岩性组合序列，解决制约膏盐岩层引起的一系列速度异常、钻井工程安全和其他油气勘探开发的难题，最大限度的减少安全事故、提高圈闭落实精度，从而提高探井成功率。在我国西部如塔里木盆地，一口井的钻井成本在一亿元左右，这一方法的提出为提高含膏盐岩的前陆盆地钻探成功率奠定了坚实基础，为我国石油工业带来巨大经济效益。

附图说明

图1为实施例中的前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法流程图。

图2为实施例中的五曲线识别膏盐岩层的各种岩性图。

图3为实施例中的挤压型复杂膏盐岩层的岩性段精细对比图。

图4为实施例中的复杂膏盐岩层内部的岩性序列重建图。

图5为实施例中的研究区膏盐岩层原始沉积剖面。

图6a和图6b为实施例中的研究区膏盐岩层平面分布图。

具体实施方式

为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。

塔里木盆地某工区内，古近系发育了3000m以上的巨厚膏盐岩层，受到新近纪以来强烈的南北向构造挤压影响，埋藏较深的膏盐岩层具有较大塑性，发生蠕变、岩性反转和聚集，其岩性序列分布和岩石厚度分布难以预测，成为了地层层速度模拟和钻井工程安全的重要难题。以此前陆盆地古近系膏盐岩层为例，对本发明的技术方案作进一步详细说明。

本实施例提供了一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法，其包括以下步骤，如图1所示：

(1)获取前陆盆地深埋挤压型复杂膏盐岩层的露头数据和/或测井数据和/或钻井数据以及三维重磁电震数据；

获取露头数据的步骤包括：在研究区内膏盐岩层露头出露地区，利用工具(包括皮尺、放大镜、锤头等)测量及描述膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化；

获取钻井数据的步骤包括：根据钻井过程中返回的岩屑确定该井点处膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化；

获取测井数据的步骤包括：放置测井仪在井筒内开展自上而下的测试，收集反映岩性及其垂向变化的测井数据；

获取三维重磁电震数据的步骤包括：设计网格以足够覆盖研究区，然后在网格节点处放置震源和电源激发器和接收器，通过激发器释放震源和电源，在接收器同时接收地下返回的信息，这些信息通过计算机程序(如中国石油集团公司的Geoeast软件)处理后，可以直接加载在landmark工作站上使用。

(2)泥膏盐岩的精细识别与岩性段的逐级精细对比

从地震剖面选择没有构造变形或构造变形弱地区的钻井，根据该钻井的测井数据和钻井岩性特征，通过对比相同深度的伽马曲线，电阻率曲线，声波曲线，钻时曲线和井径曲线(即五曲线，通过常规方法获得)的数值和形态，准确识别和区分泥岩、膏岩、盐岩及过渡岩性。五曲线识别膏盐岩层的各种岩性图如图2所示。其中，地震剖面可以是通过以下步骤获得的：将地震数据加载至Landmark工作站，通过SeisWork模块提取常规地震剖面；所述地震数据可以是通过以下步骤获得的：设计网格以足够覆盖研究区，然后在网格节点处放置震源激发器和接收器，通过激发器释放震源，在接收器同时接收地下返回的信息。

在岩性识别的基础上，通过云灰岩标志层控制开展岩性段与亚段精细划分，将该套膏盐岩层所在层序可划分为泥岩段、膏盐岩段和膏泥岩段三个岩性段，以第一套碳酸盐岩为标志，结合碎屑岩和化学岩的纵横向发育特征，将膏泥岩段划分为3-5个亚段。该挤压型复杂膏盐岩层的岩性段精细对比图如图3所示。

(3)消除构造挤压及断层重复的影响并恢复原始沉积序列及古地貌

膏盐岩段在构造挤压下发生盐岩流动增厚和断层重复增厚，而构造位置的细微差异都会造成岩性序列和地层厚度的差异。

首先，通过对比钻井数据中的岩石类型、厚度及其垂向变化，并且对比钻井处地震数据中的地震相(通过常规方法获得)，判断此处是否有断裂或者地层变形；其次，结合利用2D-Move软件开展的构造建模，计算出构造变形中挤压增厚或者断层重复增厚引起的岩性变化和厚度变化及其形成的岩性序列剖面的变化，得到复杂膏盐岩层内部的岩性序列重建图如图4所示(现今构造模型也由图4体现)；最后，保留钻井数据中因为沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，从而恢复从盆地边缘到盆地中心的原始膏盐岩沉积序列，建立原始膏盐岩沉积剖面，如图5所示，通过拉平膏盐岩层原始沉积序列及沉积剖面的顶界面恢复古地貌(图5)。

(4)开展三维重磁电震联合反演

用钻井数据标定相同深度的三维重磁电震数据，利用landmark工作站中的Seiswork模块开展三维重磁电震联合反演，寻找计算结果中的重磁密度异常、电法电阻率异常高的数据范围，以此作为依据圈定膏盐岩层的范围和厚度；再根据步骤(2)中得到的膏盐岩沉积剖面和现今的构造模型，比对不同构造部位的地震相反射特征，开展钻井周围地区膏盐岩层的岩性岩相解释，识别岩性，圈定各岩性的边界，建立现今岩性序列分布规律，精细刻画不同构造位置膏盐岩层的沉积厚度及岩性序列，得到的研究区膏盐岩层平面分布图(含不同部位沉积序列)如图6a和图6b所示。其中，电法测得的电阻率，以及地震法得到的地震相(地震数据)均通过常规方法获得。

Claims (9)

1.一种前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法，其包括以下步骤：

(1)获取前陆盆地深埋挤压型复杂膏盐岩层的露头数据和/或测井数据和/或钻井数据以及三维重磁电震数据；

(2)通过对比露头数据的岩性垂向变化规律、和/或对钻井数据和/或测井数据中的岩性垂向变化进行判别，保留因为沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，最终建立研究区膏盐岩的原始岩性序列，通过该序列恢复研究区的古地貌模型，并恢复现今的构造模型；

(3)用钻井数据标定相同深度的三维重磁电震数据，开展三维重磁电震联合反演，依靠步骤(2)中得到的膏盐岩的原始岩性序列和现今的构造模型，在地层框架内通过电法和/或地震相和/或密度异常，开展钻井周围地区膏盐岩层的岩性岩相解释，识别岩性，圈定各岩性的边界，确定岩性序列平面分布规律，完成前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测。

2.根据权利要求1所述的方法，其中，在步骤(1)中，获取露头数据的步骤包括：在研究区内膏盐岩层露头出露地区，利用工具测量及描述膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化。

3.根据权利要求1所述的方法，其中，在步骤(1)中，获取钻井数据的步骤包括：根据钻井过程中返回的岩屑确定该井点处膏盐岩层的岩石类型和/或颜色和/或厚度及这些数据的垂向变化。

4.根据权利要求1所述的方法，其中，在步骤(1)中，获取测井数据的步骤包括：放置测井仪在井筒内开展自上而下的测试，收集反映岩性及其垂向变化的测井数据。

5.根据权利要求1所述的方法，其中，在步骤(1)中，获取三维重磁电震数据的步骤包括：设计网格以足够覆盖研究区，然后在网格节点处放置震源和电源激发器和接收器，通过激发器释放震源和电源，在接收器同时接收地下返回的信息。

6.根据权利要求1所述的方法，其在步骤(1)之后，还包括步骤(1)-1：从地震剖面选择研究区没有构造变形或构造变形较弱地区的钻井，根据该钻井的测井数据和/或钻井岩性特征，通过对比相同深度的伽马曲线，电阻率曲线，声波曲线，钻时曲线和井径曲线中的一种或几种的数值和形态，识别研究区膏盐岩的岩性，并划分岩性段。

7.根据权利要求1所述的方法，其中，步骤(2)包括：通过对比露头数据的岩性垂向变化规律、和/或对钻井数据和/或测井数据中的岩性垂向变化进行判别，并且对比钻井处地震数据中的地震相，判断此处是否有断裂或者地层变形；结合构造建模，计算出构造变形中挤压增厚或者断层重复增厚引起的岩性变化和厚度变化及其形成的岩性序列剖面的变化；保留因为沉积环境变化引起的岩性序列，消除构造挤压与断层重复引起的膏盐岩层变形和叠置的序列，最终建立研究区膏盐岩的原始岩性序列和原始沉积剖面，通过该序列和剖面恢复研究区的古地貌模型，并恢复现今的构造模型。

8.根据权利要求1所述的方法，其中，在步骤(3)中，开展三维重磁电震联合反演是利用landmark工作站中的Seiswork模块进行的。

9.根据权利要求1所述的方法，其中，在步骤(3)中，在地层框架内通过电法和/或地震相和/或密度异常开展钻井周围地区膏盐岩层的岩性岩相解释包括：寻找重磁密度异常和/或电法电阻率异常高的数据范围，以此作为依据圈定膏盐岩层的范围和厚度，和/或比对不同构造部位的地震相反射特征，以开展钻井周围地区膏盐岩层的岩性岩相解释。