具体实施方式

下面将参考若干示例性实施方式来描述本发明的原理和精神。应当理解，给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。

本领域技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。

根据本发明的实施方式，提出了一种地质时期断层封闭性评价方法及系统。

下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。

图1是本发明一实施例的地质时期断层封闭性评价方法流程示意图。如图1所示，该方法包括：

步骤S1，采集石油地震勘探数据；

步骤S2，根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；

步骤S3，根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；

步骤S4，根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；

步骤S5，根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；

步骤S6，根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；

步骤S7，根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数P_f；

步骤S8，统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性P_f值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。

为了对上述地质时期断层封闭性评价方法进行更为清楚的解释，下面结合每一步骤进行详细说明。

步骤S1：

采集石油地震勘探数据。

步骤S2：

根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；

步骤S3：

根据石油地震勘探和钻井数据建立地层模型；

将所述地层模型与地下断层模型耦合，确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点。

步骤S4：

根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子。

参考图2，为本发明一具体实施例的断层倾角变化对断层倾角校正因子的影响关系示意图。

断面上每个点的断层倾角校正因子计算式为：

其中，f(θ)为断层倾角校正因子，无量纲；

θ为断层倾角值，°。

步骤S5：

根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值，包括：

其中，ESGR为有效断层岩的泥质含量，％；

i为滑过断点的第i层岩层；

T_i为滑过断点的第i层岩层厚度，m；

V_shi为滑过断点的第i层岩层泥质含量，％；

D为断层的垂直断距，m。

步骤S6：

根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率，包括：

根据一地质时期断层的断距，计算断层活动后的累加断距与对应活动时期的时间长度的比值，所述比值为一地质时期的断层活动速率；

归一化活动速率的计算式为：

其中，V_f为归一化活动速率，无量纲，％；

V为一地质时期的断层活动速率，m/Ma；

V_max和V_min分别为待评价断层在历史活动时期所具备的最大活动速率和最小的活动速率，m/Ma。

步骤S7：

根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数P_f，包括：

代表地质时期的断层封闭能力的参数P_f的计算式为：

其中，P_f为代表地质时期的断层封闭能力的参数。

步骤S8：

统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性Pf值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能，包括：

将P_f与封闭下限阈值比较；其中，若P_f比封闭下限阈值大，断层侧向封闭，反之断层侧向开启；

封闭能力通过P_f与封闭下限阈值的差值表征，若差值越大，断层封闭能力越强，反之越弱。

需要说明的是，尽管在上述实施例及附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。

为了对上述地质时期断层封闭性评价方法进行更为清楚的解释，下面结合一个具体的实施例来进行说明，然而值得注意的是该实施例仅是为了更好地说明本发明，并不构成对本发明不当的限定。

以苏丹Muglad盆地Kaikang凹陷西斜坡某断层为例，采用本发明提出的地质时期断层封闭性评价方法进行断层封闭性评价。

参考图3A，示例性展示了断层的归一化活动速率；图3B示例性展示了白垩纪末期油气关键成藏时刻断层的P_f分布情况；图3C示例性展示了古近纪末期时断层的P_fP_f的分布情况。

该断层的倾角在40-60度之间，对应的断层倾角校正因子f(θ)位于1.12-1.08；

另一参数SGR在下白垩统、上白垩统和古近系目的层段分别为24％～29％、28％～31％、34％～38％。

计算目的层段V_f所需的下白垩统、上白垩统和古近系所对应的地质历史时窗分别取值为17、20.6、44.8Ma，断距是根据构造图计算而来的。

该断层是北西-南东方向延展的，自北而南长约120km，在横向上具有分段生长的特征，见图3A，在平面上可分为三段：其中在第二分段内发现上白垩统的构造油藏；在第三分段中上白垩统圈闭钻探失利，在古近系中发现油气显示；第一分段尚未有钻井。

该断层地质时期历经3个演化时期：早白垩世、晚白垩世和古近纪活动时期。图3A表明整体上该断层在早白垩世末期和晚白垩世末期的活动速率普遍高于古近纪末期，这说明在该区的关键运移成藏期(晚白垩世末期)，该断层是活动较为强烈的，相应的封闭性能弱一些，有利于下白垩统源岩生成的油气沿着断面垂向运移到上白垩统的主要目的层；但是在横向分段的不同部位，其封闭能力有差异。

