CN115391739A - 一种裂缝渗透率的定量计算方法和系统 - Google Patents

Abstract

本发明涉及一种裂缝渗透率的定量计算方法和系统，包括以下步骤：基于获取的待分析油气田裂缝性储层各井的测井数据，确定各井中裂缝性储层的岩石矿物组分含量；基于待分析油气田裂缝性储层各井的岩石矿物组分含量以及斯通利波和电成像测井数据，计算得到待分析油气田中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；基于待分析油气田新钻井中的电成像测井数据，以及所述函数关系，计算得到该新钻井中裂缝性储层的裂缝渗透率。本发明可以广泛应用于油气勘探开发技术领域。

Description

一种裂缝渗透率的定量计算方法和系统

技术领域

本发明涉及一种裂缝渗透率的定量计算方法和系统，特别是关于一种裂缝性储层裂缝渗透率的定量计算方法和系统，属于油气勘探开发技术领域。

背景技术

向地球深部进军，拓展深层、超深层油气资源，对保障国家能源安全具有重要的现实与战略意义。目前，深层、超深层油气资源达671亿吨油当量，占中国资源总量的34％，也是未来油气勘探开发的现实领域。随着油气勘探开发不断向深层、超深层发展，裂缝性储层逐渐成为增储上产的主力，其中，裂缝是油气渗流的主要通道和产能贡献的重要因素。裂缝渗透率是评价裂缝有效性、预测裂缝性储层产能的关键参数，但裂缝渗透率的计算方法一直以来是业界研究的重点和难点。

裂缝渗透率的计算方法有两大类：第一类是根据裂缝宽度的大小，直接赋予裂缝渗透率一个定值，但该方法未能定量计算裂缝渗透率。第二类是通过岩石物理实验，通过带有裂缝的岩心，实验测量岩心的裂缝宽度和裂缝渗透率，并建立岩心的裂缝宽度与裂缝渗透率之间的函数关系；但该方法主要是根据实验手段获得，不同尺度的裂缝宽度与裂缝渗透率之间的函数关系不尽相同，尤其是裂缝渗透率计算结果的准确性未能得到实际资料的验证，从而给裂缝渗透率的计算结果带来较大的不确定性。

深层、超深层的作业成本非常高，裂缝性储层的准确评价，不仅可以大幅节约作业费用，而且可以优化井位部署和井位数量，从而大幅提高经济效益。由于没有明确裂缝渗透率的评价和检验方法，目前业界尚未形成裂缝渗透率的定量计算方法。

发明内容

针对上述问题，本发明的目的是提供一种裂缝渗透率的定量计算方法和系统，能够准确计算深层、超深层中裂缝性储层的裂缝渗透率，为评价裂缝有效性、预测裂缝性储层产能奠定基础，从而有效地提高深层、超深层裂缝性油气藏勘探开发的效益。

为实现上述目的，本发明采取以下技术方案：

第一方面，本发明提供一种裂缝渗透率的定量计算方法，其包括以下步骤：

基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及得到的函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

进一步，所述基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量的方法，包括以下步骤：

获取待分析油气田裂缝性储层在各井中预设深度下的岩心样品，开展岩心样品的物性分析实验，测量得到各岩心样品的孔隙度，作为各井的岩心孔隙度；

基于获取的待分析油气田裂缝性储层的一系列测井数据，计算得到与岩心孔隙度误差最小的测井孔隙度，并基于该测井孔隙度得到裂缝性储层在该井处的岩石矿物组分含量。

进一步，所述基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系的方法，包括以下步骤：

基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度；

基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

进一步，所述基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率的方法，包括以下步骤：

基于得到的待分析油气田裂缝性储层在各井处的岩石矿物组分含量，计算得到裂缝性储层在各井处的基质渗透率；

基于待分析油气田裂缝性储层在各井处的斯通利波数据，计算得到裂缝性储层在各井处的总渗透率；

