[发明专利]一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法

说明书

本发明公开了一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法，主要步骤：S1、构建未加载闭合应力的天然裂缝初始缝宽；S2、建立碳酸盐岩裂缝受压变形模型，并计算闭合应力下的裂缝受压变形后的缝宽；S3、在步骤S2计算的加载闭合应力的缝宽基础上，建立裂缝酸化的流动反应模型，并计算该时间节点处的酸蚀缝宽；S4、在步骤S3计算得到的酸蚀缝宽的基础上，代入步骤S2，再进行步骤S3重新更新应力平衡后的酸蚀缝宽，得到下一个时间节点的酸蚀缝宽，循环迭代步骤S2和S3直至设定的注酸时间t，最终获取注酸结束时酸蚀缝宽分布。本发明方法能定量评价考虑闭合应力的裂缝岩样酸蚀缝宽，为超深井裂缝型碳酸盐岩储层酸化设计和效果评价提供指导。

技术领域

本发明涉及石油工程技术领域，特别涉及一种碳酸盐岩油气藏酸化增产过程中的考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法。

背景技术

近年来，在深层、超深层海相碳酸盐岩相继发现了多个裂缝型大气藏：如川西观雾山组气藏(平均埋深7500m)、下二叠统气藏(平均埋深7400m)、川东元坝长兴组气藏(平均埋深7000m)等，展现出巨大的勘探开发前景。这类气藏具有埋藏深(7000m)、天然裂缝发育、基质渗透率低(1mD)、破裂压力高(90MPa)等特点，受制于现有工程技术条件(井口装置承压受限)，无法实施大型酸压改造，因此酸化成为确保此类储层获得油气增产的关键技术。

深层、超深层裂缝型碳酸盐岩酸化过程中，受制于极高的施工压力，施工排量受限，导致酸液不足以完全开启天然裂缝，从而形成酸液在部分闭合天然裂缝中的流动反应模式。其酸化机理是：由于基质渗透率极低，天然裂缝将主导酸液流动，酸液进入天然裂缝增加缝内净压力，同时溶解、破坏原有接触支撑点，破坏缝内应力平衡，导致裂缝内的支撑应力重新分布；此外，在闭合应力作用下，酸液在天然裂缝内发生流动反应，非均匀刻蚀裂缝面岩石，最终形成具有一定油气流动能力的通道。详细论述见文献：郭建春,苟波,王坤杰,等.川西下二叠统超深气井网络裂缝酸化优化设计[J].天然气工业,2017,37(6):34-41.)，其中的记载。

而现有技术中忽略了“酸蚀裂缝受压变形”与“裂缝内酸流动反应”之间的动态耦合与相互竞争共同作用的影响，容易导致预测的酸蚀缝宽失真，影响深层、超深层海相裂缝型碳酸盐岩储的层酸化设计和效果评价。

发明内容

本发明的目的是提供一种考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法。

本发明提供的考虑闭合应力的裂缝岩样酸蚀缝宽数值计算方法，具体步骤如下：

S1、构建未加载闭合应力的天然裂缝初始缝宽。具体包括下列步骤：

S11、将采集自筛选目的层井下岩心或露头，沿天然裂缝延伸方向加工成标准API小岩柱或岩板。

S12、将岩柱或岩板放置于岩心剖缝器上，采用巴西劈裂方式沿着天然裂缝剖缝，获取两个天然裂缝面。

S13、获取两个裂缝面高程特征数据，再将获取的裂缝面高程数据的最小值作为基底位置，去掉裂缝面高程数据的基底，仅保留粗糙部分的特征，最后利用数据网格化软件将高程数据进行标准网格化。

S14、沿着裂缝面长度、宽度、高度方向建立三维直角坐标系，构成裂缝面形貌高程矩阵，接着，将建立的两个裂缝面形貌高程矩阵在三维空间内重新错位组合，建立起天然裂缝空间，得到每个网格单元内的初始缝宽b。

S2、考虑碳酸盐岩弹性-弹塑性变形特征，建立碳酸盐岩裂缝受压变形模型，并计算闭合应力下的裂缝受压变形后的缝宽。包括以下步骤：

S21、将缝宽分布数据存储为ASCII文件，利用程序读取并存储为缝宽分布矩阵b_i,j。

S22、输入包括裂缝面岩石的杨氏模量E、泊松比v、最大抗压强度σ_max在内的岩石力学参数以及法向加载的闭合应力σ_c。