图3B中表明在关键油气成藏时期(晚白垩世末期)第二分段封闭性参数P_f整体小于第三分段，表明第二分段封闭性明显弱于第三分段，这说明该段是利于油气运移的，这与实际发现的油气聚集在第二分段的事实是符合的，相反，在油气运移的关键时期，第三分段封闭性高于第二分段，不利于当时油气的运移。

在图3C上整体表现为较强的封闭性，在第二分段P_f整体表现为高值，说明在古近纪该段断层具备良好了封闭性，这有利于油气的保存；但在第三分段P_f值呈现低值(根据第二分段油气分布情况，统计此处取P_f的阈值为5.2)，说明第三分段封闭性差，造成该段封闭性差的原因在于该段在古近纪的活动速率变大(见图3A)；由此可能会造成早期形成的油气藏沿着断层封闭性差的部位向上调整，最终致使油气散逸地表或者分散在古近系的圈闭中，这与在第三分段针对上白垩统断块圈闭目标钻探失利、在古近系中发现油气显示的事实是相吻合的，说明本发明的效果符合地质规律，具有可靠性。

在介绍了本发明示例性实施方式的方法之后，接下来，参考图4对本发明示例性实施方式的地质时期断层封闭性评价系统进行介绍。

地质时期断层封闭性评价系统的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的术语“模块”或者“单元”，可以是实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

基于同一发明构思，本发明还提出了一种地质时期断层封闭性评价系统，如图4所示，该系统包括：

数据采集模块410，用于采集石油地震勘探数据；

模型建立模块420，用于根据所述石油地震勘探数据，建立地下断层模型，并计算出断层的断距、断面上每个点的断层倾角值和垂直断距；

计算点选取模块430，用于根据地下断层模型确定待评价的目的层段在断面上的空间投影窗口位置，并选取断面上要进行封闭性评价的计算点；

校正因子计算模块440，用于根据断面上每个点的断层倾角值，计算得到断面上每个点的断层倾角校正因子；

泥质含量计算模块450，用于根据石油地震勘探数据及垂直断距，计算得到所述计算点的有效断层岩的泥质含量ESGR值；

活动速率计算模块460，用于根据断层的断距得到断层活动速率，并根据所述断层活动速率计算得到在一地质时期的断面上每个点的归一化活动速率；

封闭能力计算模块470，用于根据断面上每个点的断层倾角校正因子、计算点的有效断层岩的泥质含量ESGR值及在一地质时期的断面上每个点的归一化活动速率，计算得到代表地质时期的断层封闭能力的参数P_f；

封闭性评价模块480，用于统计勘探区域上已发现油藏断层封闭下限阈值，将待评价断层的地质时期封闭性P_f值与封闭下限阈值比较，根据比较结果定量评价断层在一地质时期的封闭性能。

应当注意，尽管在上文详细描述中提及了地质时期断层封闭性评价系统的若干模块，但是这种划分仅仅是示例性的并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多模块的特征和功能可以在一个模块中具体化。反之，上文描述的一个模块的特征和功能可以进一步划分为由多个模块来具体化。

基于前述发明构思，如图5所示，本发明还提出了一种计算机设备500，包括存储器510、处理器520及存储在存储器510上并可在处理器520上运行的计算机程序530，所述处理器520执行所述计算机程序530时实现前述地质时期断层封闭性评价方法。

基于前述发明构思，本发明提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现前述地质时期断层封闭性评价方法。

本发明提出的地质时期断层封闭性评价方法及系统引入断层活动速率这一与断层历史演化相关联的因素，反映了某地质时期的断层封闭特征；而且考虑了断层倾角变化的因素给断层岩泥质含量计算带来的影响，本发明所带来的评价效果更符合实际地质规律，有效提高了断层封闭性评价的准确度，降低断层圈闭钻探风险，为开发阶段优化断块油气藏开发方案提供地质评价参考。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。