基于待分析油气田裂缝性储层在各井处的总渗透率和基质渗透率，得到该油气田裂缝性储层在各井处的裂缝渗透率。

进一步，所述基质渗透率的计算公式为：

式中，K_int为基质渗透率；A_f为长石结构参数；

为孔隙度；B_i为第i种矿物的常数；M_i为第i种矿物的含量。

进一步，所述总渗透率的计算公式为：

K＝f(Δf,Δt)

式中，K为总渗透率，mD；Δf为斯通利波的频移，kHz；Δt为斯通利波的时滞，ms。

进一步，所述裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系为：

式中，K_f为裂缝渗透率，mD；

为裂缝孔隙度，无量纲；w_f为裂缝宽度，μm。

进一步，所述基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及所述函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率的方法，包括以下步骤：

基于待分析油气田裂缝性储层在新钻井处的电成像测井数据，处理解释得到裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度；

基于得到的裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度，以及所述函数关系，计算得到裂缝性储层在该新钻井处的裂缝渗透率。

第二方面，本发明提供一种裂缝渗透率的定量计算系统，包括：

组分含量确定模块，用于基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

函数关系确定模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及上述函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

进一步，所述函数关系确定模块，包括：

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

电成像数据处理模块，用于基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度；

参数拟合模块，用于基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

本发明由于采取以上技术方案，其具有以下优点：本发明通过对已有测井数据和岩心物理实验数据进行结合，得到油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度的函数关系，在对新钻井处裂缝渗透率进行计算时，可以避免开展大量的岩石物理实验和采集一系列高端的测井系列数据，能够有效地节约成本，提高裂缝性油气藏勘探开发的效益，具有较强的经济性，同时提高了计算深层、超深层的裂缝性储层裂缝渗透率的准确性。因此，本发明可以广泛应用于油气勘探开发技术领域。

附图说明

通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。在整个附图中，用相同的附图标记表示相同的部件。在附图中：

图1是本发明实施例提供的裂缝渗透率的定量计算方法流程图；

图2是Z气田裂缝性储层岩心铸体薄片显示所发育的裂缝；

图3是Z气田裂缝性储层电成像测井显示所发育的裂缝；

图4是Z气田A井裂缝性储层壁心岩心分析的孔隙度、基质渗透率与测井计算的孔隙度、基质渗透率；

图5是Z气田B井裂缝性储层全直径岩心分析的孔隙度、总渗透率与测井计算的孔隙度、总渗透率；

图6是Z气田不同井的裂缝孔隙度、裂缝宽度与裂缝渗透率之间的关系；

图7是Z气田F井的裂缝孔隙度、裂缝宽度与裂缝渗透率之间的关系。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例的附图，对本发明实施例的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于所描述的本发明的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。

需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。

本发明的一些实施例中，提供一种裂缝渗透率的定量计算方法，首先通过对某油气田裂缝性储层的多口井中的测井数据进行分析，得到该油气田中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；然后利用该函数关系，以及新钻井的电成像测井数据，即可得到该新钻井的裂缝渗透率。该方法，避免开展大量的岩石物理实验和采集一系列高端的测井系列数据，能够有效地节约成本，提高裂缝性油气藏勘探开发的效益，具有较强的经济性。

与之相对应地，本发明的一些实施例中还提供一种裂缝渗透率的定量计算系统。

实施例1

如图1所示，本实施例提供一种裂缝渗透率的定量计算方法，包括以下步骤：

1)基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

2)基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

3)基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及步骤2)中得到的函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

优选地，上述步骤1)包括以下步骤：

1.1)获取待分析油气田裂缝性储层在各井中预设深度下的岩心样品，开展岩心样品的物性分析实验，测量得到各岩心样品的孔隙度，作为各井的岩心孔隙度；

1.2)基于获取的待分析油气田裂缝性储层的一系列测井数据，计算得到与岩心孔隙度误差最小的测井孔隙度，并基于该测井孔隙度得到裂缝性储层在该井处的岩石矿物组分含量。

