确定致密砂岩微裂缝发育程度的方法、装置和存储介质与流程

1.本发明涉及勘探测井技术领域，尤其涉及一种确定致密砂岩微裂缝发育程度的方法、装置、存储介质和计算机设备。


背景技术：

2.针对裂缝的识别研究，测井技术通常是利用地层微电阻率扫描测井和井下声波电视成像测井，根据裂缝在声、电成像图上的特征识别出裂缝，并计算出裂缝的一系列参数。同时，根据声、电成像图上识别出的裂缝，研究裂缝在常规测井曲线上的测井响应特征，进而利用常规测井曲线识别裂缝。但有些裂缝非常微小，在声、电成像图上反映不出来，仅可从薄片上窥见其真容，此时无法利用成像资料对常规测井资料进行微裂缝标定识别。
3.目前，对裂缝的研究很多，包括裂缝的测井响应特征、裂缝的表征、参数计算及识别方法研究等，但对微裂缝的研究相对较少。有部分研究者对页岩地层微裂缝进行了研究。如有研究者选取声波时差、自然伽马和补偿中子等三条对储层微裂缝敏感度较高的常规测井曲线，运用人工神经网络和小波变换对下寺湾油气田长7段页岩储层的微裂缝进行了识别；还有的研究者通过计算实测地层岩石体积压缩系数clog与基于页岩矿物成分计算的理论体积压缩系数cth之间的差异、原状地层电阻率rt与地层冲洗带电阻率rxo之间的差异，得到页岩气裂缝发育指数fi，以此实现对裂缝发育段的识别。有少数研究者对非页岩储层微裂缝进行了识别研究，研究者根据岩石基质的体积模量、剪切模量和双孔隙结构参数，建立干岩石骨架弹性模量比随孔隙度变化的关系式以及双孔dem解析模型；配置多组软硬孔的孔隙体积百分比，根据上述解析模型，得到多个干岩石骨架的体积模量和剪切模量；根据多个干岩石骨架的体积模量和剪切模量，应用gassmann方程得到多个饱和岩石的体积模量和剪切模量；分别计算每个应用gassmann方程得到的饱和岩石的体积模量和剪切模量与根据实际测井数据得到的饱和岩石的体积模量和剪切模量之间的误差，当误差最小时，得到最优的软孔孔隙体积百分比；根据最优的软孔孔隙体积百分比，识别储层微裂缝。目前对微裂缝的研究均侧重于某个具体微裂缝指示参数的计算，其必须借助于特定的实验测量参数或特定的测井曲线，不能直观、方便地利用常规测井曲线对微裂缝进行识别及确定微裂缝发育程度，应用范围有限。目前利用常规测井曲线识别与确定微裂缝发育程度的研究成果未见之于任何专利及文献。
4.致密砂岩储层中的微裂缝与储层射孔及压裂改造效果关系密切，微裂缝发育与否与储层是否具备经济产能息息相关，它与储层甜点的分布关系在专利及文献中鲜有讨论，因此有必要深化研究。


技术实现要素：