实际上，测井孔隙度计算的方法有多种，本实施例选择测井孔隙度的误差最小时对应的岩石矿物组分含量作为裂缝性储层的岩石矿物组分含量；一般情况下，测井孔隙度计算的误差需要小于8％。

优选地，上述步骤1.2)中，获取的待分析油气田裂缝性储层的一系列测井数据包括自然伽马、自然电位、深电阻率、中电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差、横波时差、斯通利波时差、电成像数据等。

其中，计算裂缝性储层的岩石矿物组分含量和测井孔隙度时，基于自然伽马、深电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差等测井数据计算，计算方法为本领域技术人员公知技术，本发明在此不再赘述。

优选地，上述步骤2)中，基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系的方法，包括以下步骤：

2.1)基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

2.2)基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度等参数；

2.3)基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

优选地，上述步骤2.1)中，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率的方法，包括以下步骤：

2.1.1)基于得到的待分析油气田裂缝性储层在各井处的岩石矿物组分含量，计算得到裂缝性储层在各井处的基质渗透率。

具体地，本实施例中采用Herron模型计算裂缝性储层的基质渗透率，计算公式为：

式中，K_int为基质渗透率，mD；A_f为长石结构参数，无量纲；

为孔隙度，小数；B_i为第i种矿物的常数，无量纲；M_i为第i种矿物的含量，小数。

2.1.2)基于待分析油气田裂缝性储层在各井处的斯通利波数据，计算得到裂缝性储层在各井处的总渗透率。

本实施例中，将斯通利波的主频和到时进行匹配，采用最小二乘法构建裂缝性储层总渗透率的计算公式：

K＝f(Δf,Δt) (2)

式中，K为总渗透率，mD；Δf为斯通利波的频移，kHz；Δt为斯通利波的时滞，ms。

2.1.3)基于待分析油气田裂缝性储层在各井处的总渗透率和基质渗透率，得到该油气田裂缝性储层在各井处的裂缝渗透率，计算公式为：

K_f＝K-K_int (3)

式中，K_f为裂缝渗透率，mD；K为总渗透率，mD；K_int为基质渗透率，mD。

优选地，上述步骤2.3)中，建立的裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系为：

式中，K_f为裂缝渗透率，mD；

为裂缝孔隙度，无量纲；w_f为裂缝宽度，μm。其中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系可以根据数值的分布情况，采用指数或对数等函数形式进行拟合。

优选地，上述步骤3)中，基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及步骤2)中得到的函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率的方法，包括以下步骤：

3.1)基于待分析油气田裂缝性储层在新钻井处的电成像测井数据，处理解释得到裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度；

3.2)基于得到的裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度，以及步骤2)中的函数关系(即公式4)，计算得到裂缝性储层在该新钻井处的裂缝渗透率。

实施例2

本实施例以Z气田裂缝性储层为例，结合附图详细阐述本发明实施例1中的裂缝渗透率的定量计算方法的具体过程。

1)Z气田裂缝性储层岩心铸体薄片和电成像测井显示所发育的裂缝分别见图2和图3；基于Z气田裂缝性储层一定深度下的岩心样品，开展岩心样品的物性分析实验，从而可以测量得到每一块岩心样品的孔隙度，即岩心孔隙度。

2)针对Z气田的裂缝性储层，测井仪器可以测量得到一系列的测井数据；基于自然伽马、深电阻率、浅电阻率、体积密度、中子孔隙度、纵波时差等测井数据，可以计算得到裂缝性储层的岩石矿物组分含量和测井孔隙度，并使测井计算的孔隙度与步骤1)实验测量得到的岩心孔隙度之间的误差最小。

Z气田A井裂缝性储层壁心岩心分析的孔隙度与测井计算的孔隙度见图4中的第5道，B井裂缝性储层全直径岩心分析的孔隙度与测井计算的孔隙度见图5中的第5道，岩心分析的孔隙度与测井计算的孔隙度吻合程度均较高，绝对误差仅为0.42％。

3)基于步骤2)所得到Z气田裂缝性储层的岩石矿物组分含量，主要为石英、长石和云母，采用式1中的Herron模型可以计算得到裂缝性储层的基质渗透率。

Z气田A井裂缝性储层壁心渗透率的分析结果即为基质渗透率。如图4所示，在Z气田A井裂缝性储层孔隙度和基质渗透率定量计算成果图中，第1道是地层测量深度；第2道是自然伽马，表示该地层的岩性特征；第3道是深、浅电阻率测井曲线，刻画该地层的电性特征；第4道是体积密度、中子孔隙度和声波时差，反映该地层的物性特征；第5道是壁心岩心分析的孔隙度与测井计算的孔隙度；第6道是裂缝性储层壁心分析的基质渗透率与测井计算的基质渗透率，二者的吻合程度较高，相对误差为28.5％，验证了测井计算基质渗透率计算结果的合理性。

4)基于Z气田一口井中裂缝性储层的斯通利波，计算该井中裂缝性储层总渗透率。

Z气田全直径岩心渗透率的分析结果即为总渗透率，图5中的第6道是Z气田B井裂缝性储层全直径岩心分析的总渗透率与测井计算的总渗透率，二者的吻合程度较高，相对误差为15.56％，验证了测井计算的总渗透率计算结果的合理性。

5)基于Z气田一口井中裂缝性储层的总渗透率，与步骤3)所计算裂缝性储层的基质渗透率，得到该油气田一口井中裂缝性储层的裂缝渗透率。

Z气田A井、B井、C井、D井、E井5口井裂缝性储层的孔隙度、基质渗透率、总渗透率和裂缝渗透率见表1。

表1 Z气田5口井裂缝性储层测井计算结果

6)基于Z气田一口井中裂缝性储层的电成像测井数据，处理、解释得到裂缝倾角、裂缝密度、裂缝孔隙度、裂缝宽度等参数。

7)重复上述步骤3)到步骤6)，从而可以得到该油气田多口井中裂缝性储层的裂缝渗透率、裂缝孔隙度和裂缝宽度，见表2。

表2 Z气田5口井裂缝性储层测井计算的裂缝参数

井号	裂缝渗透率/mD	裂缝孔隙度	裂缝宽度/μm
				A	2.3	0.5％	189.0
B	2.0	0.6％	168.0
				C	1.2	0.5％	60.8
D	0.5	0.2％	34.9
				E	0.2	0.2％	57.7

8)基于得到的该油气田多口井中裂缝性储层的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系，如图6所示。

9)当Z气田新钻井F中未采集斯通利波，但采集电成像测井时，采用步骤6)计算该新钻井的裂缝孔隙度为0.30％、裂缝宽度为62.4μm，基于步骤8)中的式(4)可以计算这口新钻井中裂缝性储层的裂缝渗透率为0.72，如图7所示。

实施例3

上述实施例1提供了一种裂缝渗透率的定量计算方法，与之相对应地，本实施例提供一种裂缝渗透率的定量计算系统。本实施例提供的系统可以实施实施例1的一种裂缝渗透率的定量计算方法，该系统可以通过软件、硬件或软硬结合的方式来实现。例如，该系统可以包括集成的或分开的功能模块或功能单元来执行实施例1各方法中的对应步骤。由于本实施例的系统基本相似于方法实施例，所以本实施例描述过程比较简单，相关之处可以参见实施例1的部分说明即可，本实施例提供的系统的实施例仅仅是示意性的。

本实施例提供的一种裂缝渗透率的定量计算系统，包括：

组分含量确定模块，用于基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

函数关系确定模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及上述函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

优选地，函数关系确定模块，包括：

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

电成像数据处理模块，用于基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度等参数；

参数拟合模块，用于基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求保护范围之内。

上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (10)

1.一种裂缝渗透率的定量计算方法，其特征在于，包括以下步骤：

基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及得到的函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