5.针对上述问题，本发明提供了一种基于测井曲线确定致密砂岩微裂缝发育程度的方法，包括以下步骤：
6.分析致密砂岩样本的微裂缝面孔率，用于表征致密砂岩的微裂缝地层的发育程
度；
7.根据致密砂岩样本的深电阻率和冲洗带电阻率，确定致密砂岩的微裂缝地层的电阻率差异度；
8.基于测井曲线，利用测井声波和测井密度计算微裂缝地层的波阻抗；
9.统计分析微裂缝地层的发育程度、电阻率差异度和波阻抗之间的关系,建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式；
10.利用所述计算式，根据致密砂岩待测微裂缝地层的电阻率差异度和波阻抗估算该微裂缝地层的微裂缝发育程度。
11.在本发明的一个实施例中，分析致密砂岩样本的微裂缝面孔率，包括以下步骤：
12.获取致密砂岩岩心薄片,作为致密砂岩样本；
13.对岩心分析数据进行归位深度校正，并利用岩心深度校正量对薄片深度进行校正；
14.利用深度校正后的薄片数据确定微裂缝面孔率。
15.在本发明的一个实施例中，根据下式确定致密砂岩的微裂缝地层的电阻率差异度：
16.rd＝lg(rt/rxo)
17.式中，rd为微裂缝地层的电阻率差异度，rt为微裂缝地层的深电阻率；rxo为微裂缝地层的冲洗带电阻率。
18.在本发明的一个实施例中，根据下式计算微裂缝地层的波阻抗：
19.ai＝100*den/ac
20.式中，ai为微裂缝地层的波阻抗，ac为测井声波时差，den为测井密度。
21.在本发明的一个实施例中，统计分析微裂缝地层的发育程度、电阻率差异度和波阻抗之间的关系,建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式，包括：
22.通过统计分析确定微裂缝地层的发育程度与电阻率差异度之间的第一关系；
23.通过统计分析确定微裂缝地层的发育程度与波阻抗之间的第二关系；
24.利用多元回归法，根据所述第一和第二关系建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式。
25.在本发明的一个实施例中，所述第一关系为正相关系，所述第二关系为负相关系。
26.在本发明的一个实施例中，所述计算式为：
27.mfdd＝a*ai*rd b*rd c*rd/ai d
28.式中，mfdd为微裂缝发育程度，rd为电阻率差异度、ai为波阻抗，a、b、c、d为系数。
29.此外，本发明还提供一种确定致密砂岩微裂缝发育程度的装置，其特征在于，包括：
30.发育分析模块，用于分析致密砂岩样本的微裂缝面孔率，用于表征致密砂岩的微裂缝地层的发育程度；
31.阻率分析模块，用于根据致密砂岩样本的深电阻率和冲洗带电阻率，确定致密砂岩的微裂缝地层的电阻率差异度；
32.阻抗分析模块，用于基于测井曲线，利用测井声波和测井密度计算微裂缝地层的
波阻抗；
33.统计分析模块，用于统计分析微裂缝地层的发育程度、电阻率差异度和波阻抗之间的关系，建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式；
34.发育估算模块，用于利用所述计算式，根据致密砂岩待测微裂缝地层的电阻率差异度和波阻抗估算该微裂缝地层的微裂缝发育程度。
35.另外，本发明还提供一种存储介质，其中存储有计算机程序，其特征在于，所述计算机程序被处理器执行时，实现上述的确定致密砂岩微裂缝发育程度的方法的步骤。
36.本发明还提供一种计算机设备，包括存储器和处理器，其中所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，实现上述的确定致密砂岩微裂缝发育程度的方法的步骤。
37.与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
38.目前研究者们对声、电成像图上可识别的明显裂缝研究较多，对声、电成像图上无法反映出来的微裂缝研究较少。有部分研究者对页岩地层微裂缝进行了研究，且研究方法需要借助特定的实验测量参数或特定的测井曲线，应用范围有限。对砂岩微裂缝的研究相对较少，关于砂岩微裂缝发育程度的研究更是少人问津。微裂缝发育程度对致密砂岩储层的改造、产能的提高都有着重要作用。本发明提供的利用测井曲线计算砂岩储层微裂缝发育程度的方法，在没有薄片分析情况下的，利用常规测井曲线即可确定微裂缝地层微裂缝发育程度，为储层工程工艺的选择提供较为准确的依据。该方法易于实施，可操作性强，在现场实际中应用范围很广，应用前景十分广阔。
39.本发明的其他优点、目标，和特征在某种程度上将在随后的说明书中进行阐述，并且在某种程度上，基于对下文的考察研究对本领域技术人员而言将是显而易见的，或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书，权利要求书，以及附图中所特别指出的结构来实现和获得。
附图说明
40.附图用来提供对本技术的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本技术实施例的附图与本技术的实施例一起用于解释本技术的技术方案，但并不构成对本技术技术方案的限制。
41.图1是根据本发明实施例二的确定致密砂岩微裂缝发育程度的方法的流程示意图；
42.图2是根据本发明实施例三的确定致密砂岩微裂缝发育程度的方法的流程示意图；
43.图3是根据本发明实施例四的微裂缝发育程度(mfdd)与电阻率差异度(rd)的关系示意图；
44.图4为根据本发明实施例四的微裂缝发育程度(mfdd)与波阻抗(ai)的关系图的示意图；
45.图5是根据本发明实施例四的方法计算微裂缝发育程度与实测微裂缝发育程度的对比表。
具体实施方式
46.针对目前因缺少直观、方便、有效方法判断致密砂岩储层微裂缝发育程度而成为制约判断储层是否具有甜点的瓶颈因素，导致储层压裂施工成功率不高或识别不出有经济产能的油气层的问题，本发明提出一种利用测井曲线确定致密砂岩储层微裂缝发育程度的方法。
47.具体地，由于微裂缝增加了地层的渗透性，令声波值升高，同时也使得地层易扩径，造成冲洗带电阻率与密度曲线很大程度上反映泥浆滤液信息，测量值明显降低。电阻率差异度与地层渗透性相关，为微裂缝发育程度的反映；波阻抗源自声波和密度，声波、密度均受微裂缝影响，故而波阻抗也与微裂缝发育程度有关。该方法通过分析微裂缝发育程度与电阻率差异度、波阻抗的关系，建立了微裂缝发育程度计算公式,据此即可估算出微裂缝地层微裂缝发育程度。
48.以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成相应技术效果的实现过程能充分理解并据以实施。本技术实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本发明的保护范围之内。
49.实施例一
50.本实施例提供一种确定致密砂岩微裂缝发育程度的装置，其包括：
51.发育分析模块，用于分析致密砂岩样本的微裂缝面孔率，用于表征致密砂岩的微裂缝地层的发育程度；
52.阻率分析模块，用于根据致密砂岩样本的深电阻率和冲洗带电阻率，确定致密砂岩的微裂缝地层的电阻率差异度；
53.阻抗分析模块，用于基于测井曲线，利用测井声波和测井密度计算微裂缝地层的波阻抗；
54.统计分析模块，用于统计分析微裂缝地层的发育程度、电阻率差异度和波阻抗之间的关系，建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式；
55.发育估算模块，用于利用所述计算式，根据致密砂岩待测微裂缝地层的电阻率差异度和波阻抗估算该微裂缝地层的微裂缝发育程度。
56.实施例二
57.如图1所示，本实施例提供一种确定致密砂岩微裂缝发育程度的方法，包括以下步骤：
58.s10，分析致密砂岩样本的微裂缝面孔率，用于表征致密砂岩的微裂缝地层的发育程度；
59.s20，根据致密砂岩样本的深电阻率和冲洗带电阻率，根据下式确定致密砂岩的微裂缝地层的电阻率差异度：
60.rd＝lg(rt/rxo)
61.式中，rd为微裂缝地层的电阻率差异度，rt为微裂缝地层的深电阻率；rxo为微裂缝地层的冲洗带电阻率。
62.s30，基于测井曲线，利用测井声波和测井密度，根据下式计算微裂缝地层的波阻抗：
63.ai＝100*den/ac
64.式中，ai为微裂缝地层的波阻抗，ac为测井声波时差，den为测井密度。
65.s40，统计分析微裂缝地层的发育程度、电阻率差异度和波阻抗之间的关系,建立用于根据电阻率差异度和波阻抗确定微裂缝发育程度的计算式：
66.mfdd＝a*ai*rd b*rd c*rd/ai d
67.式中，mfdd为微裂缝发育程度，rd为电阻率差异度、ai为波阻抗，a、b、c、d为系数。
68.s50，利用所述计算式，根据致密砂岩待测微裂缝地层的电阻率差异度和波阻抗估算该微裂缝地层的微裂缝发育程度。
69.实施例三
70.如图2所示，本实施例优选在薄片资料深度校正的基础上，以薄片分析微裂缝面孔率表征微裂缝发育程度，通过研究微裂缝发育程度与电阻率差异度、波阻抗的关系，建立了微裂缝发育程度计算式，据此即可估算出微裂缝发育程度，为生产施工提供依据。其具体做法是：
71.(1)薄片深度校正：首先对常规岩心分析数据进行归位深度校正，然后依据岩心深度校正量对薄片深度进行校正。
72.(2)依据薄片分析的微裂缝面孔率来表征微裂缝发育程度(mfdd)。
73.(3)微裂缝地层电阻率差异度(rd)计算：根据深电阻率(rt)及冲洗带电阻率(rxo)，计算微裂缝地层之电阻率差异度，即rd＝lg(rt/rxo)。
74.(4)微裂缝地层波阻抗(ai)计算：根据测井声波时差(ac)及测井密度(den)，计算微裂缝地层波阻抗，即ai＝100*den/ac。
75.(5)微裂缝发育程度(mfdd)计算式确定：统计分析微裂缝发育程度与电阻率差异度(rd)、波阻抗(ai)之间的关系，通过拟合分析建立微裂缝发育程度计算式，即mfdd＝a*ai*rd b*rd c*rd/ai d
76.(6)根据(5)确定的计算式，即可估算出微裂缝发育程度。
77.实施例四
78.下面结合上述方法在某地某气田微裂缝地层的应用对本发明作进一步详细描述：
79.具体地，应用过程包括如下步骤:
80.(1)薄片深度校正。
81.薄片源自岩心，而岩心深度与测井深度常常不一致，在使用薄片资料前，需先对其进行深度校正。首先对常规岩心分析数据进行归位深度校正，然后依据岩心深度校正量对薄片深度进行校正。
82.(2)依据薄片分析的微裂缝面孔率来表征微裂缝发育程度(mfdd)。
83.根据薄片分析中微裂缝面孔率的大小来表征微裂缝的发育程度，微裂缝面孔率越大，微裂缝越发育。
84.(3)根据深电阻率(rt)及冲洗带电阻率(rxo)，计算微裂缝地层之电阻率差异度(rd)。
85.rd＝lg(rt/rxo)
86.深电阻率(rt)与冲洗带电阻率(rxo)比值的对数值即反映了两者之间的曲线幅度差，其代表了地层渗透性的强弱。由于微裂缝的存在，增加了地层的渗透性。
87.(4)根据声波(ac)及密度(den)，计算微裂缝地层波阻抗(ai).
88.ai＝100*den/ac。
89.随着微裂缝发育程度的增加，声波值升高，密度值降低，导致地层波阻抗降低。
90.(5)分别统计分析微裂缝发育程度与电阻率差异度、波阻抗之间的关系，再利用多元回归建立微裂缝发育程度计算式。
91.电阻率差异度、波阻抗都与微裂缝地层发育程度相关，均可反映微裂缝发育程度。图3为微裂缝发育程度(mfdd)与电阻率差异度(rd)的关系图。图4为微裂缝发育程度(mfdd)与波阻抗(ai)的关系图。由图3和图4可知，微裂缝发育程度与电阻率差异度为正相关关系，与波阻抗为负相关关系，且均具有一定的相关度。因此，对微裂缝发育程度(mfdd)与电阻率差异度(rd)、波阻抗(ai)进行了统计回归，得到了微裂缝发育程度(mfdd)的计算公式：
92.mfdd＝a*ai*rd b*rd c*rd/ai d
93.(6)利用(5)中的计算式，对微裂缝地层的微裂缝发育程度进行计算，并将其与实际分析数据进行对比。
94.图5为微裂缝地层实测微裂缝发育程度与通过本发明提供的方法计算微裂缝发育程度的对比表。由图中数据可以看出，微裂缝发育程度计算结果与实测结果很接近。
95.本发明提供的微裂缝发育程度计算方法，易于实施，可操作性强。该方法在某地某气田微裂缝地层应用效果好，符合率达到75％，估算的微裂缝发育程度平均绝对误差仅为0.019。
96.实施例四
97.此外，为解决现有技术中存在的上述技术问题，本发明实施例还相应地提供一种存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述的确定致密砂岩微裂缝发育程度的方法的步骤。
98.实施例五
99.此外，为解决现有技术中存在的上述技术问题，本发明实施例还相应地提供一种计算机设备，包括存储器和处理器，其中所述存储器存储有计算机程序，所述计算机程序被所述处理器执行时，实现上述的确定致密砂岩微裂缝发育程度的方法的步骤。
100.本发明的目的是克服现有的确定致密砂岩微裂缝发育程度的方法的处理繁琐、分析复杂、成本较高，周期过长等不足，采用完全不同于现有的确定致密砂岩微裂缝发育程度的方法的技术思路，提供一种对微裂缝发育程度进行快速有效估算的简易方法，从而满足油气勘探开发的现场实际需求。
101.需要说明的是，本发明实施例的方法可以由单个设备执行，例如一台计算机或服务器等。本实施例的方法也可以应用于分布式场景下，由多台设备相互配合来完成。在这种分布式场景的情况下，这多台设备中的一台设备可以只执行本发明实施例的方法中的某一个或多个步骤，这多台设备相互之间会进行交互以完成所述的方法。
102.流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本发明的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本发明的实施例所属技术领域的技术人员所理解。
103.应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
104.本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤是可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
105.此外，在本发明各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。所述集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
106.上述提到的存储介质可以是只读存储器，磁盘或光盘等。
107.虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的保护范围，仍须以所附的权利要求书所界定的范围为准。