2.如权利要求1所述的一种裂缝渗透率的定量计算方法，其特征在于，所述基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量的方法，包括以下步骤：

获取待分析油气田裂缝性储层在各井中预设深度下的岩心样品，开展岩心样品的物性分析实验，测量得到各岩心样品的孔隙度，作为各井的岩心孔隙度；

基于获取的待分析油气田裂缝性储层的一系列测井数据，计算得到与岩心孔隙度误差最小的测井孔隙度，并基于该测井孔隙度得到裂缝性储层在该井处的岩石矿物组分含量。

3.如权利要求1所述的一种裂缝渗透率的定量计算方法，其特征在于，所述基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系的方法，包括以下步骤：

基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度；

基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

4.如权利要求3所述的一种裂缝渗透率的定量计算方法，其特征在于，所述基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率的方法，包括以下步骤：

基于得到的待分析油气田裂缝性储层在各井处的岩石矿物组分含量，计算得到裂缝性储层在各井处的基质渗透率；

基于待分析油气田裂缝性储层在各井处的斯通利波数据，计算得到裂缝性储层在各井处的总渗透率；

基于待分析油气田裂缝性储层在各井处的总渗透率和基质渗透率，得到该油气田裂缝性储层在各井处的裂缝渗透率。

5.如权利要求4所述的一种裂缝渗透率的定量计算方法，其特征在于，所述基质渗透率的计算公式为：

式中，K_int为基质渗透率，mD；A_f为长石结构参数，无量纲；

为孔隙度，小数；B_i为第i种矿物的常数，无量纲；M_i为第i种矿物的含量，小数。

6.如权利要求4所述的一种裂缝渗透率的定量计算方法，其特征在于，所述总渗透率的计算公式为：

K＝f(Δf,Δt)

式中，K为总渗透率，mD；Δf为斯通利波的频移，kHz；Δt为斯通利波的时滞，ms。

7.如权利要求4所述的一种裂缝渗透率的定量计算方法，其特征在于，所述裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系为：

式中，K_f为裂缝渗透率，mD；

为裂缝孔隙度，无量纲；w_f为裂缝宽度，μm。

8.如权利要求1所述的一种裂缝渗透率的定量计算方法，其特征在于，所述基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及所述函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率的方法，包括以下步骤：

基于待分析油气田裂缝性储层在新钻井处的电成像测井数据，处理解释得到裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度；

基于得到的裂缝性储层在新钻井处的裂缝孔隙度和裂缝宽度，以及所述函数关系，计算得到裂缝性储层在该新钻井处的裂缝渗透率。

9.一种裂缝渗透率的定量计算系统，其特征在于，包括：

组分含量确定模块，用于基于获取的待分析油气田裂缝性储层在各井处的测井数据和岩心样品物性分析实验数据，确定待分析油气田裂缝性储层在各井处的岩石矿物组分含量；

函数关系确定模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及对应的斯通利波和电成像测井数据，计算得到待分析油气田裂缝性储层中裂缝渗透率与裂缝孔隙度和裂缝宽度之间的函数关系；

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在新钻井中的电成像测井数据，以及上述函数关系，计算得到待分析油气田裂缝性储层在新钻井处的裂缝渗透率。

10.如权利要求9所述的一种裂缝渗透率的定量计算系统，其特征在于，所述函数关系确定模块，包括：

裂缝渗透率计算模块，用于基于待分析油气田裂缝性储层在各井处的岩石矿物组分含量以及斯通利波，计算得到待分析油气田裂缝性储层在各井处的裂缝渗透率；

电成像数据处理模块，用于基于待分析油气田裂缝性储层在各井处的电成像测井数据，处理、解释得到裂缝孔隙度和裂缝宽度；

参数拟合模块，用于基于得到的裂缝性储层在各井处的裂缝渗透率、裂缝孔隙度和裂缝宽度，建立待分析油气田裂缝性储层中，裂缝渗透率与裂缝孔隙度、裂缝宽度之间的函数关系。

Images